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J, 

Predictive  heat  transfer  equations  have  been  developed  to  estimate  the 

i 

local  supplementary  heating  necessary  to  support  a  resting  diver.  These 

formulations  consider  the  diver's  passive  thermal  protection  in  the  form  of  a 

j 

i  diving  dress  and  the  effects  of  his  ambient  environment  upon  the  dress' 

I 

thermal  resistance.  The  mathematical  model  performs  a  numerical  analysis  on 
discrete  composite  layers  comprised  of  the  skin,  dry  protective  suit, 
breathing  gas  mixture,  and  the  surrounding  fluid  environment.  The  parametric 
properties  associated  with  the  heat  transfer  correlations  were  obtained  from 
validated  experimental  heat  flux  and  skin  temperature  data  collected  under  j 
specific  environmental  conditions.  j 

i 

The  regional  skin  temperatures  of  this  model  are  predicted  by  using  a 
series  of  linear  algebraic  equations  which  have  been  derived  from  those  of 
Kerslake  and  which  predict  the  segmental  skin  temperatures  of  a  i 

quasi-euthermlc  diver  haying  pronounced  vasoconstriction.  The  equations  were 
developed  by  numerically  analyzing  twelve  segmental  skin  temperatures  recorded  j 

f 

from  five  resting  subjects  exposed  to  a  20  °C,  95%  helium  -  5%  oxygen 
j  hyperbaric  chamber  environment  pressurized  to  the  equivalent  of  200  meters  of 
J  seawater.  The  independent  variable  of  the  basic  correlation  is  an  arbitrarily  " 


defined  mean  skin  temperature.  The  developed  equations  were  authenticated  by 
comparing  the  predicted  normalized  segmental  temperature  with  the  respective 
experimental  normalized  temperature  obtained  from  several  sets  of 
physiological  data  which  were  collected  during  the  evaluation  of  the  Naval 
Coastal  System  Center's  prototype  Diver  Thermal  Protection  garment.  The 
results  of  these  comparisons  of  non-dimensional  1  zed  temperatures  Indicated 
that  the  derived  correlations  should  accurately  predict  the  skin  temperature 
of  the  principal  segments  as  a  function  of  mean  skin  temperature  with  a 
nominal  error  of  no  more  than  15X. 

The  supplementary  heating  that  is  required  to  maintain  a  diver  Is 
obtained  by  determining  the  heat  flux  that  would  be  lost  were  he  not  heated 
and  subtracting  from  this  the  normal  heat  flux  expected  from  a  resting 
euthermic  subject.  Applying  this  principle  to  each  of  the  regions.  It  Is 
possible  to  predict  the  supplementary  heating  necessary  to  compenstate  for  an 
excessive  flux. 

vThe  mathematical  model  was  written  as  an  Interactive  computer  program. 
The  program  collects  information  on  the  diver's  posture  and  the  ambient 
environment,  and  then  generates  estimates  of  the  required  regional 
supplementary  heating.  The  program  also  allows  the  comparison  of  predicted 
temperatures,  heat  fluxes,  overall  heat  transfer  coefficients,  normalized 
temperatures,  and  Blot  numbers  to  values  computed  from  experimental  data. 
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Portions  of  this  report  have  been  excerpted  from  the  original  ONR  Contract 
proposal  of  Dr.  Charles  E.  Johnson  and  that  contract's  annual  reports  for  1980 
and  1981.  Host  of  Chapters  III-V  and  part  of  Chapter  VII  have  been  presented 
at  scientific  meetings  of  the  Undersea  Medical  Society  and  the  American  Society 
of  Mechanical  Engineers.  The  predictive  skin  temperature  model  (Chapter  IV) 
was  published  in  Hay  of  1982  in  Undersea  Biomedical  Research  and  the  overall 
heat  loss  model  is  summarized  in  a  paper  accepted  for  publication  in  1982  by 
the  Journal  of  Biomechanical  Engineering,  a  transaction  of  the  American  Society 
of  Mechanical  Engineers.  This  document,  with  the  exception  of  the  Addendum: 
Mixed  Breathing  Gas  Interactive  Model  Program,  was  presented  as  the 
dissertation  of  James  D.  Collins  in  partial  fulfillment  of  the  requirements  for 
the  degree  of  Doctor  of  Philosophy  In  the  Department  of  Mechanical  Engineering 
and  Materials  Science  in  the  Graduate  School  of  Duke  University. 

Several  terms  which  are  used  throughout  this  document  Have  specially 
Intended  connotations.  These  terms  ire  defined  on  the  following  pages  for  the 
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benefit  of  the  reader. 


Segment 


Segment  or  segments  Is  used  when  describing  the  twelve  Hody  segments 
(Fig.  2).  Whenever  segment  appears  in  the  text,  It  refers  to  these  twelve 

j  Kody  segments. 

i 

|  Region 

I 

Region  is  used  to  describe  the  nine  model  man  regions  (Fig.  1)  throughout 
all  chapters,  except  Chapter  IV  where  it  Is  also  used  to  describe  Kerslake's 
eight  regions.  In  general,  any  references  to  regions(s)  refer  to  the  nine 
region  model  man  of  Figure  1. 


Rate  of  Heat  Loss/Flow 

Rate  of  heat  loss  or  flow  is  used  to  imply  watts  (joules/sec).  Any 
reference  to  a  rate  of  heat  loss  should  be  associated  with  the  units  of  watts. 


Heat  Flux 


Heat  flux  refers  to  the  rate  of  heat  loss  per  unit  of  area  (W/H^). 


Time  Record 


!  Time  record  is  used  to  refer  to  the  time  moment  associated  with  a  single 

I 

I  data  event  in  the  experimental  data  files.  Each  time  record  contains  the 
|  event  time,  twelve  segmental  rates  of  heat  loss  or  heat  fluxes,  twelve 

I 

1 

i  associated  temperatures,  and  the  ambient  and  rectal  temperatures.  (See 
j  Appendix  A  for  details.). 

I 

File 

File  refers  to  the  individual  experimental  data  files  of  Tables  2-4 
(Chapter  II).  Each  file  contains  multiple  time  records  and  represents  a 
cos^lete  experiment.  File  is  used  to  refer  to  both  the  individual  subject 
and  the  composite  subject  data  files. 

We  have  conducted  a  comparative  mathematical  investigation  of  limited 
experimental  data  of  skin  temperature  and  heat  flux  obtained  from  subjects  in 
j  various  environments  and  states  of  thermal  equilibrium.  Except  to  delete 
obviously  incongruous  data  points,  the  data  that  we  used  was  accepted  as  being 
correct.  The  reader  is  referred  to  the  Fiaval  Medical  Research  and  Development 
Command's  Scientific  Protocol  for  the  Conduct  of  Thermal  Physiology 
Experiments  on  Divers  for  a  discussion  of  the  potential  error  associated  with 
diver  thermal  physiology  experiments.  Our  investigation  resulted  in  the 
proposed  correlations  to  predict  supplemental  heating  requirements  to  maintain 
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a  resting  subject  In  a  state  of  euthermia  in  any  ambient  environment.  The 
work  is  based  upon  the  hypothesis  of  extrapolating  empirical  heat  transfer 
formuations  that  are  developed  from  known  experimental  data  of  one  environment 
to  another  environment.  This  is  an  accepted  modeling  technique  utilized  with 
dimensionless  heat  transfer  correlations  that  are  derived  from  experimental 
studies  and  is  a  simpler  alternative  to  modeling  the  entire  thermal  regulatory 
system. 

This  is  not  a  definitive  treatise  of  the  topic,  rather  it  is  an 
introduction  to  the  concept  of  developing  empirical  supplementary  heating 
!  correlations  for  divers  from  experimental  studies  of  diver  heat  loss  and 

i 

:  protective  garment  evaluations.  We  expect  that  the  correlations  suggested 
'  herein  will  be  improved  in  accuracy  as  additional  experimental  data  becomes 
available  and  that,  ultimately,  the  technique  will  be  extended  for  working 

i 

divers  when  suitably  controlled  experiments  have  been  conducted.  Until  these 
additional  experimental  studies  are  conducted,  we  recommend  that  the  equations 
presented  in  this  report  be  used  by  an  engineer  who  is  designing  a  regionally 
heated  thermal  protection  garment  for  a  resting  diver. 


Charles  E.  Johnson 
James  0.  Collins 


(vii) 


ACKNOWLEDGMENTS 


The  research  reported  herein  has  been  supported  under  Office  of  Naval 
Research  Contract  NQQQlWG-C-OSTjLwith  funds  provided  by  the  Naval  Medical 
Research  and  Development  Command  (NMRDC).  Personal  thanks  are  extended  to  the 
consecutive  Submarine  and  Diving  Medicine  Program  Managers  of  NMRDC 

:  Captain  James  Vorosmarti,  Jr.,  (MC)USN,  and  Captain  Kristopher  Greene,  (MC)USN, 

| 

!  for  their  interest  and  support  of  this  research  topic. 

It  would  not  have  been  possible  to  conduct  the  research  reported  herein 
without  the  assistance  and  cooperation  of  a  number  of  people  and 

i 

I  organizations.  Among  the  many  people  who  have  provided  assistance,  we  would 
'  like  to  extend  special  thanks  to  the  following: 

Dr.  Marshall  L.  Nuckols,  Naval  Coastal  Systems  Center,  for  the  time  he 
|  spent  arranging  our  acquisition  of  the  prototype  Diver  Thermal  Protection 
!  garment  evaluation  studies  and  for  personal  consultations  about  the  execution 
j  of  the  experimental  studies.  We  also  thank  Dr,  Nuckols  for  arranging  an  August 
!  1980  visit  of  James  D,  Collins  to  the  Navy  Expc  imental  Diving  Unit  of  Panama 
City,  Florida  for  the  purpose  of  witnessing  a  Diver  Thermal  ProgecMon  garment 
evaluation  study. 


r 


(viii) 


;  Or.  Claude  Piantadosi  for  providing  the  respiratory  heat  loss  studies 
j  that  he  performed  while  attached  to  the  Navy  Experimental  Diving  Unit  of 
Panama  City,  Florida.  Or.  Piantadosi  was  also  extremely  helpful  in  explaining 
the  physiological  phenomenon  which  are  associated  with  this  recorded  skin 
temperature  and  heat  flux  data. 

Mr.  Shiangtai  Tuan  of  the  Duke  University  Computation  Center,  Mr.  Steve 
Vickery  of  the  Navy  Experimental  Diving  Unit,  and  Dr.  Larry  Burton  and 
Ms.  Karen  Shephard  of  the  Department  of  Anesthesiology,  Duke  University 
Medical  Center,  for  their  help  in  the  transfer  of  data  files  from  the  Navy 

I 

!  Experimental  Diving  Unit's  Hewlett-Packard  21MX  computer  to  the  Department  of 

J 

i 

Anesthesiology's  Digital  Equipment  Corporation  11/60  computer.  We  extend 
:  additional  thanks  to  Ms.  Karen  Shephard  for  her  help  in  our  Interacting  with 

!  the  Department  of  Anesthesiology's  computer  system. 

i 

! 

) 

Ms.  Judy  Milne  and  Ms.  Melodie  Feather  for  their  time  spent  in  entering 
;  and  verifying  the  prototype  Diver  Thermal  Protection  garment  data. 

i 

■  Ms.  Melodie  Feather  also  typed  and  proofed  this  document  and  contributed  many 
'  valuable  editorial  suggestions, 

j  Cdr.  John  Zumrick  (MC)USN  of  the  Navy  Experimental  Diving  Unit  and 

i  Mr.  Max  Lippitt  of  the  Naval  Coastal  Systems  Center  for  the  information  they 
provided  during  the  visit  to  the  avy  Experimental  Diving  Unit  in  August  1980. 
,  Additional  thanks  are  extended  to  Dr.  Zumrick  for  acting  as  the  official 
j  sponsor  of  this  visit. 


We  are  also  very  appreciative  of  the  continued  support  of  the  following 
organi zations: 

The  staff  of  the  F.G.  Hall  Environmental  Laboratory,  Duke  University 
Medical  Center 

The  staff  of  the  Departure  of  Anesthesiology,  Duke  University  Medical 
Center,  e$oeci>lly  tre  Anesthesiology  computer  group 

The  staff  ol  the  Department  of  Mechanical  Engineering  and  Matierals 
Science,  especially  Mri.  °aul  i»c*  .  Roberts  for  editorial  assistance  in 
the  preparation  of  reports  and  manuscripts 


(*) 


TABLE  OF  CONTENTS 


Page 


ABSTRACT 

ii 

FOREWORD 

iv 

ACKNOWLEDGMENTS 

vi  i  i 

LIST  OF  TABLES 

xiii 

LIST  OF  FIGURES 

xvi  i  i 

CHAPTERS 

I.  INTRODUCTION  2 

Problem  Background,  2 
Objective,  4 
Approach,  4 

II.  AVAILABLE  EXPERIMENTAL  DATA  AND  METHODS  Of  ANALYSIS  APPLIED  14 

Generation  of  an  Example  Average  Subject  file,  21 
Description  of  Data  Groups,  75 
Siswary,  S3 

fit.  DEVELOPMENT  Of  THE  DIVER  MODEL  35 

The  Physical  Model,  R5 
The  Model  and  Its  Solution,  56 

IV.  DEVELOPMENT  Of  PREDICTIVE  SKIN  TEMPERATURE  CORRELATION 

Development  of  Ecuat^ons.  9; 

Authentication  of  the  Model,  U? 


«*n 


94 


(Table  of  Contents,  continued) 


Page 

V.  DEVELOPMENT  OF  OVERALL  HEAT  TRANSFER  COEFFICIENT  132 

Resistance  of  the  Outer  Boundary  Layer,  133 
Undergarment  Thermal  Resistance,  142 

VI.  THE  INTERACTIVE  MODEL  PROGRAM  152 

Main  Program,  156 
REINIT,  157 
OHTC,  157 
F(T),  158 

CriTC,  159  I 

TR,  163  i 

PRED,  163 
TC,  163 
PLOTD,  164 

!  PLOTN,  164  j 

!  i 

|  VII.  MODEL  VERIFICATION,  CONCLUSIONS,  AND  RECOMMENDATION  166  j 

i 

Empirical  Correlation  of  the  Model,  166  J 

Discussion  and  Recommendation,  174  J 


APPENDICES 


A. 

ORGANIZATION  OF  THE  DIRECT  ACCESS  EXPERIMENTAL  DATA  FILES 

182 

B. 

COMPOSITE  DATA  FILE  LISTINGS  FOR  FILES  RHL1AR1,  DTP1ART, 

DTP2AR1,  AND  DTP3AR1 

184 

C. 

DETERMINATION  Of  THE  VISCOSITY  (u),  THERMAL  CONDUCTIVITY  (k), 

AND  OTHER  PROPERTIES  FOR  A  95%  HELIUM  -  5%  OXYGEN  GAS  MIXTURE 

205 

D. 

INTERACTIVE  PROGRAM  DOCUMENTATION 

213 

E. 

INTERACTIVE  PROGRAM  FLOW  CHART 

269 

F. 

INTERACTIVE  PROGRAM  LISTING 

303 

G, 

SAMPLE  CALCULATION 

344 

REFLRE 

NCES 

366 

ADDENDUM:  MIXED  BREATHING  GAS  INTERACTIVE  MODEL  PROGRAM 

371 

m&msmamw  map 


1 


;** 


vn'i 


«  ‘  'ta 

A** 


Number 


LIST  OF  TABLES 


Title 


Indices,  Nomenclature,  Abbreviations,  and  Surface  Area 
Fractions  for  Hody's  Twelve  Segment  Map 

Respiratory  Heat  Loss  Studies  (RHL) 

Diver  Thermal  Protection  Garment  Evaluations 

Diver  Thermal  Protection  Garment  Evaluations 

Model  Ma  in  nsions 

Comparison  of  Model  Man  and  Hody's  Man  Regional  Surface 
Area  Percents 

Kerslake's  Curve  Fitting  Coefficients  A^ 

Comparison  of  Composite  Experimental  Temperatures 
Extracted  From  Figures  52-63  With  Those  Predicted  With 
Kerslake's  Equation  (Eq.  7) 

Comparison  of  Composite  Experimental  Temperatures 
Estimated  by  Means  of  Respective  Linear  Regression 
Equations  Which  Were  Derived  From  Figures  52-63  With 
Those  Predicted  With  Kerslake's  Equation  (Eq.  7) 

Distribution  Fractions  for  Predictive  Skin  Temperature 
Correlations 

Percentage  of  the  Total  Number  of  Compared  Points  (56) 
Over  the  Span  of  35  to  90  Minutes  Which  Have  Absolute 
Temperature  Differences  <  0.50  °C  and  <  0.25°C 


Page 


118  ; 


(xiii) 


(List  of  Tables,  continued) 


Number  Title  Page 


12  Time  Averaged  (Mean)  Normalized  Temperatures  for  Each  124 

of  Hody's  Twelve  Segments  Computed  From  the  Mean 

Temperature  Profiles  of  Experimental  File  RHL1AR1  and 
the  Equivalent  Predicted  Normalized  Temperature  Which 
was  Estimated  by  Equations  10  and  11 

13  Time  Averaged  (Mean)  Normalized  Temperatures  for  Each  125  i 

of  Hody's  Twelve  Segments  Computed  From  the  Mean  j 

Temperature  Profiles  of  Experimental  File  DTP1ART  and  ! 

the  Equivalent  Predicted  Normalized  Temperature  Which 

was  Estimated  by  Equations  10  and  11  j 

14  Time  Averaged  (Mean)  Normalized  Temperatures  for  Each  126  J 

of  Hody's  Twelve  Segments  Computed  From  the  Mean 

Temperature  Profiles  of  Experimental  File  DTP2AR1  and 
the  Equivalent  Predicted  Normalized  Temperature  Which 
was  Estimated  by  Equations  10  and  11 

15  Tima  Averaged  (Mean)  Normalized  Temperatures  for  Each  127 

of  Hody's  Twelve  Segments  Computed  From  the  Mean 

Temperature  Profiles  of  Experimental  File  DTP3AR1  and 
the  Equivalent  Predicted  Normalized  Temperature  Which 
was  Estimated  by  Equations  10  and  11 

16  Model  Man  Region  Names,  Indices,  and  Corresponding  Hody  131 

Segments 

17  Empirical  Convective  Heat  Transfer  Correlations  141 

18  Regional  Allowable  Heat  Fluxes  154 

19  a  as  a  Function  of  Temperature  161 

20  Thermal  Conductivity  of  Water  as  a  Function  of  162 

Temperature 

21  Time  Averaged  (Mean)  Biot  Numbers  for  Each  Model  Region  170  ! 

(Fig.  1)  Computed  From  the  Composite  Profiles  of  ; 

Experimental  File  0TP1ART  and  the  Equivalent  Predicted  j 

Biot  Number 

i 

22  Time  Averaged  (Mean)  Biot  Numbers  for  Each  Model  Region  171  I 

(Fig.  1)  Computed  From  the  Composite  Profiles  of  . 

Experimental  File  DTP2AR1  and  the  Equivalent  Predicted  i 

Biot  Number 


(xiv) 


(List  of  Tables,  continued) 


Number 

Title 

Pa^e 

23 

Time  Averaged  (Mean)  Blot  Numbers  for  Each  Model  Region 
(Fig.  1)  Computed  From  the  Composite  Profiles  of 
Experimental  File  DTP3AR1  and  the  Equivalent  Predicted 

Biot  Number 

172 

24 

Mean  Experimental  and  Predicted  Biot  Numbers  Averaged 

From  the  Respective  Biot  Numbers  of  Tables  21,  22,  and 

23  From  Files  DTP1ART,  DTP2AR1,  and  DTP3AR1, 

Respectively 

173 

25 

Geometrical  Data  for  the  Thigh 

177 

B.l,  A 

Temperature  C  RHL1AR1 

185 

B.l.B 

Rate  of  Heat  Loss  (Watts)  RHL1AR1 

188 

B.2.A 

Temperature  C  DTP1ART 

191 

B.2.B 

R<,te  "  Heat  Loss  (Watts)  DTP1ART 

193 

3.3.A 

Temperature  C  nTP2ARl 

195 

B.  3.8 

Ra*e  of  Heat  Loss  (We^ts)  DTP2AR1 

197 

B.4.A 

Temperat  .re  C  DTP3AR1 

199 

B.4.B 

Rate  of  Heat  Loss  (Watts)  DTP3AR1 

202 

ca 

Critical  anti  Reduced  Properties,  and  Compressibility 

Factor  Values  for  a  95%  Helium  -  5%  Oxygen  Mixture  at 
295.30  °K  and  20.8?  ATA 

206 

C.2 

Properties  Substituted  into  Equations  C.5  and  C.6  for  a 

95%  Helium  -  5%  Oxygen  Mixture  at  295.3  °K  and 

307.0  PSIA 

209 

C.3 

Properties  for  Equations  C.7  and  C.8  for  a  95%  Helium  - 
5%  Oxygen  Mixture  at  29o.3  °K  > -i  307,0  PSIA 

211 

C.4 

Properties  Necessafy  to  Detennir<>  drashof,  Prandtl,  and 
Nusselt  Numbers  for  a  95%  Helium  -  5%  Oxygen  Mixture  at 
22.15  °C  and  200  MSW 

212 

0.1 

Production  Mode  Question  Flow 

233 

0.2 

Sample  Production  Mode  Output 

233 

(xv) 


I  (List  of  Tables,  continued) 


Number 

Title 

0.3 

Research  Mode  Question  Flow 

234 

D.4 

Sample  Research  Mode  Output 

235 

0.5 

Experimental  (T-EXPT)  and  Predicted  (T-PREO) 

Temperature  Values  (°C)  as  a  Function  of  Time  as 

Plotted  in  Figure  D.l 

238 

0.6 

Experimental  (Q-EXPT)  and  llncorrected  Predicted 
(Q-PRED)  Heat  Flux  Values  (W/M^)  as  a  Function  of  Time 
as  Plotted  in  Figure  D.2 

241 

0.7 

Experimental  (ll-EXPT)  and  Uncorrected  Predicted 
(U-PRE0)  Overall  Heat  Transfer  Coefficient  Values 
(W/m2°C)  as  a  Function  of  Time  as  Plotted  in 

Figure  0.3 

244 

D.8 

Experimental  (TN-EXPT)  and  Predicted  (TN-PRED) 

Normalized  Temperature  Values  (Dimensionless)  as  a 

Function  of  Time  as  Plotted  in  Figure  0.4 

247 

D.9 

Experimental  (8I-EXPT)  and  Uncorrected  Predicted 
(BI-PRED)  Biot  Humber  Values  (Dimensionless)  as  a 

Function  of  Time  as  Plotted  in  Figure  0.5 

250 

0.10 

Sample  Research  Mode  Output  With  Corrected  Overall  Heat 
Transfer  Coefficients 

253 

0.11 

Experimental  (T-EXPT)  and  Predicted  (T-PRED) 

Temperature  Values  (°C)  as  a  Function  of  Time  as 

Plotted  in  Figure  0.6 

256 

0.12 

Experimental  (Q-EXPT)  and  Corrected  Predicted  (Q-PRED) 

Heat  Flux  Values  (W/M^)  as  a  Function  of  Time  as 

Plotted  in  Figure  0.7 

259 

0.13 

Experimental  (U-EXPT)  and  Corrected  Predicted JU-PRtO) 
Overall  Heat  Transfer  Coefficient  Values  (W/M^C)  as  a 
Function  of  Time  as  Plotted  in  Figure  D.8 

26  2 

0.14 

Experimental  (TN-EXPT)  and  Predicted  (TN-PRED) 

Normalized  Temperature  Values  (Dimensionless)  as  a 

Function  of  Time  as  Plotted  in  Figure  0.9 

265 

(*vi) 


j  (List  of  Tables,  continued) 


1  Number  Title 


Page 


0.15  Experimental  (BI-EXPT)  and  Corrected  Predicted  268 

(BI-PRED)  Biot  Number  Values  (Dimensionless)  as  a 
Function  of  Time  as  Plotted  in  Figure  D.10 

G.l  A  Production  Mode  Output  Enhanced  With  Debug  345 

Information 

6.2  Question  Flow  for  Generation  of  Table  G.l  346 


•  ^ 


(svi i ) 


LIST  OF  FIGURES 


Title 


Average  Diver  Modeled  as  Connected  Cylinders  and  3 

Spheres 

Hody's  Twelve  Segments  11 

Segmental  Subject  Temperature  Profiles  for  the  Head  24 

(HEAD) 

Segmental  Subject  Temperature  Profiles  for  the  Chest  25 

(CHST) 

Segmental  Subject  Temperature  Profiles  for  the  Abdomen  26 

(ABD) 

Segmental  Subject  Temperature  Profiles  for  the  Upper  27 

Back  (UBAC) 

Segmental  Subject  Temperature  Profiles  for  the  Lower  28 

Back  (L6AC) 

Segmental  Subject  Temperature  Profiles  for  the  Arm  29 

(ARM) 

Segmental  Subject  Temperature  Profiles  for  the  Hand  30 

(URST) 

Segmental  Subject  Temperature  Profiles  for  the  Front  31 

Thigh  {F-TH} 

Segmental  Subject  Temperature  Profiles  for  the  Front  32 

Calf  (FCIFJ 

Segmental  Subject  Temperature  Profiles  for  the  Rear  33 

Thigh  (S-TH) 


(xvi i  1 ) 


(List  of  Figures,  continued) 


Number 

Title 

fisi 

13 

Segmental  Subject  Temperature  Profiles  for  the  Rear 

Calf  (RCLF ) 

34 

14 

Segmental  Subject  Temperature  Profiles  for  the  Foot 
(FOOT) 

35 

15 

Segmental  Subject  Sate  of  Heat  Loss  Profiles  for  the 

Head  (HEAD) 

37 

16 

Segmental  Subject  Rate  of  Heat  Loss  Profiles  for  the 

Chest  (CHST) 

38 

17 

Segmental  Subject  Rate  of  Heat  Loss  Profiles  for  the 
Abdomen  (ABO) 

39 

18 

Segmental  Subject  Rate  of  Heat  Loss  Profiles  for  the 

Upper  Back  (UBAC) 

40 

19 

Segmental  Subject  Rate  of  Heat  Loss  Profile-':  for  the 

Lower  Back  (L8AC) 

41 

20 

Segmental  Subject  Rate  of  Heat  Loss  Profiles  for  the 

Arm  (ARK) 

42 

21 

Segmental  Subject  Rate  cf  Heat  Loss  Profiles  for  the 

Hand  (WST) 

43 

22 

Segmental  Subject  Rate  of  Heat  Loss  Profiles  for  the 

Front  Thigh  (F-TH) 

44 

23 

Segmental  Subject  Rate  of  Heat  Loss  Profiles  for  the 

Front  Calf  (FCLf) 

45 

24 

Segmental  Subject  Rate  of  Heat  Loss  Profiles  for  the 

Rear  Thigh  (R-TH) 

46 

25 

segmental  Subject  Rate  of  Heat  Loss  Profiles  for  the 

Rear  Calf  (RCLF) 

47 

26 

Segmental  Subject  Rate  of  Heat  Loss  Profiles  for  the 

Foot  (FOOT) 

48 

2? 

Segmental  Subject -Averaged  Temperature  Profile  for  the 
Head  (HEAD) 

50 

(*i*) 


i 


|  (List  of  Figures,  contl out,  , 

i  Number  Title  Page 


28 

Segmental  Subject-Averaged 
Chest  (CHST) 

25 

Segmental  Subject -Averaged 
Abdomen  (ABO) 

30 

Segmental  Subject -Averaged 
Upper  Back  (UBAC) 

31 

Segmental  Subject-Averaged 
Lower  Back  (LBAC) 

32 

Segmental  Subject-Averaged 
Arm  (ARM) 

33 

Segmental  Subject -Averaged 
Hand  (WRST) 

34 

Segmental  Subject -Averaged 
Front  Thigh  (f-TH) 

35 

Segmental  Subject -Averaged 
Front  Calf  (FCtF) 

36 

Segmental  Subject -Averaged 
Rear  Thigh  (R-TH) 

37 

Segmental  Subject-Averaged 
Rear  Calf  (RCL?) 

38 

Segmental  Subject -Averaged 
Foot  (FOOT) 

39 

Segmental  Subject-Averaged 
for  the  Head  (HEAD) 

40 

Segmental  Subject -Averaged 
for  the  Chest  (CHST) 

41 

Segmental  Subject -Averaged 
for  the  ftdomn  (ABO) 

42 

Segmental  Subject -Averaged 
for  the  Upper  Back  (OSAC) 

Temperature  Profile  for  the  51 
Temperature  Profile  for  the  52 
Temperature  Profile  for  the  53 
Temperature  Profile  for  the  54 
Temperature  Profile  for  the  55 
Temperature  Profile  for  the  56 
Temperature  Profile  for  the  57 
Temperature  rrofile  for  the  58 
Temperature  Profile  for  the  59 
Temperature  Profile  for  the  60 
Temperature  Profile  for  the  61 
Rate  of  Heat  Loss  Profile  63 
Rate  of  Heat  Loss  Profile  64 
Rate  of  Heat  Loss  Profile  65 
Rate  of  Heat  Loss  Profile  66 


{«) 


(List  of  Figures,  continued) 


Number  Title 


43  Segmental  Subject-Averaged  Rate  of  Heat  Loss  Profile 
for  the  Lower  Back  (LBAC) 

44  Segmental  Subject-Averaged  Rate  of  Heat  Loss  Profile 
for  the  Arm  (ARM) 

45  Segmental  Subject-Averaged  Rate  of  Heat  Loss  Profile 
for  the  Hand  (WRST) 

46  Segmental  Subject -Averaged  Rate  of  Heat  Loss  Profile 
for  the  Front  Thigh  (F-TH) 

4?  Segmental  Subject-Averaged  Rate  of  (teat  Loss  Profile 

for  the  Front  Calf  (FCLF) 

48  Segmental  Subject -Averaged  Rate  of  Heat  Loss  Profile 
for  the  Rear  Thigh  (R-TH) 

49  Segmental  Subject-Averaged  Rate  of  Heat  Loss  Profile 
for  the  Rear  Calf  (RCLF) 

50  Segmental  Subject -Averaged  Rate  of  Heat  Loss  Profile 
for  the  Foot  (FOOT) 

51  Vertical  Cylindrical  Region  Divided  into  Five  Nodes 

(>$) 

52  Segmental  Temperature  Profile  for  the  Head  (HEAD)  and 
the  Nody  Mean  Skin  Temperature  (HMST) 

53  Segmental  Temperature  Profile  for  the  Che^t  (CHST)  and 
the  Hody  Mean  Skin  Temperature  (HHST) 

54  Segmental  Temperature  Profile  for  the  Abdomen  (A8D)  and 
the  Hody  Ksan  Skin  Temperature  (HMST) 

55  Segmental  Temperature  Profile  for  the  Upper  Back  (U8AC) 
and  the  Hody  Kean  Skin  Temperature  (KH$T) 

56  Segmental  Temperature  Profile  for  the  Lower  Back  (L6AC) 
and  the  Hody  Mean  Skin  Temperature  (HHST) 

5?  Segmental  Temperature  Profile  fo r  the  Arm  (AJ?H)  and  the 

Hody  Kean  Skin  Temperature  (HHST) 


(List  of  Figures,  continued) 


Number  Title  Page 


58  Segmental  Temperature  Profile  for  the  Hand  (WRST)  and  106  ! 

the  Hody  Mean  Skin  Temperature  (HMST)  j 

59  Segmei.Ul  Temperature  Profile  for  the  Front  Thigh  107  ! 

(F-TH)  and  the  Hody  Mean  Skin  Temperature  (HMST) 

60  Segmental  Temperature  Profile  for  the  Front  Calf  (FCLF)  108 

and  the  Hody  Mean  Skin  Temperature  (HMST) 

i 

61  Segmental  Temperature  Profile  for  the  Rear  Thigh  (R-TH)  109 

and  the  Hody  Mean  Skin  Temperatur0  (HMST) 

62  Segmental  Temperature  Profile  for  the  Rear  Calf  (RCLF)  110  ' 

and  the  Hody  Mean  Skin  Temperature  (HMST)  ! 

63  Segmental  Temperature  Profile  for  the  Foot  (FOOT)  and  111  j 

the  Hody  Mean  Skin  Temperature  (HMST) 

64  Linearized  Abdominal  (ABD)  and  Hody  Mean  Skin  (HMST)  113  j 

Temperatures  i 

65  Log  of  Nusselt  Number  Versus  Log  of  Grashof  Number  and  136  s 

Prandtl  Number  Product  for  Vertical  Cylinder 

66  Log  of  Nusselt  Number  Versus  Log  of  Grashof  Number  and  137  , 

Prandtl  Number  Product  for  Horizontal  Cylinder 

£/  Log  d  Nusselt  Number  Versus  Log  of  Grashof  Number  and  138 

Prandtl  frsftbe*  Product  for  Sphere 

68  Linear  Relationship  Describing  the  Undergarment's  145 

Specific  Thermal  Resistance  Change  as  a  Function  of 

Change  in  Thickness 

69  Percent  Change  in  Thickness  as  a  Function  of  148 

Compressive  Pressure 

70  Percent  Change  in  Thermal  Resistance  as  a  Function  of  149 

Compressive  Pressure 

71  Lateral  Mid-lines  for  Horizontal  Cylinder  and  Sphere  153 

72  Single  Nodal  Element  for  a  Vertical  Body  Region  and  155 

Kulti-nodal  Element  for  a  Vertical  Body  Region 

73  Representative  Region  Presented  as  a  Truncated  Cone  178 

(xxii) 


(List  of  Figures,  continued) 


Number 

Title 

Paje 

D.l 

Plots  of  Experimental  (T-EXPT)  and  Predicted  (T-PRED) 
Temperatures  (°C)  for  the  Abdomen 

237 

0.2 

Plots  of  Experimental  (Q-EXPT)  and  Uncorrected 

Predicted  (Q-PRED)  Heat  Flux  (W/H^)  for  the  Abdomen 

240 

D.3 

Plots  of  Experimental  (U-EXPT)  and  Uncorrected 

Predicted  (U-PRED)  Overall  Heat  Transfer  Coefficient 
(OHTC  in  W/m2°C)  for  the  Abdomen 

243 

0.4 

Plots  of  Experimental  Normalized  (TN-EXPT)  and 

Predicted  Normalized  (TN-PRED)  Temperatures  for  the 

Abdomen 

246 

j 

0.5 

Plots  of  Experimental  Biot  (BI-EXPT)  and  Uncorrected 
Predicted  Biot  (BI-PRED)  Numbers  for  the  Abdomen 

249  ' 

0.6 

Plots  of  Experimental  (T-EXPT)  and  Predicted  (T-PRED) 
Temperatures  (°C)  for  the  Abdomen 

255 

0.7 

Plots  of  Experimental  (Q-EXPT)  and  Corrected  Predicted 
(Q-PRED)  Heat  Flux  (W/H*)  for  the  Abdomen 

258 

0.8 

Plots  of  Experimental  (U-EXPT)  and  Corrected  Predicted 
(U-PRED)  Overall  Heat  Transfer  Coefficient  (OHTC  in 

W/fi2°C)  for  the  Abdomen 

261 

0.9 

Plots  of  Experimental  Normalized  (TN-EXPT1  and 

Predicted  Normalized  (TN-PRED)  Temperatures  for  the 

Abdomen 

T84 

D.10 

Plots  of  Experimental  Biot  (31-EXPT)  and  Corrected 
Predicted  Biot  (BI-PRED)  Numbers  for  the  Abdomen 

267 

(xxiii ) 


THE  DEVELOPMENT  OF  PREDICTIVE  ENGINEERING  FORMULATIONS  FOR  DIVER  HEATING 


CHAPTER  I 


INTRODUCTION 


Problem  Background 

There  is  virtually  no  situation  in  manned  hyperbaric  or  hydrostatic 
operations  where  the  effects  of  metabolic  heat  loss  are  not  present.  The 
extent  and  effect  of  metabolic  heat  loss  from  divers  in  various  operational 
environments  from  temperate  to  arctic  water  conditions  are  highly 
documented  [l-5j.  The  problems  associated  with  cold  and  its  physiological  and 
psychological  impairments  are  common  to  all  forms  of  diving:  salvage, 
free-swimming  special  warfare,  tethered  deep  ocean,  and  surface  supplied 
shallow  water.  Fatal  mishaps  in  diving  operations  such  as  the  SeaLab  III  and 
the  Link  tragedies  only  remind  us  of  the  insidious  influence  of  cold 
temperatures  in  the  alien  sea. 

Scientists  have  been  studying  these  problems  for  *any  years.  As  early  as 
1966,  attempts  were  made  to  develop  analytical  models  of  the  buwan  thermal 
regulatory  system  [6,7?.  In  1971  Rayne  and  Austin  [8]  proposed  a  set  of 
design  guidelines  for  the  thermal  protection  of  divers.  In  J9?6  the  advisory 
and  review  group  which  evaluated  thermal  problems  for  the  national  Rian  for 


{2} 


r 


3  1 

the  Safety  and  Health  of  Divers  in  the  Quest  for  Subsea  Energy  [9]  pi aced  the  ; 

I 

the  highest  priority  on  research  involved  in  the  prediction  of  and 

j 

compensation  for  diver  heat  loss  under  various  marine  conditions.  The  primary  ' 
recommendation  of  this  committee  was  to  establish  f  set  of  interim  guidelines 
■  for  human  tolerance  in  hyperbaric  conditions  based  upon  a  review  of  the  j 

[  j 

I  existing  literature.  This  was  accomplished  by  Webb  et  al .  [10]  in  1976,  with  j 

j  a  revision  appearing  in  1980  [11],  However,  both  of  these  proposed  guidelines  ' 

! 

for  diver  safety  are  expressed  in  purely  physiological  terms  of  maximum  . 

allowable  heat  loss  quanta  and  minimum  allowable  skin  temperatures  and,  thus, 
do  not  readily  lend  themselves  to  use  by  a  designer  of  diver  thermal  > 

protection  equipment. 

I  The  second  recommendation  of  the  thermal  committee  for  the  National  ! 

! 

j  Plan  [9]  suggested  a  complete  experimental  and  analytical  study  which  would 
consist  of: 

A.  A  serious  study  of  skin  and  respiratory  heat  exchange  !n  divers  to 
provide  data  that  would  be  useful  both  operationally  and  in  the  design  of 
thermal  protection  equipment.  This  data  would  perhaps  enable  the 
establishment  of  thermal  limits  for  depth  or  duration  of  exposure  either  in 
water  or  in  dry  habitats. 

8,  An  analytical  approach  to  devise  predictive  equations  based  on 
accurate  heat  transfer  data  so  that  body  heat  loss  and  temperature  can  be 
predicted  for  any  environment.  A  well -organized  and  comprehensive 
experimental  program  to  develop  and  test  such  predictive  equations  is  needed. 


4 


1  Objective 

The  objective  of  this  research  has  been  to  develop  formulations  for  j 

predicting  the  regional  supplementary  heating  required  to  maintain  a  resting, 

i 

nearly  euthermic  diver.  Assuming  that  the  respiratory  heat  loss  is  conserved, 

I 

I 

we  hypothesized  that  the  regional  supplementary  heating  requirement  could  be  j 
estimated  as  the  difference  between  the  mandatory  regional  heat  flux1  due  to 
metabolism  and  that  predicted  for  an  unheated  immersed  diver.  Furthermore,  we 

I 

assumed  that  the  predictive  correlations  for  the  regional  heat  flux  could  be  j 

I 

developed  by  utilizing  an  overall  heat  transfer  coefficient  to  describe  the 
|  thermal  behavior  of  the  garment  ensemble.  We  hypothesized  that  this 

!  coefficient  could  be  developed  by  considering:  a)  the  response  of  the 

i  1 

|  garment's  insulating  materials  to  compression  [12]  and  permeation  [13]  with 

i 

I  conventional  breathing  gases  and  b)  the  geometrical  parameters  of  individual 

i 

j  body  regions  associated  with  three  specific  postures:  standing,  sitting,  and 

i 

!  prone.  Pursuing  the  expressed  hypothesis,  we  postulated  that  the  resultant 

i  ; 

!  predictive  correlations  for  diver  regional  heat  flux  could  be  corrected  by 

i 

|  using  an  empirically  der tved  biasing  factor.  This  factor  would  be  obtained  by 

» 

comparing  the  respective  predicted  overall  heat  transfer  coefficient  with  mean 
I  experimental  values  obtained  from  the  U.S.  Navy's  prototype  Diver  Thermal 
•  Protection  (DTP)  garment  evaluation  studies. 


Approach 

Traditional ’y,  mathematical  models  attempting  to  predict  physiological 
heat  loss  consider  segments  of  the  human  body  as  cylindrical  elements  (Head, 


‘See  Foreword  for  connotations  of  region  and  Heat  flux. 


torso,  upper  arm,  lower  arm,  hand,  thigh,  etc).  In  the  case  of  the  Immersed 
diver,  heat  is  considered  to  be  lost  from  the  cylindrical  surfaces, 
predominantly  by  convection.  Metabolic  heat  is  generated  within  the  mass  of 
each  segment  and  is  exchanged  among  the  elemental  cylinders  by  means  of  a 
modeled  circulatory  system.  A  portion  of  this  heat  is  transported  to  the 
cylindrical  surface  by  conduction  through  the  tissue  and  by  convective 
transport  through  a  local  circulatory  system  [6,7],  The  various  modes  of 
transport  approximating  the  physiological  heat  exchange  process  are 
interrelated  by  »at  and  mass  transfer  balances.  The  coefficients  of  local 
metabolic  heat  generation,  tissue  conductance,  and  circulatory  transport  are 
either  arbitrarily  assumed  or  empirically  derived  in  such  equations.  Some 
progress  has  been  made  in  identifying  these  parameteric  coefficients  by 
comparing  the  model  to  experimental  data  [14],  Body  temperature  and  heat  flux 
may  then  be  predicted  for  conditions  analogous  to  the  experimental  situation. 

The  accuracy  of  the  discussed  type  of  physiological  heat  exchange  model 
is  highly  dependent  upon  the  selection  of  heat  distribution,  transport,  and 
generation  parameters.  However,  these  parameters  vary  quite  significantly  as 
a  function  of  the  state  of  thermal  regulation,  metabolic  activity,  and  body 
region.  For  instance,  tissue  thermal  conductance  may  vary  by  a  factor  of 
fifteen  between  the  states  of  vasodilation  and  vasoconstriction  [15].  Even 
the  application  of  hot  and  cold  stimuli  to  the  peripheral  limbs  has  been 
observed  to  cause  changes  in  the  vasomotor  state  of  the  uppep  respiratory 
system  [16-13].  Thus,  we  infer  that  the  variation  of  tissue  conductance  and 
the  transport  of  heat  by  the  local  inter-regional  circulation  are  quite 
difficult  to  predict  as  a  function  of  either  the  thermal  environment  or  some 
metabolic  indicator  (e.g.,  skin  temperature  or  local  heat  flux). 

A  simpler  alternative  to  modeling  the  entire  thermal  regulatory  system 
and  the  thermal  resistance  imposed  by  the  garment  is  to  model  only  the  heat 


flux  through  the  Insulatlve  garment.  The  thermal  resistance  imposed  by  this 
garment  is  a  function  of  the  ambient  environment  and  the  thermal  state  of  the 

subject  as  reflected  at  the  skin  surface  of  each  region.  This  concept  relies 

i 

heavily  upon  existing  experimental  data  of  regional  skin  temperature  and  heat 
flux  obtained  from  subjects  In  various  amoient  environments  and  states  of 
|  thermal  equilibrium.  Through  a  comparative  mathematical  analysis,  the 
regional  heat  flux  correlations  would  be  developed.  The  predicted  regional 
heat  flux  would  then  be  used  to  estimate  the  regional  supplementary 
heating  required  to  maintain  a  subject  in  a  state  of  euthermia.  These  heat 
fluxes  could  also  be  employed  to  predict  the  most  prudent  distribution  of  a 
i  limited  amount  of  available  supplementary  heating  while  considering  the 

I  limited  heat  conserving  mechanisms  of  the  human. 

! 

i  As  noted  above,  this  simpler  approach  relies  heavily  upon  existing 

l 

j  experimental  data.  An  extensive  search  was  performed  to  Identify  sources  of 

I 

;  experimental  regional  skin  temperature  and  heat  flux  data.  Unfortunately,  the 

i 

i  data  available  in  the  open  literature  was  determined  to  be  of  limited  value 
|  since  most  of  the  data  is  presented  is  terms  of  mean  skin  temperatures  rather 

i 

i 

than  the  desired  regional  values  or  Is  presented  in  terras  of  experimental 
1  .  regional  temperatures  and  computed,  not  experimental,  heat  fluxes.  This  lack 
[  of  suitable  regional  experimental  temperatures  and  heat  fluxes  prompted  our 
requests  to  the  U.S.  Navy  for  diver  experimental  data  recorded  with  calibrated 
instrumentation  certified  by  the  principal  investigator.  This  request  for 

< 

only  investigator  certified  instrumentation  resulted  from  questions  concerning 
the  accuracy  of  certain  commonly  used  heat  flux  transducers  when  only  the 
manufacturer's  calibration  was  employed  [19,20].  Our  request  to  the  U.S.  Navy 
identified  two  types  of  experimental  studies.  The  first  set  of  studies  was 
perforraed  by  the  Navy  Experimental  Diving  Unit  to  investigate  respiratory  heat 
loss  in  a  hyperbaric  environment  [21],  Sitting,  nearly  nude  subjects  were 
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exposed  to  hyperbaric  helium-oxygen  environments  at  a  depth  equivalent  to 
200  meters  of  seawater  (MSW).  These  experiments  consisted  of  resting  dives  at 
15  °C  and  20  °C  ambient  gas  temperatures  and  working  dives  in  a  20  °C 
environment. 

The  second  series  of  studies  was  performed  by  the  Naval  Coastal  Systems 
Center  (NCSC),  in  conjunction  with  NEDU,  to  evaluate  the  Navy's  prototype 
Diver  Thermal  Protection  (OTP)  garment  [22].  Data  of  segmental ^  skin 
temperature  and  heat  flux  were  recorded  from  prone  immersed  subjects  breathing 
compressed  air.  This  data  was  collected  at  two  different  depths,  3  and 
21  meters  but  at  essentially  the  same  ambient  temperature,  5  *0.  During  all 
of  the  experiments,  the  dry  DTP  garment  was  pressure  compensated  with 
compressed  air. 

The  experimental  data  has  been  utilized  to  develop  a  model,  in  accordance 
with  the  stated  hypothesis,  for  predicting  regional  heat  flux  and  required 
Supplementary  heating.  The  modeled  diver  Is  comprised  of  connected  cylinders 
for  all  regions  except  the  head  and  hands  which  are  treated  as  spheres 
(Fig.  1).  The  dimensions  of  this  nine  region  man  were  selected  to  approximate 
those  of  the  standard  U.S.  Air  Force  man  [23]  and  to  produce  surface  area 
fractions  equivalent  to  the  respective  coefficients  of  the  Hody  mean  skin 
temperature  correlation  [24,25]. 

Equation  1  estimates  the  regional  rate  of  heat  loss  based  on  the 
regional  skin  surface  area  (Si),  overall  heat  transfer  coefficient  (U^),  and 
temperature  difference  between  the  skin  surface  (T^)  and  ambient  (Ta) 
temperatures. 

$1  •  UiMTi  -  V  O) 


^See  Foreword  for  connotation  of  segment. 


The  regional  skin  temperature,  T^,  Is  predicted  with  the  Hedy  mean  skin 
temperature  [24,25]  as  the  independent  argument  In  a  predictive  correlation. 
This  correlation  was  developed  by  adapting  Kerslake's  eight  region  predictive 
formula  [2]  to  Hody's  twelve  segments  by  using  the  data  of  the  first 
experimental  series  discussed  [21].  The  skin  surface  area  1*.  determined  from 
the  dimensions  assigned  to  the  cylindrical  or  spherical  region  of  interest. 

The  ambient  temperature  must  be  specified.  A  mean  skin  temperature  is  used  as 
the  independent  argument  of  the  predictive  temperature  formulations  since  the 
limits  of  diver  thermal  equilibrium  are  usually  specified  in  terms  of  a 
minimum  mean  skin  temperature,  and  the  temperatures  of  the  hand  and  foot  [11]. 

The  concept  of  mean  skin  temperature  is  predicated  on:  a)  The  true  mean 
temperature  is  the  mean  of  an  infinite  number  of  body  temperature  points,  and 
the  weighted  mean  of  some  finite  number  of  these  points  provides  a  good 
estimate  of  the  true  value;  and  b)  Local  areas  of  the  body  may  be  assumed  to 
be  nearly  homogeneous  with  respect  to  temperature.  The  sum  of  mean  wighted 
regional  temperatures  is  then  assumed  to  be  characteristic  of  the  mean  surface 
temperature  of  the  human  body.  The  weighting  factors  used  with  the  various 
regions  of  the  body  are  normally  equivalent  to  the  fraction  of  the  total 
surface  area  comprised  by  the  respective  region  [26],  Teichner  [26]  compared 
various  methods  which  were  cited  in  the  literature  for  calculating  mean  skin 
temperature.  He  found  that  a  six  or  seven  point  weighted  system  (Palmes  and 
Park,  and  Burton,  respectively)  and  an  unweighted  ten  point  mean  estimated  the 
mean  sktn  temperature  of  man  as  well  as  more  sophisticated  schemes  did. 
Teichner  concluded  that  no  more  than  six  values  are  required  to  provide  a  good 
estimate  of  mean  skin  temperature  and  that  little  is  gained  when  more  than  six 
are  utilised.  However,  the  statistical  reliability  of  an  estimate  of  true 
mean  temperature  is  improved  when  additional  points  are  considered. 


10 

The  mean  skin  temperature,  used  as  the  independent  variable  in  our 
predictive  skin  temperature  correlation,  is  calculated  by  using  the  Hody 
correlation  [24,25].  This  twelve  segment  model  was  chosen  over  other  formulae 
[26-29]  since  the  segmental  sensor  sites  specified  in  Body’s  map  (Fig.  2  and 
Table  1)  matched  those  used  by  the  Navy's  thermal  instrumentation  system  [30] 
in  acquiring  the  experimental  data. 

A  correlation  for  the  overall  heat  transfer  coefficient,  u,  of 
Equation  1,  was  developed  by  modeling  a  composite  garment  ensemble  after  the 
Navy's  prototype  OTP  garment.  It  will  be  shown  that  this  formulation  taxes 
into  account  the  garment  entrapped  gas,  ambient  temperature  and  depth,  the 
effects  of  local  hydrostatic  compression  of  the  garment,  and  the  garment's 
external  film  coefficient  of  heat  transfer.  Experimental  values  of  U,  to 
which  the  predicted  values  will  be  compared  were  obtained  by  averaging 
segmental  skin  temperatures  and  heat  fluxes  collected  during  the  senes  of 
studies  performed  to  evaluate  the  Navy's  prototype  Diver  Thermal  Protection 
garment. 

The  diver  heat  loss  model  has  been  written  as  an  interactive  FORTRAN 
computer  program  which  has  two  modes  of  execution:  production  and  research. 
The  program  allows  the  user  in  either  mode  to  vary  the  garment’s  thermal 
properties  and  the  assessed  values  of  allowable  regional  heat  flux.  The  user 
is  also  permitted  to  vary  the  assumed  overall  heat  transfer  coefficient  (U) 
biasing  factors  in  either  mode.  These  biasing  factors  wilt  be  shown  to 
correct  the  non-dimensional i zed  values  of  predicted  U  such  that  in  a 
comparison  with  non-dimensional  1  zed  experimental  values  the  percent  difference 
will  be  £  15i.  These  options  are  provided  to  facilitate  comparative 
investigation  without  requiring  experimental  studies. 


Table  1 


Indices,  Nomenclature,  Abbreviations,  and  Surface  Area  Fractions 
for  Hody's  Twelve  Segment  Map.* 


Segmental 

Index 

Segmental 

Name 

Abbreviation 

Surface  Area 
Fraction 

.1 

Head 

HEAD 

0.070 

2 

Chest 

CHST 

0.085 

3 

Abdomen 

ABO 

0.085 

4 

Upper  Back 

UBAC 

0.090 

5 

Lower  Back 

LBAC 

0.090 

6 

Arm 

ARM 

0.140 

7 

Hand 

URST 

0.050 

8 

Front  Thigh 

F-TH 

0.095 

9 

Front  Calf 

FCLF 

0.065 

10 

Rear  Thigh 

R-TH 

0.095 

ii 

Rear  Calf 

RCIF 

0.065 

12 

Foot 

FOOT 

0.070 

*Data  adapted  from  Zuarick  [30] 
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The  production  mode  is  designed  to  produce  estimates  of  the  regional  heat 
flux  and  required  regional  supplementary  heating  based  on  user  supplied 
information.  For  this  mode,  the  user  must  specify  a  mean  skin  temperature  and 
the  ambient  temperature  and  depth.  The  value  specified  for  the  mean  skin 
temperature  indicates  the  desired  thermal  state  of  the  subject.  The  research 
mode  allows  the  comparison  of  experimental  and  predicted  values  of  regional 
temperatures,  heat  fluxes,  and  overall  heat  transfer  coefficients;  the 
experimental  mean  skin  and  ambient  temperatures  are  used  in  calculating  the 
predicted  data. 

The  computer  code  has  been  documented  and  flow  charted,  as  well  as 
internally  documented  with  comment  statements.  The  manner  in  which  the 
interactive  program  has  been  written  should  allow  a  scientist  versed  both  in 
programming  with  the  FORTRAN  language  and  in  diver  thermal  protection 
principles  to  modify  the  program  to  consider  new  garment  materials  and  suit 
configurations. 


CHAPTER  II 


AVAILABLE  EXPERIMENTAL  DATA  AND  METHODS  OF  ANALYSIS  APPLIED 


As  noted  earlier,  the  experimental  data  available  in  the  open  literature 
was  found  to  be  of  limited  value.  Most  of  the  data  is  presented  in  terms  of 
mean  skin  temperatures  rather  than  the  desired  regional  values  or  is  presented 
in  terms  of  experimental  regional  temperatures  and  computed,  not  experimental, 
heat  fluxes.  As  a  result,  we  requested  diver  experimental  data  which  was 
recorded  with  investigator  certified  instrumentation  from  the  U.S.  Navy.  Two 
types  of  studies  were  made  available  for  our  use.  The  first  type  of  study  was 
performed  by  the  Navy  Experimental  Diving  Unit  (NEDU)  to  investigate 
respiratory  heat  loss  in  hyperbaric  environments  [21}.  Sitting,  nearly  nude 
subjects  were  exposed  to  hyperbaric  helium-oxygen  environmenvs  at  a  depth  of 
200  meters  of  seawater  (NSW)  and  ambient  temperatures  of  15  UC  and  20  °C.  The 
second  type  of  Study  was  performed  by  the  Naval  Coastal.  Systems  Center  (NCSC), 
in  conjunction  with  NEDU.  to  evaluate  the  Navy's  prototype  Div®r  Thermal 
Protection  (OTP)  garment  [13,22,31].  Prone,  OTP  clad  subjects  were  exposed  in 
water  environments  at  depths  of  3  and  21  meters  and  at  an  ambient  temperature 
of  5  °C.  In  both  cases  studies  were  performed  with  both  resting  and  wording 


divers  [21,13,22,31]. 
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Table  2  describes  the  various  experiments  performed  during  the 
respiratory  heat  loss  study  while  Tables  3  and  4  describe  the  OTP  garment 
evaluation  studies.  Each  table  is  organized  with  the  resting  dives  listed 
first,  followed  by  the  working  dives.  If  only  a  part  of  the  resting  subjects 
were  used  in  the  verification  phases  of  our  work,  these  subjects  are  listed 
first. 

The  usefulness  of  each  of  the  experimental  data  files-*  shown  in 
Tables  2-4  was  evaluated  by  first  determining  the  conditions  under  which  the 
data  was  recorded.  By  examining  the  reports  which  describe  the  execution  of 
the  experimental  studies  [21,22],  the  data  files  were  determined  to  be  from 
either  resting  or  working;  prone,  sitting,  or  standing;  swim  suited  or  DTP 
clad;  and  dry  or  immersed  subjects  who  were  breathing  either  a  helium-oxygen 
mixture  or  air.  From  this  information,  the  ambient  temperature  and  pressure 
and  the  date  of  the  study's  execution,  the  data  files  were  divided  into 
similar  groups. 

To  minimize  the  effects  of  individual  subjects  within  a  group  ana  to  give 
a  degree  of  generality  to  our  conclusions,  it  was  decided  to  create  average 
subject  (composite)  files.  The  average  subject  data  was  generated  by 
averaging  the  data  at  each  time  increment  from  all  of  the  selected  subjects. 
This  process  was  performed  separately  for  each  segment,  thus  producing  twelve 
composite  segmental  temperature  and  twelve  composite  segmental  rate  of  heat 
toss*5  profiles. 

To  reduce  the  possibility  of  basing  our  conclusions  on  unrealistic 
information,  all  of  the  subjects'  segaentdl  data  was  pessed  through  filters 
before  being  plotted.  The  filter  applied  to  the  segmental  temperature  data 

-*See  Foreword  for  the  connotation  of  file. 

^See  Foreword  for  the  connotation  of  rate  of  heat  loss. 
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Diver  Thermal  Protection  Garment  Evaluations  [13,22] 
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Rest  =  Resting;  W/C/50  =  Working  at  a  Constant  rate  of  50  watts;  W/I/50  =  Working  Intermittently  (work  6  min 
rest  4  min.)  at  a  constant  work  rate  of  50  watts;  W/I/Var  =  Working  Intermittently  (work  6  min., 
rest  4  min.)  at  an  increasing  work  rate  of  50,  100,  150  watts;  R/LW  =  Rest/Light  Work; 

LW/MW  =  Light  Work/Moderate  Work;  LW/R  =  Light  Work  followed  by  Rest 
Sit  =  Sitting,  P  »  Prone,  S1t/45  =  Sitting  with  legs  at  45°  angle 


stipulated  that  the  segmental  temperature  (°C)  could  not  change  by  more  than 
1.0  °C  during  one  time  increment.  This  was  applied  at  each  time  increment 
after  the  first  ten.  The  filter  used  with  the  segmental  rate  of  heat  loss 
data  stipulated  that  the  segmental  change  in  the  rate  of  heat  loss  (watts) 
should  not  be  more  than  15  watts  during  one  time  increment.  This  filter  was 
used  only  after  the  fourth  time  increment.  The  delayed  use  of  these  filters 
was  to  allow  for  sudden  changes  in  the  temperatures  ana  rates  of  heat  loss  at 
entry  into  the  cold  environment. 

The  selected  subjects  (mentioned  above)  were  chosen  from  the  available 
pool  of  subjects  by  producing  a  plot  for  each  segment  which  superimposes  the 
corresponding  segmental  profiles  of  each  subject.  Examination  of  the  24 
segmental  temperature  and  rate  of  heat  loss  plots  indicated  if  any  subject 
seemed  to  be  significantly  different  from  the  rest  of  the  subjects  plotted. 

If  an  experimental  subject  appeared  to  be  consistently  different  across  all 
segments,  then  that  subject  was  excluded  from  the  segmental  subject  averages. 

Using  only  the  selected  subjects,  segmental  subject  averages,  M(i),  were 
produced  (at  each  time  increment)  with  Equation  2. 

I  Sj(1) 

N(t)  ■  iil -  (2) 

n 

where:  M(i)  =  Segmental  subject  average  temperature  or  rate  of  heat 
flow  at  time  1 

S-j(i)  =  Segmental  temperature  or  rate  heat  flow  of  subject  j  at 
time  i 

n  a  Number  of  selected  subjects 

i  =  Time  increment  index 

j  =  Subject  number  index 
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A  standard  deviation  [32],  V(i)  (Eq.  3),  was  also  calculated  for  each  segment 
at  each  time  increment  as  the  subject  average  was  being  determined. 


1/2 


V(i)  = 


l  (Sj{i))2  -  (  l  3j(1)}2 
j-!  _ 

_  n 

n  -  1 


(3) 


where:  V(i)  =  Segmental  standard  deviation  about  the  segmental  subject 
average,  M(i } 

Sj(i)  =  Segmental  temperature  or  rate  heat  flow  of  subject  j  at 
time  i 

n  =  Number  of  selected  subjects 

i  =  Time  increment  index 

j  =  Subject  number  index 


The  subject-averaged  segmental  skin  temperatures  and  rates  of  heat  loss 
:  were  then  plotted  one  segment  per  figure.  By  examining  these  mean  plots,  it 
;  was  found  that  some  of  the  subject  data  had  been  recorded  out  of  time 
sequence.  For  example,  a  typical  subject  could  have  properly  sequenced  data 
;  recorded  at  3  minute  time  increments  from  time  zero  to  27  minutes  into  the 
,  test.  Then  the  next  time  record^  might  occur  at  a  time  of  29  minutes  rather 

i 

|  than  the  expected  time  of  30  minutes.  From  this  point  on,  all  of  the  data 
could  remain  out  of  sequence  or  another  incorrect  time  increment  could  occur 
j  to  cause  the  data  to  return  to  the  initial  time  sequence.  In  either  case,  the 

I 

subjects'  mean  plots  would  show  a  sawtooth  appearance  due  to  the  time  records 
I  of  one  subject  being  out  of  sequence  with  the  other  subjects.  To  eliminate 
this  problem,  all  of  the  offset  data  was  either  arbitrarily  shifted  back  into 
sequence  or  deleted  from  the  averages.  The  choice  of  corrective  method  was 

j  ~ 

!  ^See  Foreword  for  connotation  of  time  record. 
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I 
l 

!  based  on  the  amount  of  data  out  of  synchroni zation  as  well  as  a  consideration 
I  of  the  limited  number  of  subjects  available  in  any  group. 

The  mean  (subject-averaged)  segmental  temperatures  and  rates  of  heat  loss 
i  were  then  reproduced  and  stored  in  a  new  file.  These  stored  values  are 
!  representative  of  the  data  that  would  be  expected  from  an  average  diver  having 
a  total  body  surface  area  similar  to  the  average  of  the  body  surface  areas  of 
the  selected  subjects.  Subsequently,  all  of  the  experimental  data  used  to 
|  create  and  evaluate  this  model  was  taken  from  these  average  diver  files.  The 
decision  to  use  average  cjivar  data  is  based  on  the  desire  to  have  the  model 
|  represent  a  range  of  divers  and  to  give  the  results  a  degree  of  statistical 
significance.  Also,  since  the  model  predicts  supplementary  heating  values 
based  on  average  values  of  allowable  regional  heat  flux,  the  use  of  this 
average  man  data  should  not  cause  any  significant  errors. 


Generation  of  an  Example  Average  Subject  File 

As  an  example  of  the  procedure  described  in  the  previous  section,  the 
development  of  the  average  man  file  from  the  four  subjects  defined  in  files 
i  DTP21R1,  DTP22R1,  07P23R1,  and  DTP24R1  will  be  shown.  These  experiments  were 
performed  by  NCSC  to  evaluate  the  performance  of  a  prototype  DTP  garment  worn 
by  resting,  mostly  prone  divers  in  5  °C  water  at  a  pressure  equivalent  to 
3  MSU, 

From  the  report  describing  these  dives  [22],  it  was  discovered  that 
subject  number  one's  suit  flooded  early  in  the  experiment.  Since  a  wet 
subject  is  expected  to  have  thermal  parameters  which  are  quite  different  from 
those  of  a  dry  diver,  we  decided  to  exclude  the  data  of  subject  1  (found  in 
file  DTP21R1 )  from  the  average  man  calculations.  The  remaining  subjects  were 


then  plotted  one  segment  at  a  time  on  the  same  set  of  axes.  Figures  3-14  show 
the  segmental  subject  temperature  profiles,  and  Figures  15-26  display  the 
segmental  subject  rate  of  heat  loss  plots.  Each  figure  displays  three 
segmental  profiles,  one  for  each  of  the  selected  subjects.  By  examining  tha 
twelve  temperature  and  twelve  heat  loss  figures,  it  can  be  seen  tnat,  for  the 
most  part,  the  segmental  profiles  which  were  taken  from  three  different 
subjects  are  quite  similar.  It  was  found  that  the  data  recorded  from  the 
head,  hand,  and  toot  was  generally  quite  anomalous  when  compared  to  the  other 
segments,  but,  for  consistency,  these  segments  were  averaged. 

The  average  subject  data  is  now  generated  from  the  segmental  data  of  the 
selected  subjects.  Figures  2/-3S  show  the  subject-averaged  segmental 
temperature  profiles  which  correspond  to  the  segmental  subject  temperature 
data  of  Figures  3-14,  and  Figures  39-5C  display  the  subject-averaged  segmental 
rate  of  heat  loss  profiles  derived  from  me  segmental  subject  rate  of  heat 
loss  data  of  Figures  15-26.  All  of  the  mean  subject  data  was  calculated  with 
Equation  2.  The  dashed  lines  enclosing  each  of  the  segmental  subject-averaged 
profiles  represent  t  1  standard  deviation  (Eq.  3)  about  the  mean  value  that  is 
calculated  at  each  time  increment.  The  mean  pro! lies  of  Figures  27-50  did  not 
exhibit  any  unsynchronous  times,  thus  no  corrective  measures  were  required. 

The  segmental  subject-averaged  (composite)  temperatures  (Fig.  27-38)  and  rates 
of  heat  loss  (Fig.  39-50)  are  stored  in  file  DTP2AR1  in  the  format  as 
described  in  Appendix  A.  The  data  stored  In  file  DTP2AR1  represents  a 
composite  subject  which  was  generated  frost  the  data  of  three  different 


divers. 


Figures  3-14.  Segmental  Subject  Temperature  Profiles.  Each  figure  displays 

the  three  subject's  temperature  profiles  for  a  specific  segment. 
Figures  3-14  correspond  to  segments  1-12  of  Figure  2  and 
Table  1. 
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Segmental  Subject  Temperature  Profiles  for  the  Head  (HEAD). 


Figure  4.  Segmental  Subject  Temperature  Profiles  for  the  Chest  (CHST). 


Segmental  Subject  Temperature  Profiles  for  the  Upper  Back  (UBAC). 
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Figure  7.  Segmental  Subject  Temperature  Profiles  for  the  Lower  Back  (LBAC). 


Figure  8.  Segmental  Subject  Temperature  Profiles  for  the  Arm  (ARM). 


Figure  9.  Segment  il  Subject  Temperature  Profiles  for  the  Hand  (HRST) 
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Figure  10.  Segmental  Subject  Temperature  Profiles  for  the  Front  Thigh  (F-TH). 
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Figures  15-26.  Segmental  Subject  Rate  of  Heat  Loss  Profiles.  Each  figure 

displays  the  three  subject's  rate  of  heat  loss  profiles  for  a 
specific  segment.  Figures  15-26  correspond  to  segments  1-12 


of  Figure  2  and  Table  1 
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Figure  15.  Segmental  Subject  Rate  of  Heat  Loss  Profiles  for  the  Head  (HEAD). 
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Segmental  Subject  Rate  of  Heat  Loss  Profiles  for  the  Upper  Back  (UBAC). 
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Figure  23.  Segmental  Subject  Rate  of  Heat  Loss  Profiles  for  the  Front  Calf  (FCLF). 


Figures  27-38.  Segmental  Subject-Averaged  Temperature  Profiles.  Each  figure 
displays  a  specific  segment's  mean  temperature  profile. 
Figures  27-38  *ere  developed  from  data  of  Figures  3-14  and 
correspond  to  segments  1-12  of  Figure  2  and  Table  1. 
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figure  27.  Segmental  Subject-Averaged  Temperature  Profile  for  the  Head  (HEAD). 
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Figure  28.  Segmental  Subject-Averaged  Temperature  Profile  for  the  Chest  (CHST). 


figure  2 9.  Segmental  Subject -Averaged  Temperature  Profile  for  the  Abdomen  (ABD 
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•jure  JO.  Segment <sl  Subject -Averaged  Temperature  Profile  for  the  Upper  Back  (UBAC). 


figure  32.  Segmental  Subject-Averaged  Temperature  Profile  for  the  Arm  (ARM) 
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fUjurc  33.  Segmental  Subject-Averaged  Temperature  Profile  for  the  Hand  (W?ST). 
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jed  Temperature  Profile  for  the  Front  Calf  (FCLF) 


Figure  37.  Segmental  Subject-Averaged  Temperature  Profile  for  the  Rear  Calf  (RCLF) 


Figure  28.  Segmental  Subject-Averaged  Temperature  Profile  for  the  Foot  (FOOT). 


Figures  39-50.  Segmental  Subject-Averaged  Rate  of  Heat  Loss  Profiles.  Each 
figure  displays  a  specific  segment's  mean  rate  of  heat  loss 
profile.  Figures  39-50  wei e  developed  from  data  of 
Figures  15-26  and  correspond  to  segments  1-12  of  Figure  2  and 
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Figure  39.  Segmental  Subject-Averaged  Rate  of  Keat  Loss  Profile  for  the  Head  (HEAD). 


Figure  40.  Segmental  Subject-Averaged  Rate  of  Heat  Loss  Profile  for  the  Chest  (CHST). 


Segmental  Subject-Averaged  Rate  of  Heat  Loss  Profile  for  the  Abdomen  (ABO). 
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Segmental  Subject-Averaged  Rate  of  Heat  Loss  Profile  for  the  Upper  Back  (UBAC). 


Averaged  Rate  of  Heat  Loss  Profile  for  the  Lower  Back  (LBAC). 
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Segment jl  Subject-Averaged  Rate  of  Heat  Loss  Profile  for  the  Front.  Thigh  (F-TK 


it ject -Averaged  Rate  of  Heat  Loss  Profile  for  the  Front  Calf  (FCLF). 
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gure  49.  Segmental  Subject-Averaged  Rate  of  Heat  Loss  Profile  for  the  Rear  Calf  (RCLF) 
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Description  of  Data  Groups 

Tables  2-4  list  all  of  the  data  files  we  acquired  which  had  been  recorded 
with  investigator  calibrated  instrumentation.  As  mentioned  earlier,  we 
employed  subject-averaged  (composite  subject)  data  files  in  all  of  our 
developments  to  permit  a  general  application  to  divers.  Each  of  the  data  file 
groupings,  which  represent  a  composite  subject,  are  listed  below  with  a  brief 
description  deliniating  the  files  included,  the  source  of  the  data,  any 
problems  encountered,  and  whether  or  not  we  used  the  data  (and  if  not,  why). 

Dashed  lines  denote  a  group  change.  Each  group  has  been  assigned  a  file 
for  the  storage  of  its  composite  data  even  though  some  of  the  groups  were 
eliminated  before  any  composite  subject  data  was  computed.  These  storage 
files  were  given  the  names  which  appear  in  the  first  column  (Composite  File 
Name)  of  the  tables.  These  mean  file  names  are  used  below  to  identify  the 
groups,  and  the  groups  are  listed  in  the  order  of  appearance  in  Tables  2-4. 

RHL1AR1 

The  six  subjects  in  this  group  were  exposed  to  a  20  °C  hyperbaric 
helium-oxygen  environment  at  a  pressure  of  200  MSW  while  wearing  only  swim 
suits  and  resting  in  a  sitting  position  [21].  These  subjects  participated  in 
a  respiratory  heat  loss  study  performed  by  the  Navy  Experimental  Diving  Unit. 
The  segmental  temperature  and  rate  of  heat  loss  plots,  which  compare  the 
segmental  profiles  from  each  of  the  six  subjects,  revealed  that  the  data 
recorded  from  subject  3  (RHL13R1)  was  uncharacteristic  when  compared  to  the 
five  other  subjects.  Thus  the  data  file  RHL13R1  was  excluded  from  the 
development  of  the  composite  subject.  The  mean  subject  profiles  which  are 
stored  in  RHL1AR1  were  computed  from  subjects  1,  2,  4,  5,  and  6,  whose  data  is 
stored  in  files  RH111R1,  RHL12R1,  RHL14R1,  RHL15R1,  and  RHL16R1,  respectively. 
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The  experimental  heat  flux  data  contained  in  these  files  was  recorded  in  watts 
per  meter  squared.  As  a  matter  of  accuracy,  all  of  the  segmental  heat  flux 
data  from  the  individual  subjects  was  converted  to  watts  before  beginning  the 
subject  averaging  process.  This  was  accomplished  by  multiplying  the  segmental 
data  in  W/M?  by  the  segmental  surface  area.  Each  subject's  segmental  surface 
area  was  determined  by  multiplying  his  recorded  total  body  surface  area  by  the 
respective  Hody  segmental  surface  area  fraction  (Table  1). 

The  RHL1AR1  data  was  employed  in  developing  our  segmental  skin 
temperature  predictive  correlation  for  the  twelve  Hody  segments  (Fig.  2). 

This  correlation  is  adapted  from  the  eight  region  model  proposed  by 
Kerslake  [2]  and  is  described  in  detail  in  Chapter  IV.  This  composite  data 
was  also  used  to  develop  convective  heat  transfer  coefficients  for  the 
immersed  diver  subjects.  This  development  was  performed  by  extrapolating  data 
from  one  environment  to  another  by  use  of  dimensionless  heat  transfer 
correlations. 

RHL1AR2 

The  six  subjects  in  this  grouping  were  exposed  to  a  15  °C  hyperbaric 
helium-oxygen  environment  at  a  pressure  of  200  MSW  while  wearing  only  swim 
suits  and  resting  in  a  sitting  position  [21].  These  subjects  were 
participants  in  the  same  respiratory  heat  loss  study  as  the  subjects  in  the 
RHL1AR1  grouping.  The  segmental  temperature  and  rate  of  heat  loss  plots, 
which  compare  the  segmental  profiles  from  each  of  the  six  subjects,  revealed 
no  characteristic  pattern  among  the  subjects.  The  lack  of  similarity  among 
these  subjects  suggested  that  the  data  produced  by  subject  averaging  would  be 
unrepresentative  of  the  group.  Thus  the  data  in  files  RHL12R2,  RHL13R2, 
RHL14R2,  RHL15R2,  RHL16R2,  and  RHL17R2  were  not  used  for  any  of  our 
developments. 
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RHL1A62 

The  two  divers  in  this  group  were  exposed  to  a  15  °C  hyperbaric 
helium-oxygen  environment  at  a  depth  of  200  MSW  while  clad  in  swim  suits  and 
resting  in  a  sitting  position  [21].  The  difference  between  this  set  of  dives 
and  the  RHL1AR2  dives  is  that  these  subjects  were  breathing  a  heated 

helium-oxygen  mixture  at  20  °C,  and  the  RHL1AR2  divers  were  breathing  an  ; 

I 

unheated,  15  °C  helium-oxygen  mixture.  Since  these  subjects  were  breathing  \ 

! 

warm  gas,  they  could  not  be  used  in  comparisons  with  the  other  helium-oxygen  j 
data.  Thus,  because  of  the  comparison  difficulties  and  the  fact  that  the 
15  °C  ambient  temperature  was  probably  too  cold,  we  decided  not  to  use  this 
data  in  any  of  our  developments.  < 

i 

I 

RHL1AW1 

The  five  swim  suited  subjects  of  this  group  were  exposed  to  a  20  °C 
hyperbaric  helium-oxygen  environment  at  a  depth  of  200  MSW  while  sitting  on  a 
bicylcle  ergometer  and  working  at  a  fixed  rate  of  50  watts.  These  divers  were  j 

i 

also  participants  in  the  respiratory  heat  loss  study  performed  by  NEDU.  1 

i 

The  plan  for  the  development  of  the  diver  heat  loss  model  first  called  I 

i 

for  the  development  and  verification  of  a  resting  model,  followed  by  (data  ' 

i 

permitting)  development  of  a  numerical  working  diver  model.  Due  to  the  small  j 
number  of  working  studies  performed  under  the  same  environmental  conditions  | 

and  at  the  same  rates  of  work,  we  decided  (in  consultation  with  our  Navy  j 

Medical  Research  and  Development  Command  sponsor)  to  defer  the  development  of  : 

! 

the  working  model  until  more  experimental  data  was  available.  As  a  result, 
our  analysis  of  the  working  data  has  been  limited  to  a  basic  decision  of  Its 
usefulness. 


The  data  found  in  RHL11W1,  RHL13W1,  RHL14W1 ,  RHL15W1,  and  RHL16W1  could 


easily  be  converted  to  an  average  record. 


This  composite  data  file  could  then 
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1 

be  employed  for  the  purposes  of  adapting  the  resting  segmental  predictive  skin 
temperature  correlation  and  the  convective  heat  transfer  coefficients  to  the 
working  case. 

DTP1ART 

This  mean  data  file  was  generated  from  files  DTP11RT,  DTP12RT,  DTP13RT, 
and  DTP14RT.  These  studies  were  performed  by  the  NCSC  to  evaluate  their  : 

i 

thermal  data  acquisition  system.  Each  diver  was  dressed  in  the  prototype  DTP 
garment  and  submerged  to  a  depth  equivalent  to  3  MSW  in  5  °C  water  [13,22]. 

These  air-breathing  divers  were  then  requested  to  remain  at  rest  and  in 
roughly  a  prone  position  for  the  duration  of  the  dive.  The  plots  which 

display  the  superimposed  segmental  profile  of  each  diver  showed  a  consistent  1 

j 

pattern  as  a  function  of  time  among  all  the  subjects.  Initially  the  composite  ■ 

i 

profiles  displayed  a  sawtooth  pattern  which  is  characteristic  of  i 

I 

unsynchronized  time  records6.  Examination  of  this  composite  segmental  data 
indicated  that  not  all  of  the  constituent  files  were  synchronized  with  a  three  [ 
minute  time  increment.  Since  the  number  of  time  records  which  were  not  ■ 

synchronous  was  small  compared  to  the  total  number  of  time  Increments,  the 

I 

unsynchronous  time  records  were  deleted  from  all  computations.  The  corrected  j 
rate  of  heat  loss  data  stored  In  DTP1ART  is  in  the  units  of  watts.  Since  the  | 

j 

original  data  files  were  recorded  in  KCAL/M^.HR,  the  data  was  converted  to 

i 

watts  before  beginning  the  averaging  operation.  Readily  available  conversion 
factors  and  the  individual  subject's  seumental  surface  area  were  used  to  ! 

change  KCAL/M^.HR  to  watts.  These  segmental  surface  areas  are  computed  in  the  | 
same  manner  as  described  for  the  RHL1API  file.  This  data  was  subsequently 
used  to  verify  the  heat  loss  model  for  a  prone,  resting,  DTP  clad  diver. 


6See  Foreword  for  the  connotation  of  time  record. 
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DTP2AR1 

The  four  subjects  of  this  group  were  participants  in  a  study  performed  by 
NCSC  to  evaluate  a  prototype  DTP  garment  during  air  diving  [22].  Each  subject 
was  clad  in  the  prototype  DTP  garment  and  submerged  in  5  °C  water  to  a 
pressure  equivalent  to  3  MSW.  The  divers  were  to  rest  in  a  prone  position 
until  instructed  otherwise.  From  the  information  provided  about  this  series 
of  dives  [22],  it  was  learned  that  subject  number  one's  suit  flooded  very  soon 

after  immersion.  Thus  all  data  taken  from  subject  1  and  contained  in  DTP21R1  ! 

I 

was  deleted.  The  segmental  plots,  which  superimposed  the  segmental  profile  of  ! 
each  subject,  did  not  indicate  any  irregularities  among  the  divers.  Thus  the 
three  subject's  data  that  are  contained  in  DTP22R1,  DTP23R1,  and  DTP24R1  was 
averaged  through  use  of  Equation  2.  As  was  the  case  with  the  original  source  ! 

i 

files  which  make  up  the  DTP1ART  group,  the  constituent  source  files  of  the  j 

i 

DTP2AR1  group  were  recorded  in  KCAL/M2*HR.  As  indicated  earlier,  all  heat  i 

i 

flux  data  was  converted  to  units  of  watts  before  beginning  the  averaging  ! 

process.  For  this  group,  the  applied  conversion  factors  were  created  from 
commonly  available  conversion  factors  and  the  segmental  surface  area  of  the  J 

appropriate  subject.  These  segmental  surface  areas  were  computed  from  the  | 

i 

total  body  surface  area  of  each  subject  by  using  the  Hody  segmental  surface 

j 

area  fractions  (Table  1),  as  was  the  case  in  the  RHL1AR1  section.  This  data 

i 

was  used  to  verify  the  heat  loss  model  for  a  resting,  DTP  clad,  air-breathing,  | 

| 

submerged  diver  in  water  at  a  pressure  and  temperature  of  3  MSW  and  5  °C, 
respectively.  The  data  presented  here  and  in  the  DTP1ART  file  gave  us  two 
sets  of  similarly  conditioned  studies  with  which  to  verify  our  predictive 
diver  heat  loss  model,  therefore  improving  the  statistical  authenticity  of  the 
model.  1 
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DTP2AW2 

The  four  subjects  constituting  this  group  participated  in  the  same  series 
of  garment  evaluation  studies  as  the  previous  group  and  were  exposed  to  the 
same  environmental  conditions.  The  only  difference  is  that  these  divers 
worked  on  a  vertical  ergometer  at  a  rate  of  50  watts  for  six  minutes  followed 
by  four  minutes  of  rest.  This  sequence  of  work  followed  by  rest  was  repeated 
throughout  the  experiment.  The  data  found  in  DTP21W2,  DTP22W2,  DTP23W2,  and 
DTP24W2  could  be  easily  combined  into  a  mean  (composite)  subject  file.  As  was 
the  case  with  the  corresponding  resting  studies,  the  heat  fluxes  were  recorded 
in  KCAl/M2»HR  and  must  be  converted  to  watts  before  being  processed  into 
averages.  This  data  would  be  useful  in  verifying  a  working  model  for  a 
sitting  subject  who  is  intermittently  working  at  a  constant  rate  of  50  watts. 

I 

DTP2AW3 

I 

The  four  files  in  this  group  constitute  the  final  dives  performed  in  this 
series  of  NCSC  garment  evaluation  studies  [22].  Here  the  subjects  worked 
under  the  same  immersed  conditions  as  the  previous  groups  involved  in  the 
garment  evaluations.  The  only  difference  between  this  group  and  the  DTP2AW2 

i 

grouping  is  that,  instead  of  working  at  a  fixed  rate,  the  subjects  started  at  ' 

i  i 

|  i 

j  50  watts  but  incrementally  Increased  the  ergometer  work  rate  by  50  v»tts 

!  between  successive  cycles  of  work  for  six  minutes  and  rest  for  four  minutes 

! 

until  150  watts  had  be  completed.  At  that  point,  the  subjects  ceased  working  j 
and  rested  for  the  duration  of  the  study.  Since  the  information  provided  on  ' 

files  0TP21W3,  DTP22W3,  DTP23W3,  and  DTP24W3  [22]  indicated  that  no  problems  j 

occurred  during  the  dive,  these  files  could  be  easily  applied  to  generate  an  I 

average  (composite)  record.  As  seen  with  the  other  groups  Involved  in  these 

I 

i  garment  evaluation  studies,  the  heat  loss  data  was  recorded  in  KCAl/M^.HR  and  ! 

i 

|  therefore  must  be  converted  to  watts  before  being  processed  into  averages. 
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These  studies  will  not  be  very  useful  in  developing  an  initial  working  model. 
This  statement  is  made  since  the  subjects  performed  at  a  higher  rate  of  work 
than  any  other  available  files,  and  since  these  subjects  quit  working  after 
reaching  their  maximum  rate.  This  pattern  is  unique  to  this  set  of  studies. 

DTP3AR1 

This  set  of  dives  was  performed  in  a  hyperbaric  water  environment  at  a 
depth  of  21  MSW  to  further  evaluate  the  prototype  DTP  garment.  These  files 
were  provided  to  us  by  NEDU  although  the  study  was  performed  jointly  by  NEDU 
and  NCSC  [22].  The  data  in  files  DTP33R1,  DTP34R1,  DTP35R1,  and  DTP36R1  was 
all  taken  from  air-breathing,  DTP  clad  subjects  who  were  exposed  to  5  °C  water 
at  a  hyperbaric  pressure  equivalent  to  21  MSW  and  resting  in  a  prone  position, 
j  The  segmental  plots,  which  show  the  superimposed  segmental  profile  of  each 

i 

subject,  indicated  that  the  data  of  subject  4  was  consistently  different  from 

j  the  other  subjects.  Thus  the  data  in  file  DTP34R1  was  excluded  from  the 

| 

!  production  of  the  composite  files.  When  the  segmental  subject-averaged 

i 

!  profiles  were  plotted,  a  great  number  of  unsynchroncus  times  were  discovered. 

i 

i 

|  Due  to  the  large  number  of  data  points  and  the  number  of  files  Involved,  the 

I  unsynchronized  time  records  were  realigned  by  decrementing  each  unsynchronous 

1 

» 

time  value  by  one  minute.  The  resulting  files  still  had  a  few  incongruous 

time  records  which  were  deleted.  The  subject-averaged  values  were  computed  by 

i 

using  the  corrected  files  and  are  stored  in  DTP3AR1.  As  was  the  case  with  the 
original  source  files  which  make  up  the  RHtlARl  group,  the  constituent  source 
files  of  the  DTP3AR1  group  were  recorded  in  W/M^.  Thus,  the  D1P3AR1  group's 
|  unit  conversions  were  performed  in  the  same  manner  described  for  the  RHL1AR1 

'  data.  This  data  was  used  to  verify  the  heat  loss  model  for  the  prone,  resting 

1 

■  diver  at  21  MSW  even  though  it  was  necessary  to  correct  the  time  records  to 
|  obtain  a  representative  data  file.  We  used  this  average  data  in  conjunction 


t 
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with  the  other  two  DTP  files  (DTP1ART  and  DTP2AR1)  to  enhance  the  statistical  | 

i 

significance  of  our  model  verification.  j 

DTP3AW1  ! 

- '  j 

The  five  subjects  of  this  composite  file  were  also  part  of  the  21  SW  j 

i 

garment  evaluation  study  described  in  the  DTP3AR1  section.  These  j 

air-breathing,  seated  subjects  were  exposed  to  the  5  °C  wet  hyperbaric 

i 

environment  at  a  depth  of  21  MSW.  The  difference  between  these  dives  and 

I 

those  of  DTP3AR1  is  that  this  data  was  taken  from  subjects  who  worked  at  a  j 

i 

rate  of  50  watts  for  six  minutes  and  then  rested  for  four  minutes.  These 
sitting  subjects  repeated  the  work/rest  sequence  for  the  duration  of  the 
study.  Due  to  the  large  number  of  corrections  that  had  to  be  applied  to  the 
resting  studies  (DTP3AR1)  of  this  series,  there  is  reason  to  believe  that  the 
data  found  in  files  DTP32W1,  DTP33W1,  DTP34W1,  DTP37W1,  and  DTP38W1  will  have 
many  problems.  With  this  in  mind,  the  creation  of  the  subject-averaged  file 
(DTP3AW1)  is  likely  to  be  rather  Involved.  As  described  in  the  DTP3AR1 
section,  the  heat  flux  values  will  need  to  be  converted  from  watts  per  meter 
squared  to  watts.  The  conversion  which  was  used  with  the  DTP3AR1  data  may 
also  be  applied  here.  This  data  could  be  used.  If  acceptable,  to  verify  the 
working  model  for  a  sitting  diver  who  is  intermittently  working  at  a  fixed 
rate  of  50  watts, 

DTP4AR1 

The  three  subjects  In  this  group  were  not  resting  as  Implied  by  the  name, 
but  were  performing  light  work  of  unknown  magnitude.  This  data  was  taken  from 
DTP  clad,  submerged  subjects  in  5  °C  water  at  a  depth  of  3  MSW.  Thesr 
subjects,  on  command,  would  swim  to  a  new  position  on  the  bottom  of  the  tank 
and  perform  a  timed  manual  dexterity  test.  The  studies  in  this  group  were 
performed  by  NCSC  to  evaluate  the  DTP  garment's  ability  to  satisfy  the  mission 
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of  providing  suitable  protection  for  six  hours  [31].  These  operational  tests 
were  not  designed  to  acquire  heat  transfer  data,  but  were  designed  to  evaluate 
the  response  of  the  suit  and  djver  at  the  extremes  of  the  suit 
capabilities  [31].  Thi*  data,  most  probably,  will  not  provide  statistically 
signifies4  resuits  and  would  very  difficult  to  compare  with  the  other  data 

files  since  the  rat*  of  wc  r  not  quantized. 

0TP4XYZ 

The  remaining  eight  files  are  single  case  tests  in  which  the  subjects 
performed  specific  tasks  or  were  exposed  to  a  sudden  change  in  environment, 
such  as  total  suit  flooding  or  flooding  of  a  hand.  The  fact  that  the 
protocols  used  in  taking  this  data  are,  in  general,  different  for  each  dive 
makes  it  impossible  to  group  this  data  together  and  thus  to  obtain  any 
statistical  significance.  As  explained  under  the  DTP4AR1  section,  these 
studies  were  not  performed  to  acquire  heat  transfer  data,  but  to  evaluate  the 
response  of  the  suit  and  the  diver  to  unique  situations.  These  dives  were 
done  as  a  check  of  the  DTP  garment's  ability  to  meet  its  specified  mission 
scenario  [31]. 


Summary 

This  chapter  has  described  the  procedure  used  to  create  the  mean 
(composite)  subject  files  which  were  used  in  all  verification  and  development 
phases  of  our  resting,  submerged  diver  heat  loss  model.  As  stated  in  the 
description  of  the  fi-st  working  data  file  which  is  presented  (RBliAUl),  we 
have  evaluated  the  available  working  data  fron  the  perspective  of  its 
usefulness  in  the  developoent  of  a  working  diver  ssodel.  We  concluded  fron 
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i  this  evaluation,  with  the  concurrence  of  our  Navy  contract  sponsor,  that  the 
adaptation  of  our  resting  diver  heat  loss  model  to  the  case  of  a  working  diver 
;  should  be  deferred  until  adequate  experimental  working  data  is  available.  As 

i 

i  Indicated  previously,  this  conclusion  is  based  on  the  fact  that  there  are  only 

| 

;  two  available  working  diver  data  sets  that  were  performed  in  similar 
i  environments  and  at  the  same  work  rate  and  cycle  (rest  to  work).  However, 

l 

I  these  dives  were  conducted  at  two  different  depths  (3  MSW  in  DTP2AW2  and 
;  21  MSW  in  DTP3AW1),  thus  limiting  the  number  of  mean  (composite)  data  files 

I 

;  each  depth  to  one.  The  lack  of  statistical  significance,  when  using  su^  a 
,  small  amount  of  data,  would  make  any  findings  tenuous. 

The  four  composite  data  files  (RHL1AR1,  DTP1ART,  DTP2AR1,  and  TABARI) 
are  tabulated  in  Appendix  8.  The  data  files  used  to  create  these  files  and 
the  other  data  of  Tables  2-4  may  be  obtained  from  the  Director  of  the 
F.G.  Hall  Environmental  Laboratory,  Duke  University,  only  with  the  permission 
of  the  Naval  Research  and  Development  Command  and  the  originating  command. 


H 


T 

4* 


i 

i 

CHAPTER  III  j 

1 

DEVELOPMENT  OF  THE  DIVER  MODEL  ! 

\ 

l 

|  Our  development  of  the  diver  heat  loss  model  was  performed  in  three  ! 

distinct  steps.  The  first  task  was  to  develop  a  geometrical  representation  of 
an  average  diver.  The  second  objective  was  to  derive  a  mathematical 
correlation  for  predicting  the  regional  skin  temperatures  of  our  modeled  I 

subject.  The  final  step  was  to  develop  a  model  for  predicting  the  behavior  of 

i 

an  assumed  garment  ensemble  at  any  location  on  a  prone,  sitting,  or  standing 
diver.  The  development  of  the  physical  model  and  the  definition  of  the 
problem  and  its  solution  are  described  in  this  chapter.  The  development  of 
the  regional  predictive  skin  temperature  correlation  is  described  in 
Chapter  IV,  and  the  modeling  of  the  physical  behavior  of  the  garment  ensemble 
is  presented  in  Chapter  V. 

The  Physical  Model 

j  The  model  man,  as  shown  in  Figure  1,  is  fabricated  from  cylinders  and 
j  spheres.  Since  both  the  head  and  the  hand  (when  closed  to  conserve  heat) 
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resemble  spheres,  they  were  treated  as  such;  the  seven  remaining  regions  were  | 

] 

modeled  as  cylinders.  The  proportions  of  each  region  were  chosen  to  i 

! 

approximate  the  physical  dimensions  of  a  standard  U.S.  Air  Force  man  as  j 

i 

defined  by  Webb  [23]  and  to  produce  regional  surface  area  fractions  comparable  ! 
to  those  defined  for  the  twelve  segment  Hody  man  in  Table  1  and  Figure  2. 

The  choice  of  the  standard  U.S.  Air  Force  man  [23]  as  the  basis  of  the 
physical  model  was,  in  part,  a  result  of  the  lack  of  complete  anthropometric 

j 

data  found  in  the  open  literature.  Webb's  data  [23],  the  only  complete  source  | 

r 

found,  was  derived  from  a  large  population  of  U.S.  Air  Force  men.  We  assumed  ! 

1 

that  average  Navy  divers  would  conform  to  this  Air  Force  data  because  all 
military  personnel  would  be  expected  to  be  in  a  similar  state  of  physical 

i  conditioning.  The  large  number  of  subjects  represented  in  Webb's  standard  man 

\ 

data  suggested  that  accurate  predictions  would  result  from  its  use. 

The  dimensions  chosen  for  the  model  man  ere  displayed  in  Table  5.  The  : 
j  height  and  total  surface  area,  calculated  to  be  172.7  CM  and  1.902  M^, 

I 

|  respectiv  '"■»  ompare  favorably  to  the  standard  Air  Force  man's  height  of 
175,5  CM  and  body  surface  area  of  1.90  The  model  man's  regional  surface 

i 

|  area  percentages  are  compared  to  those  derived  from  riody's  segmental  surface 

j 

|  fractions  (Table  1,  In  Table  6.  Requiring  our  model  man's  dimensions  to  yield 

l 

j  regional  surface  area  fractions  comparable  to  those  of  Kody's  facilitated  use 
;  of  the  available  experimental  data. 

i 

i 

i 

\ 

s 

\ 

[  The  Kode1  and  Its  Solution 

Our  mathematical  model  for  predicting  diver  regional  heat  loss  and  any 
?  required  regional  supplementary  heating  for  a  resting  diver  who  is  wearing 
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Table  5 

Model  Man  Dimensions 


Region 

ID  # 

Region 

Name 

Shape 

Diameter 

(M) 

Characteristic 
Length  (M) 

1 

Head 

sphere 

0.206 

0.206 

2 

Torso 

cylinder 

0.300 

0.353 

3 

Abdomen 

cyl inder 

0.300 

0.353 

4 

Thigh 

cyl inder 

0.150 

0.383 

5 

Calf 

cyl inder 

0.120 

0.328 

6 

Foot 

cylinder 

0.150 

0.104 

7 

Upper  Arm 

cyl inder 

0.096 

0.223 

8 

Lower  Arm 

cyl inder 

0.070 

0.267 

9 

Hand 

sphere 

0.123 

0.123 

Total  Body  Surface  Area  {Sum  of  Regional  Surface  Areas) 


i 


I 

i 

i 

i 


1 


Surface  Area 

m 

0.1333 

0.3327 

0.3327 

0.3610 

0.2474 

0.1334 

0.1490 

0.1174 

0.0950 

1.902 


t Model  man  regional  surface  area  percents  were  calculated  by:  Regional  surface  aiea  (Table  5)  times  100.0 
divided  by  the  total  body  surface  area  (TBSA).  The  value  of  TBS A  was  1.902  M2  (Table  5). 

ttHo dy  man  regfonal  surface  area  percents  were  computed  by:  Summing  Hody  segment  surface  area  fractions 
(Table  1)  and  then  multiplying  by  100. 
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a  dry  thermal  protection  garment  is  formulated  as  a  function  of  the  following 
parameters: 

a.  Prescribed  mean  skin  temperature.  Regional  skin  temperatures  are 

l 

then  calculated  by  means  of  specific  predictive  correlations  which 
use  this  mean  skin  temperature  as  the  independent  argument. 

b.  Ambient  water  temperature. 

c.  Depth  of  immersion  of  the  diver's  shoulder. 

d.  Garment  thermal  resistance  characteristics. 

e.  Gas  with  which  suit  is  inflated  to  compensate  for  hydrostatic 
pressure. 

f.  Hydrostatic  suit  squeeze  due  to  local  hydrostatic  pressure, 

i.e.,  hydrostatic  pressure  relative  to  the  locus  of  lung  inflation. 

I 

For  ease  of  modeling,  this  is  arbitrarily  taken  as  the  height  of  the 
j  shoulders  of  an  erect  diver  (Fig.  1). 

j  g.  Posture:  prone,  sitting,  or  standing. 

j  h.  Allowable  regional  surface  heat  fluxes  for  resting  man  according  to 

i 

j  8urri$s  [33,34]. 

I 

I 

i.  Convective  film  coefficients  calculated  according  to  the  empirical 
correlations  which  are  discussed  in  Chapter  V. 

I 

j 

|  The  numerical  model  synthesizes  the  parameters  listed  above  to  obtain  an 

|  overall  thermal  conductance,  U^(j)  (Eq.  4)  for  each  node  j  of  region  i 

i 

i  (Fig.  51),  which  is  employed  in  Equation  5  for  estimating  the  actual  heat  loss 
I  per  unit  length  from  each  region  were  there  no  supplementary  heating 
available. 


Figure 


where: 


Mid-point  of 

aNode,A2/2 


51.  Vertical  Cylindrical  Region  Divided  into  Five  Nodes  (j=5). 

ui(j)  3  l  Mo)  (i 

kml  1 

Oi(O)  -  IH(J)S1CJKT1  -  ra)  {5 

U^j)  »  Overall  heat  transfer  coefficient  for  node  j  of  region  i 
Rk(j)  <*  Thermal  resistance  of  garment  layer  k  surrounding  node  j 
n  -  Kwber  of  garment  layers  covering  all  nodes  of  region  i 

6i(j)  3  Rate  of  heat  flow  from  node  j  of  region  i 

Sf(j)  a  S*:in  surface  area  enclosed  by  node  j  of  region  i 

T-j  -  Regional  predicted  skin  temperature 

Ta  *  Ambient  temperature 
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The  regional  predicted  skin  temperature,  (Eq.  5),  is  predicted  as  a 
function  of  the  prescribed  mean  skin  temperature  by  the  predictive  skin 
temperature  correlation  described  in  Chapter  IV.  The  overall  heat  transfer 


coefficient  (U^ ( j ) ) »  based  on  the  skin  surface  area  ( S-j ( j ) ) ,  is  represented 
for  each  jth  nodal  element  of  the  ith  cylindrical  region  by: 

(6) 

1 

U,(j)  =  - 

1  +  s,(j)  «  In  (r2/r1)  +  S^j)  •  In  (r3/r2)  +  $j  ( j ) 

h^j)  2nk1g  .  A2  2nkog  .  A2  s0(j)h0(o) 


I 


i 


I 


I 


I 


I 


1 


where:  Ui(j)  =  Overall  heat  transfer  coefficient  for  node  j  of  region  i 

h-j ( j )  =  Convective  heat  transfer  coefficient  for  the  gas 

layer  between  the  skin  surface  and  the  undergarment  of 
node  j  and  region  i 


Si(j)  *  Surface  area  of  skin  surface  encircling  node  j  of 
region  i 

rt  =  Radius  from  regional  center  line  to  skin  surface 


r2  a  Radius  from  regional  center  line  to  undergarment/outer 
garment  interface 

r3  *  Radius  from  regional  center  line  to  outside  surface  of 
outer  garment 

A2  «  Nodal  length 


k3g  =  Thermal  conductivity  of  inner  garment 

kog  3  Thermal  conductivity  of  outer  garment 

Sq(J)  *  Surface  area  of  outer  garment  encircling  node  j  of 
region  i 

h0(j)  *  Convective  heat  transfer  coefficient  for  the  outer 

boundary  layer  on  the  surface  of  the  outer  garment  of 
node  j  and  region  i 


The  first  thermal  resistance  term  (l/h^Cj))  accounts  for  the  insulative  effect 
that  might  be  provided  by  a  gas  layer  between  the  skin  surface  and  an 
;  insulating  garment.  This  would  be  considered  a  natural  convection  phenomenon 
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due  to  a  microcirculation.  However,  the  effect  of  this  term  is  negligible 
since  the  assumed  protective  dry  garment  is  subjected  to  hydrostatic  suit 
squeeze  below  the  chest  level  (In  the  case  of  our  erect  or  sitting  model  from 
the  shoulder  crest  down),  thus  minimizing  this  gas  layer.  The  second  term  of 
Equation  6  may  be  replaced  by  an  experimental  value  of  garment  thermal 

] 

resistance  that  has  been  extrapolated  to  the  ambient  conditions  of  the  suit 
i  entrapped  gas  and  then  further  reduced  for  the  local  effect  of  suit 

i 

i 

I  compression  (See  Chapter  V  for  details.).  The  third  term  of  the  denominator 
i  considers  the  contribution  of  the  outer  garment  to  the  overall  resistance. 

The  assumed  outer  garment  is  a  pre-crushed  neoprene  foam  of  0.16  CM 

(1/16  inch)  nominal  thickness,  worn  primarily  as  a  means  of  maintaining  water 

tight  Integrity.  The  garment  can  be  characterized  as  a  solid  sheet  of 
;  neoprene  rubber  of  comparable  thickness  because  of  the  effects  of  the 
pre-immersion  treatment  on  the  foam  material.  The  last  thermal  resistance 
term  of  Equation  6  is  computed  from  dimensionless  empirical  convective  heat 

i 

transfer  correlations  that  are  defined  in  Chapter  V. 

•  Qj(j),  as  defined  in  Equation  5,  is  calculated  for  discrete  nodes  of  a 

body  region  (Fig  51).  These  nodal  rate  of  heat  loss  values  are  summed  to 

obtain  an  overall  regional  rate  of  heat  loss  for  the  ambient  environment  under 
!  consideration.  This  computation  models  the  rate  of  heat  flow  from  the  passive 
garment  alone.  The  regional  supplementary  heating  required  to  sustain  a 
resting  diver  is  estimated  by  subtracting  a  permissible  regional  rate  of  heat 
loss,  as  defined  by  Burriss  [33,34],  from  the  value  calculated  in  the  manner 
described  above.  These  resultant  values  specify  the  rate  of  supplementary 
heating  required  for  each  region,  the  diver's  overall  heating  requirement 
(obtained  by  summing  the  heating  requirements  of  each  region),  and  the 
relative  heating  intensity  of  each  body  region  with  respect  to  the  total 
heating  requirement. 
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The  describe  predictive  model  has  been  developed  by  assuming  that  the 
diver  would  wear  some  form  of  pressure  compensated  dry  suit.  Since  all  of  the 
available  experimental  data  was  collected  from  immersed  subjects  clad  in  the 
Diver  Thermal  Protection  (DTP)  garment  [22]  whose  thermal  characteristics  are 
known  [12,13],  we  decided  to  develop  our  model  with  the  DTP  garment  as  the 
assumed  dry  suit  design.  Since  the  development  of  an  overall  heat  transfer 
coefficient  is  dependent  on  the  selected  dry  suit  design,  consideration  of  a 
suit  which  has  characteristics  different  from  those  of  the  DTP  garment 
requires  that  the  resistance  terms  of  Equation  6  be  re-analyzed.  This  may  be 
easily  performed  if  the  thermal  characteristics  of  the  garment  at  sea  level 

i  are  known  and  the  variation  of  these  parameters  can  be  modeled  as  a  function 

i 

!  of  temperature  and  pressure. 

I 

The  DTP  garment  Is  configured  as  two  different  layers  [22].  The 

i 

j  Innermost  layer  consists  of  a  thin,  cotton/wool  long  john  underwear  with  wool 

j 

!  socks  and  a  M-400  Thlnsulate7  thermal  undergarment  which  incorporates  double 

l 

f 

j  layer  boots  [22].  Because  of  the  significantly  larger  thermal  resistance  of 

[  the  Thinsulate  undergarment  relative  to  the  long  John's  [13],  the  long  John 

;  underwear  was  neglected  during  the  development  of  the  overall  heat  transfer 
coefficient.  The  outer  layer  is  the  crushed  neoprene  foam  ganaent. 


^H-400  Thinsrlate  is  a  commercial  trade  name  for  a  fibrous  batt  material 
■ade  by  the  3-H  Corporation. 


CHAPTER  IV 


DEVELOPMENT  OF  PREDICTIVE  SKIN  TEMPERATURE  CORRELATION 


It  Is  a  well  established  observation  that  the  preferred  mean  weighted 
skin  temperature  of  man  at  rest  in  a  euthermic  environment  Is  33  °C  ±  1  °C. 
Associated  with  this  resting  state  is  a  metabolic  rate  of  about  55  W/M^,  the 
majority  of  which  is  emitted  as  waste  heat  and  the  remainder  of  which  is 
utilized  for  metabolic  processes  [33-36].  Physiological  laboratory  studies 
have  defined  specific  regional  skin  temperatures  and  rates  of  heat  loss  that 
are  associated  with  this  preferred  mean  skin  temperature  [33,34]. 

Kerslake  [2]  developed  an  algebraic  model  for  predicting  the  temperature 
of  any  of  eight  different  regions®,  as  a  function  of  mean  skin  temperature, 
from  extensive  experimental  studies  of  men  in  a  nearly  steady  thermal  state. 
The  equations  that  he  developed  were  of  the  form: 

*  Tska  t  Ai  •  R  (7) 

where:  T|  *  Skin  temperature  for  region  i 
Tska  »  Mean  skin  temperature 

Aj  a  Distribution  fraction  for  region  i  (Table  7) 

ft  =  Temperature  difference  between  trunk,  and  hand  and  foot 

®$ee  Foreword  for  the  connotation  of  region  and  Kerslake  region. 
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Kerslake's  Curve  Fitting  Coefficients 


Region 
Index,  i 

Region 

Name 

^i 

Corresponding 
Indices  of 

1 

Head 

0.27 

1 

2 

Trunk 

0.27 

2,  3,  h,  i 

3 

Lower  Arm 

-0.365 

6 

4 

Hand 

-0.73 

7 

5 

Thigh 

0.00 

8  and  10 

6 

Calf 

-0.365 

9  and  11 

7 

Foot 

-0.73 

12 

8 

Upper  Arm 

0.00 

6 

■•■Data  adapted  from  Kerslake  [2] 
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Piironen  and  Takalo  [37]  conducted  a  series  of  physiological  experiments 
j  with  an  instrumented  suit  calorimeter  by  which  they  determined  both  regional 
heat  fluxes  and  skin  temperatures  for  a  variety  of  work  rates  and  degrees  of 

i 

!  vasoconstriction.  In  their  three  experimental  subjects,  they  found  that 
j  maximum  vasoconstriction  occurred  in  the  mean  skin  temperature  range  of 
|  31-32  °C  and  that  a  10%  increase  in  metabolic  heat  production  was  observed 
|  whenever  the  minimum  tissue  conductance  (maximum  vasoconstriction)  was 

I  attained.  This  observed  increase  in  metabolic  rate  did  not  correlate  with  the 

I 

|  onset  of  shivering  thermogenesis.  Shivering  appeared  only  with  a  further 
'  lowering  of  the  mean  skin  temperature  below  that  required  to  attain  the 
;  minimum  tissue  conductance.  When  shivering  did  occur,  an  associated  increase 
in  tissue  conductance  was  indicated  by  the  instrumented  suit.  Plots  of  the 
total  body  and  regional  experimental  conductances  from  three  resting  subjects 
indicated  that  the  minimum  tissue  conductance  occurs  at  a  mean  skin 

; 

temperature  of  31.5  °C  [37]. 

Ir.  designing  an  actively  heated  diver  thermal  protection  garment  of  a 
limited  energy  supply,  it  would  seem  reasonable  to  allow  peripheral 
vasoconstriction,  tut  to  avoid  the  onset  of  shivering  and  the  associated 
increase  in  tissue  conductance,  limits  of  thermal  equilibrium  are  usually 
specified  in  terms  of  a  minimum  mean  skin  temperature  and  minimum  temperatures 
for  the  hand  and  foot  [11].  Supplementary  heating  of  a  diver  may  be  required 
to  .maintain  a  specified  minimum  mean  skin  temperature,  but  most  probably,  at 
different  regional  rates  to  avoid  compromising  the  regional  skin  temperatures. 
Therefore,  a  correlation  for  predicting  regional  skin  temperature  in  the 
presence  of  peripheral  vasoconstriction  and  as  a  function  of  mean  skin 
temperature  is  most  desirable  for  specifying  minimus  regional  supplementary 


heating  requirements. 
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Development  of  Equations 

Experimental  data  obtained  from  five  nearly  nude,  sitting,  resting  divers 

I  exposed  to  a  20  °C  helium-oxygen  gas  mixture  at  a  pressure  of  200  MSW  was 

I  utilized  to  develop  a  Kerslake  type  algebraic  equation  for  predicting  Hody's 

i 

1  twelve  segmental  temperatures  (Fig.  2)  rather  than  Kerslake's  eight.  This 

'  data  was  collected  during  a  series  of  experimental  dives  conducted  by  the  Navy 
Experimental  Diving  Unit  to  study  respiratory  heat  loss  under  controlled 

| 

|  hyperbaric  conditions  [21].  As  described  in  Chapter  II,  the  segmental  data  ot 

!  these  subjects  was  averaged  to  produce  the  mean  subject  file  RHL1AR1.  Any 

i 

'  further  references  to  experimental  data  refer  to  data  contained  in  the  RHL1AR1 
composite  diver  file.  Hody's  [24,25]  twelve  point  mean  skin  temperature 
formula.  Equation  8,  was  used  in  developing  the  predictive  segmental  skin 
temperature  equations  because  all  of  the  experimental  data  available  for 
analysis  had  been  recorded  in  the  Kody  format.  Figure  2  illustrates  these 
twelve  sites  and  Table  i  defines  the  indices,  nomenclature,  abbreviations,  and 
surface  area  fractions  f  each  of  the  twelve  secpnents.  The  subscripts  of 
Equation  8  correspond  to  the  segmental  indices  appearing  in  Figure  2  and 
Table  I. 

Tskta  ■  Q.OTQTj  ♦  0.085T,  ♦  0.085T3  *  0.090T4  ♦  0.09QT5  ♦  Q.idOTg 

(8) 

*  0,0 SOT 7  *■  0.095Tg  *  Q.065Tg  ♦  0«095Tjq  ♦  0. 065T|j  ♦  0. 070Tj? 

Figures  52-63  each  display  a  segment's  te«perature  profile  from  the 
RHUARi  data  and  the  Hody  *ean  skin  temperature  calculated  with  Equation  8. 
Each  graph  is  labeled  with  the  abbreviated  se$»ent  naae  {Table  \)  that 
corresponds  to  the  segment  plotted;  the  abbreviation  KMST  indicates  the  Kody 
«ean  skin  temperature.  Each  experiaental  curve  was  analyzed  by  the  method  of 
Least  Squares  [3S],  after  incongruous  data  points  were  eliminated,  to  obtain 


a  linear  regression  equation  for  the  time  span  of  35  to  90  minutes.  This  was 
selected  as  a  representative  quasi-steady  time  period.  Figure  64  is  a 
representative  sample  of  the  generated  linear  regression  curves,  which 
displays  the  linear  regression  lines  for  the  abdomen  (ABD)  and  the  Hody  mean 
skin  temperature  (HMS1).  Subsequently,  these  linear  regression 
representations  of  the  experimental  data  were  employed  to  facilitate  the 
numerical  development  of  the  predictive  segmental  skin  temperature 
correlations. 

Our  first  task  was  to  evaluate  the  suitability  of  Kerslake's  equation  for 
predicting  regional  skin  temperatuares  under  hyperbaric  conditions.  To  obtain 
values  uf  experimental  temperature  which  correspond  to  the  eight  Kerslake 
regions,  segmental  values  at  the  84th  minute  of  Figures  52-63  were  weight 
averaged  by  using  the  Hody  surface  area  fractions  (Table  1)  of  the 
corresponding  Hody  segments.  Table  8  illustrates  the  very  good  comparison 
between  the  respective  regional  weight  averaged  experimental  temperatures  and 
the  temperatures  that  are  predicted  when  the  experimental  Hody  mean  skin 
temperature  is  used  as  the  argument  in  Equation  ?,  Three  separate  R  values, 
defined  in  Table  3,  were  employed  in  Equation  7  to  evaluate  the  accuracy  of 
the  Kerslake  equation  for  predicting  segmental  skin  temperatures  of  resting 
man  under  hyperbaric  conditions.  Table  8  also  defines  the  weight  averaging 
correlations  for  extrapolating  twelve  site  data  to  eight  sites. 

Table  9,  similar  in  fora  to  Table  3,  was  prepared  with  data  that  was 
gathered  fro®  the  linear  regression  equations  which  were  fitted  to  the  mean 
experimental  data  of  five  subjects.  The  analyses  reflected  in  both  Tables  ? 
and  9  were  conducted  with  data  obtained  at  the  8dth  minute  when  all  of  the 
subjects  were  deemed  to  have  attained  a  quasi -steady  state  of  thermal 
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Figures  52-63.  Representative  Segmental  Temperature  and  Hody  Mean  Skin 

Temperature  (HM3T)  Profiles.  One  figure  for  each  of  the  12 
segments.  Mean  data  taken  from  200  MSW  respiratory  heat  loss 
studies  (RHL1AR1) . 
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Figure  52.  Segmental  Temperature  Profile  for  the  Head  (HEAD)  and  the  Hody  Mean  Skin  Temperature  (HMST) 
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Segmental  Temperature  Profile  for  the  Upper  Back  (UBAC)  and  the  Hody  Mean  Skin  Temperat 


0.80  20.09  40.09  60.09  80.09  108.90  128.98 

MINUTE 5 

Figure  56.  Segmental  Temperature  Profile  for  the  Lower  Back  (LBAC)  and  the  Hody  Mean  Skin  Temperature  (HMST). 
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Figure  57,  Segmental  Temperature  Profile  for  the  Arm  (ARM)  and  the  Hody  Mean  Skin  Temperature  (HMST) 


0.90  20.00  40.00  60.00  80.00  109.09  129.08 

MINUTES 

Figure  5fl.  Segmental  Temperature  Profile  for  the  Hand  (WRST)  and  the  Hody  Mean  Skin  Temperature  (HMST) 
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Segmental  Temperature  Profile  for  the  Rear  Thigh  (R-TH)  and  the  Hody  Mean  Skin  Temperature 
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Figure  64.  Representative  Linearized  Segmental  Temperature  for  the  todomen 
and  Hody  Mean  Skin  Temperature  Profiles.  Linearized  data  for  the 
abdomen  taken  from  mean  data  of  Figure  54  (RK1AR1). 
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Table  8 

Comparison  of  Composite  Experimental  Temperatures  Extracted  From 
Figures  52-63  With  Those  Predicted  With  Kerslake's  Equation  (Eq.  7).* 
Experimental  data  recorded  in  the  twelve  segment  format  of  Hody  was 
extrapolated  to  Kerslake's  format  by  the  weight  averaging 
correlations  presented  below. 


Temperatures  in  °C 

Experimental  Predicted  with  Rl,  R2,  and  R3 


Region 

Rl 

R2 

R3 

Head 

26.67 

26.60 

26.45 

26.76 

Trunk 

26.57 

26.60 

26.45 

26-76 

Lower  Arm 

26.32 

24.41 

24.61 

24.20 

Hand 

23.67 

23.15 

23.56 

22.74 

Thigh 

25.65 

25.67 

25.67 

25.67 

Calf 

24.43 

24.41 

24.61 

24.20 

Foot 

22.55 

23.15 

23.56 

22.74 

Upper  Arm** 

26.32 

25.67 

25.67 

25.67 

|  R  Values  Employed  in  Equation  1: 

|  Rl  =  Trunk  -  (Hand  +  Foot)/2.0  =  3.46 

I 

R2  =  Trunk  -  Hand  =  2.90 
R3  3  Trunk  -  Foot  =>  4.02 

i  Weight  averaging  correlations: 

j 

Trunk  »  (0.085  CHST  +  0.085  ABD  ♦  0.090  UBAC  ♦  0.090  LBAC)/0. 350 
;  Thigh  »  (0.095  F-Th  +  0.095  R-TH}/0.190 

!  Calf  3  (0.065  FCLF  +  0.065  RCLF)/0.130 


*Tbe  comparison  is  predicated  on  the  mean  of  data  collected  during  the  84th 
minute  of  the  experimental  series.  The  Hody  mean  skin  temperature  (Tskm  of 
Equation  7)  was  calculated  by  means  of  Equation  8  as  25.67  °C. 

^Experimental  data  recorded  from  the  lower  arm  was  util .zed  since  none  was 
recorded  for  the  upper  arm. 


Comparison  of  Composite  Experimental  Temperatures  Estimated  by 
Means  of  Respective  Linear  Regression  Equations  Which  Were 
Derived  From  Figures  52-63  With  Those  Predicted  With  Kerslake's 
Equation  (Eq.  7).*  Experimental  values  estimated  from  Hody's  twelve 
segment  format  were  extrapolated  to  Kerslake's  format  by  the 
weight  averaging  correlations  presented  below. 

Temperatures  in  °C 


Experimental 

Predicted 

with  Rl, 

R2,  and  R3 

Region 

R1 

R2 

R3 

Head 

26.65 

26.53 

26.39 

26.68 

Trunk 

26.51 

26.53 

26.39 

26.68 

Lower  Arm 

26.10 

24.32 

24.51 

24.12 

Hand 

23.55 

23.04 

23.43 

22.66 

Thigh 

25.65 

25.59 

25.59 

25.59 

Calf 

24.37 

24.32 

24.51 

24.12 

Foot 

22.49 

23.04 

23.43 

22.66 

Upper  Arm1"1" 

26.10 

25.59 

25.59 

25.59 

R  Values  Employed  in  Equation  1: 

R1  »  Trunk  -  (Hand  *  Foot)/2.0  =  3.49 
R2  =  Trunk  -  Hand  =  2.96 
R3  =  Trunk  -  Foot  =  4.02 


Weight  Averaging  Correlations: 

Trunk  =  (0.085  CNST  +  0.085  ABD  +  0.090  U8AC  +  0.090  LBAC)/0.350 

Thigh  =  (0.095  F-TH  +  0.095  R-TH)/0.190 

Calf  a  (0.065  FCLF  +  0.065  RCLF)/0. 130 

*Tbe  comparison  is  developed  for  the  84th  minute  of  the  experimental  series. 

The  Hody  mean  skin  temperature  (Tskm  of  Equation  7)  was  calculated  by  means 
of  Equation  8  as  25.59  °C. 

^linearized  experimental  data  of  the  lower  arm  was  utilized  since  no 
experimental  data  was  recorded  for  the  upper  arm. 
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|  The  excellent  comparison  between  the  tabulated  regional  values  in  the 

j 

|  first  column  of  Tables  8  and  9  supports  our  assumption  that  the  linear 
!  regression  equations  accurately  represent  the  experimental  profiles  for  the 
span  of  35  to  90  minutes.  Also,  we  observed  the  best  comparison  between  the 
j  experimental  and  predicted  values  of  Table  8  when  the  constant  R1  was  employed 

j  in  Equation  7  to  predict  the  Kerslake  regional  temperatures. 

| 

Eight  graphs  of  the  absolute  temperature  difference  between  the  averaged 

i 

experimental  temperature  and  that  predicted  by  Kerslake' s  correlation. 

Equation  7,  were  plotted  by  using  the  respective  linear  regression  equations 
■  t0  determine  the  segmental  skin  temperatures  as  a  function  of  time.  The  value 

j 

of  mean  skin  temperature  employed  in  Equation  7  was  determined  by  means  of 
=  Equation  s  for  each  minute  of  time  from  the  experimental  temperatures  which 
were  estimated  by  the  linear  regression  equations.  R  in  Equation  7  was 
1  defined  as  being  equivalent  to  R1  of  Tables  8  and  9  and  was  calculated  at  each 

I 

time  interval  prior  to  its  use  in  Equation  7.  The  temperature  averaging 
scheme  which  was  used  to  obtain  the  equivalent  Kerslake  regional  values  from 
,  the  linearized  experimental  data  was  the  same  as  that  defined  in  Tables  8 
i  and  9.  Review  of  these  graphs  revealed  that  for  the  entire  period  of  analysis 
:  (35  to  90  minutes)  the  experimental  temperature  of  the  trunk,  thigh,  and  calf 
agreed  to  within  ±  0.5  °C  of  that  predicted  by  Equation  7.  The  same  agreement 
was  observed  for  the  head  and  hand  for  90S  of  the-  time  span,  and  for  70S  of 

i 

the  period  for  the  foot.  The  Kerslake  predicted  upper  and  lower  arm  values 
;  were  compared  to  the  linearized  experimental  data  of  the  arm  (See  Table  7.). 

f 

The  absolute  temperature  difference  between  the  predicted  value  for  the  upper 
arm  and  the  experimental  value  for  the  arm  varied  between  0.5  and  0.8  °C 
throughout  the  time  span;  for  the  lower  arm,  the  absolute  temperature 
difference  observed  was  considerably  larger,  i.e.,  greater  than  1.50  °C 
throughout  the  entire  time  span.  We  do  not  understand  the  exact  cause  of  the 


poor  comparison  between  the  predicted  and  linearized  experimental  skin 
!  temperatures  for  the  upper  and  lower  arm  when  such  a  good  correspondence  was 
observed  for  the  other  body  regions.  We  suspect  that  the  lack  of  similar 
regions  between  which  to  compare  may  have  contributed  to  the  noted 
incongruity.  Regardless  of  the  results  of  the  upper  and  lower  arm  comparison, 
the  good  agreement  noted  for  the  principal  regions  confirms  the  hypothesis  of 
extrapolating  Kerslake's  predictive  equation  for  eight  regions  (Table  7)  to 
the  twelve  segment  mao  of  Body  (Fig.  ?.  and  Table  1). 

To  adapt  Kerslake's  algebraic  model  to  the  Hody  twelve  segment  map,  it 
was  necessary  to  compute  curve  fitting  coefficients,  A(i),  from  Equation  9, 
which  is  a  transposition  of  Equation  7. 

AO)  .W)  -™">  (9) 

R1 

where:  A(i)  =  Distribution  fraction  for  segment  i 

T(i)  »  linearized  experimental  temperature  for  segment  i 
Tskm  =  linearized  Experimental  Hody  mean  skin  temperature 
R1  -  2.946 

We  substituted  into  Equation  9  values  of  T(1),  Tskm,  and  Ri,  all  of  which  had 
been  evaluated  at  the  midpoint  of  the  selected  time  span  of  35  to  90  minutes. 
We  then  solved  for  the  respective  values  of  A(i)  which  are  tabulated  in 
Table  10. 

Although  the  numerical  analysis  was  based  upon  data  which  was  obtained 
from  shivering  men  with  mean  skin  temperatures  (Tskm)  of  approximately  26  4C, 
we  propose  that  the  modified  Kerslake  correlation.  Equation  10,  can  be  used 
with  the  tabulated  values  of  Table  10  to  estimate  the  twelve  segmental  skin 
temperatures  of  a  vasoconstricted  but  non-shivering  diver  whose  mean  skin 
temperature  is  between  31  an  32  °C. 
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I 

Table  10 

Distribution  Fractions  for  Predictive  Skin  Temperature  Correlations 

I  Time  =  63  minutes 

I  R1  =  2.946 

I 

J 

I 


Segment  Segment  Distribution  Fraction 


I 


ID  #,  i  Name  _ A(i ) 

1  Head  0.36 

2  Chest  0.27 

3  Abdomen  0.30 

4  Upper  Back  0.52 

5  Lower  Back  -0.08 

6  Arm  0.22 

7  Hand  -0.64 

8  Front  Thigh  -0.04 

9  Front  Calf  -0.16 

10  Rear  Thigh  -0.06 

11  Rear  Calf  -0.55 

12  Foot  -0.86 


J 

t 


i 
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T(i)  =  Tskm  ♦  A(i )  .  R1  (10) 

J 

where:  T(i)  =  Predicted  temperature  for  segment  i  i 

Tskm  =  Hody  mean  skin  temperature  J 

I 

A(i)  =  Distribution  fraction  for  segment  i  (Table  10) 

R1  =  2.946 

The  application  of  Equation  10  to  the  Hody  mean  skin  temperature  (Tskm) 
range  of  31  to  32  °C  is  by  extrapolation.  Ideally,  this  proposed  application 
should  be  tested  experimentally  under  physiological  conditions  analogous  to  ! 
those  for  which  it  is  intended.  Realistically,  no  such  experimental  data  is  j 
currently  available.  i 


Muthenti cation  of  the  Model 

The  mean  subject  data  which  had  been  obtained  during  several  evaluations 
of  the  Navy's  prototype  DTP  garment  (DTP1ART,  DTP2AR1,  DTP3AR1;  see 
Chapter  II.)  were  used  to  evaluate  the  accuracy  of  the  predictive  temperature 
correlation.  The  subjects  in  these  evaluations  had  mean  skin  temperatures 
which  approached  those  observed  in  the  respiratory  heat  loss  studies  (RHL1AR1) 
and  were  observed  to  be  both  shivering  and  pronouncedly  vasoconstricted.  We 
assumed  that  if  our  predicted  segmental  skin  temperatures  compared  favorably 
with  those  found  in  the  experimental  data  files,  then  we  could  presume  that 
the  correlation  would  yield  reasonably  accurate  segmental  temperatures  at  the 
more  conservative  mean  skin- temperature  range  of  31  to  32  °C. 

The  modified  Kerslake  correlation  was  derived  with  skin  temperature  data 
taken  from  five  seated,  resting  subjects  in  a  hyperbaric  helium-oxygen 
environment  at  ambient  conditions  of  20  °C  and  200  MSW  [21]  because  this  data 
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correlated  well  when  comparing  one  subject  to  another  and  when  it  was  free  of 
segmental  anomalies.  Visual  inspection  of  graphs  prepared  from  data  recorded 
during  the  three  available  immersed  studies  (DTP1ART,  DTP2AR1,  DTP3AR1) 
exhibited  aberrant  temperatures  associated  with  the  head,  hand,  and  foot  for  a 
number  of  the  experimental  subjects.  The  deviation  of  data  from  its  mean  is  a 
very  important  consideration  in  making  statistical  inferences  when  the  number 
of  experimental  subjects  is  small.  For  this  reason,  the  most  statistically 
accurate  predictive  segmental  skin  temperature  correlations  should  be  obtained 
from  the  hyperbaric  data  (RHL1AR1)  despite  significant  respiratory  heat  loss 
in  addition  to  losses  attributed  to  surface  convection  and  radiation.  We 
ignored  any  effect  that  respiratory  heat  loss  may  have  produced  on  the 
peripheral  vasomotor  state  of  these  hyperbaric  subjects  because  of  the 
observation  that  the  peripheral  thermoreceptors  play  the  major  role  in 
regulating  metabolic  heat  production  without  apparent  adjustment  for  cooling 
through  the  respiratory  tract  [21]. 

To  examine  the  accuracy  of  the  modified  Kerslake  correlation  to  represent 
the  entire  quasi -steady  period  of  thermal  equilibrium  (35  to  Q0  minutes) 
chosen  from  the  hyperbaric  records,  we  calculated  T(i)  for  each  minute  of  time 
by  means  of  Equation  10  and  the  parameters  of  Table  10.  The  independent 
variable  Tskm  was  computed  with  Equation  8  from  segmental  temperatures  which 
were  determined  from  the  experimental  linearized  equations  (derived  from  data 
in  the  RHL1AR1  file).  Finally,  each  predicted  segmental  temperature  was 
compared  to  the  equivalent  linearized  experimental  value  for  the  same  instant 
of  time  and  an  absolute  temperature  difference  was  obtained.  These 
computations  were  conducted  twice:  once  with  Rl  of  Equation  10  equal  to  a 
constant  as  defined  in  Table  10,  and  then  with  Rl  computed  at  each  instant  of 
time.  The  results  of  these  comparisons  are  summarized  in  Table  11.  Also 
tabulated  for  both  the  variable  and  constant  Rl  cases  is  the  percentage  of  the 
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Table  11 

Percentage  of  the  Total  Number  of  Compared  Points  (56)  Over  the  Span 
of  35  to  90  Minutes  Which  Have  Absolute  Temperature  Differences 
£  0.50  °C  and  _<  0.25°C.  Evaluations  for  both  R1  constant  and  R1  as 
a  function  of  time  are  presented. 


Percentages 


R1  = 

constant 

R1  » 

Rl(time) 

Segment 

Name 

<  0.25  e 

'C  <  0.50  °C 

<  0.25  ° 

C  <  0.50  °C 

Head 

100 

100 

76.8 

100 

Chest 

69.6 

100 

100 

100 

Abdomen 

100 

100 

100 

100 

Upper  Back 

80.4 

100 

100 

100 

Lower  Back 

100 

100 

100 

100 

Arm 

100 

100 

83,9 

100 

Hand 

100 

100 

100 

100 

Front  Thigh 

100 

100 

ICO 

100 

Front  Calf 

10O 

100 

100 

100 

Rear  Thigh 

69.6 

100 

62,5 

100 

Rear  Calf 

91.1 

100 

100 

100 

rOQt 

3Q.4 

63 

100 

100 
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total  number  of  compared  points  which  had  absolute  differences  of  £  0.50  °C 
and  jC  0.25  “C.  For  R1  evaluated  as  a  function  of  time,  the  predicted  skin 
temperature  for  each  segment  is  always  within  C.50  °C  of  the  corresponding 
linearized  experimental  value.  When  R1  is  constant,  all  segments  except  the 
foot  compare  within  0.50  °C;  however,  the  foot  satisfies  the  0.50  °C  criteria 
for  63i  of  the  time  considered. 

Further  authentication  of  the  modified  Kerslake  correlation  was  conducted 
with  the  actual  experimental  data  from  which  the  correlation  was  derived 
(RHL1AR1)  and  with  three  sets  of  experimental  segmental  temperatures, 
collected  as  described  by  Zumrick  [30],  from  immersed  divers  during  NCSC 
evaluations  of  the  DTP  garment  [22]  (DTP1ART,  DTP2AR1,  DTP3AR1).  DTPlART's 
and  DTP2ARl‘s  data  were  obtained  from  divers  immersed  in  water  at  a  depth 
equivalent  to  3  MSW,  and  DTP3ARl's  data  was  obtained  from  divers  submerged  to 
a  depth  equivalent  to  21  MSW.  The  number  of  experimental  subjects  in  each  of 
the  four  studies  (RHL1AR1 ,  DTP1ART,  DTP2AR1 ,  and  DTP3AR1)  were  5,  4,  3,  and  3, 
respectively. 

The  authentication  process  was  conducted  with  a  normalized  temperature 
(Eq.  11}  which  permits  comparison  of  experimental  studies  with  different  mean 
skin  (Tskm)  and  ambient  (Ta)  temperatures. 

t»(u  -  HILJjl  on 

Tskm  -  T 

a 


where: 


TN(i)  «  normalized  skin  temperature  for  segment  i 

T(i)  s  Predicted  or  experimental  skin  temperature  for  segment  i 

Ta  =  Experimental  ambient  temperature 

Tskm  »  Experimental  Hody  mean  skin  temperature 


|  Mean  segmental  temperature  profiles,  which  were  prepared  for  each  experimental 
study,  were  analyzed  for  a  period  of  75  minutes.  This  period  began  at  the 
30th  minute  of  each  experimental  record  after  the  suggestion  of 
Piantadosi  et  al .  [22]  that  transitory  effects  of  peripheral  vasoconstriction 
were,  by  then,  complete.  Each  data  file  was  analyzed  at  discrete  time 
intervals  where  a  normalized  temperature  for  each  of  Body's  twelve  segments 
i  was  calculated.  The  normalized  temperatures  were  calculated  by  Equation  11 
with  a  mean  skin  temperature  (Tskm)  value  determined  by  Equation  8  from  the 
mean  experimental  data.  With  this  same  value  of  Tskm,  segmental  predicted 
temperatures  were  calculated  for  the  same  time  inverval  by  means  of 
Equation  10  and  the  parameters  defined  in  Table  10.  The  normalized  predicted 
skin  temperature  was  then  determined  with  Equation  11.  Both  the  experimental 
and  the  predicted  normalized  temperatures  were  averaged  over' time  and  are 
tabulated  in  Tables  12-15  for  data  files  RHL1AR1,  DTP1ART,  0TP2AR1,  and 
DTP3AR1,  respectively.  The  standard  deviation,  determined  when  the  normalized 
temperatures  were  time  averaged,  is  also  listed  for  each  of  the  twelve 
segments.  To  quantify  the  accuracy  of  the  predictive  correlation  for 
estimating  the  actual  experimental  temperatures,  the  time  averaged  predicted 
TN(i  )preo  3n<*  exPerim6nta1  normalized  'TTTif i ) £ x P  femP6ratures  were  employed  to 
calculate  a  percent  difference  with  the  normalized  experimental  temperature  as 
the  computational  base  [Eq.  12).  This  percent  difference  is  also  tabulated  in 
Tables  12-15.  A  negative  percent  difference  indicates  that  the  normalized 
experimental  temperature  is  greater  than  the  predicted  value  as  shown  below. 
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Table  12 

Time  Averaged  (Mean)  Normalized  Temperatures  for  Each  of  Hody’s 
Twelve  Segments  Computed  From  the  Mean  Temperature  Profiles  of 
Experimental  File  RML1AR1  and  the  Equivalent  Predicted 
Normalized  Temperature  Which  was  Estimated  by  Equations  10  and  11. 
The  percent  difference  between  the  segmental  predicted  and 
experimental  time  averaged,  normalized  temperatures  was  computed 
by  Equation  12. 


Time  Averaged  Tskm: 

25.70  °C, 

Std. 

Dev. : 

0.2464 

Time  Averaged  Ta: 

20.03  °C, 

Std. 

Dev. : 

0.3162 

Normalized  Temperatures 
Experimental  Predicted 


Segment 

Mean 

Std.  Dev.* 

Mean 

Std.  Dev.+ 

%  Diff 

Head 

1.208 

4.24E-02 

1.188 

1.21E-02 

-1.66 

Chest 

1.146 

7.41E-02 

1.141 

9.03E-03 

-0.44 

Abdomen 

1.160 

3.63E-02 

1.157 

1.Q1E-Q2 

-0.26 

Upper  8ack 

1.268 

4.89E-02 

1.271 

1.75E-02 

0.24 

Lower  8a  ck 

0.963 

1.64E-02 

0.958 

2.59E-03 

-0.52 

Arm 

1.121 

4. 16E-02 

1.115 

7.45E-03 

-0.54 

Hand 

0.642 

7.36E-02 

0.666 

2. 15E-02 

3.74 

Front  Thigh 

0.993 

2.70E-Q2 

0.979 

1.43E-03 

-1.41 

Front  Calf 

0.924 

3.21E-02 

0.916 

5.37E-03 

-0.87 

Rear  Thigh 

0.963 

5.136*02 

0.969 

2.04E-03 

0.62 

Rear  Calf 

0.719 

4.33E-02 

0.713 

1.84E-02 

-0.83 

Foot 

0.526 

1.31E-01 

0.551 

2.S9E-02 

4.75 

*E  format  notation  is  equivalent  to  scientific  notation;  i.e.,  E-02  ■  IG~?. 


Time  Averaged  (Mean)  Normalized  Temperatures  for  Each  of  Hody's 
Twelve  Segments  Computed  From  the  Mean  Temperature  Profiles  of 
Experimental  File  DTP1ART  and  the  Equivalent  Predicted 
Normalized  Temperature  Which  was  Estimated  by  Equations  10  and  11. 
The  percent  difference  between  the  segmental  predicted  and 
experimental  time  averaged,  normalized  temperatures  was  computed 
by  Equation  12. 

Time  Averaged  Tskm:  27.77  °C,  Std.  Dev.:  0.7740 
Time  Averaged  Ta:  4.04  °C,  Std.  Dev.:  0.1316 

Normalized  Temperatures 
Experimental  Predicted 


Segment 

Mean 

Std.  0ev.T 

Mean 

Std.  Dev.' 

%  Oiff 

Head 

0.903 

6.53E-02 

1.045 

1. 61E-n3 

15.73 

Chest 

1.111 

4.03E-02 

1.034 

1.03E-03 

-  6.93 

Abdomen 

1.131 

1.84E-02 

1.037 

1.24E-03 

-  8.31 

Upper  Back 

1.155 

4.28E-02 

1.065 

2.28E-03 

-  7.79 

Lower  Back 

1.084 

1 . 0  IE- 02 

0.990 

6.18E-04 

-  8.6-7 

Arm 

1.090 

1.56E-02 

1.027 

1.03E-03 

-  5.78 

Hand 

1.116 

3.16E-02 

2.86E-03 

-17.56 

Front  Thigh 

0.712 

8.51E-Q2 

0.995 

2.76E-Q4 

39. 7S 

Front  Cal f 

0.933 

2.56E-02 

0.980 

5.52E-04 

5.04 

Rear  Thigh 

0.877 

2.28E-02 

0.993 

5.52E-04 

13.23 

Rea'-  Calf 

0.938 

* 

£- 

4 

r» 

1 

O 

fs> 

0.932 

2.49E-03 

-  0.64 

Foot 

0.907 

3.25E-02 

0.393 

3.8IE-03 

-  1.54 

f£  format  notation  is  equivalent  to  scientific  notation;  t.e.,  E-02  *  10"^. 
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Table  14 

Time  Averaged  (Mean)  Normalized  Temperatures  for  Each  of  Body's 
Twelve  Segments  Computed  From  the  Mean  Temperature  Profiles  of 
Experimental  File  DTP2AR1  and  the  Equivalent  Predicted 
Normalized  Temperature  Which  was  Estimated  by  Equations  10  and  11. 
The  percent  difference  between  the  segmental  predicted  and 
experimental  time  averaged,  normal ized  temperatures  was  computed 
by  Equation  12. 

Time  Averaged  Tskm:  28.58  °C,  Std.  Dev.:  1.2526 
Time  Averaged  Ta:  4.77  °C,  Std.  Dev.:  0.0500 

Normalized  Temperatures 
Experimental  Predicted 


Segment 

Mean 

Std.  Dev. 

Mean 

Std.  Dev.’ 

1  Diff. 

Head 

0.980 

3.87E-02 

1.045 

2.26E-03 

6.63 

Chest 

0.937 

1. 15E-02 

1.033 

1.68E-03 

10.25 

Abdomen 

0.977 

9. 368-03 

1.037 

1.77E-03 

6.14 

Upper  Back 

1.194 

3.21E-02 

1.Q65 

3.29E-03 

o 

CD 

• 

O 

fH 

t 

Lnwer  Back 

1.102 

1.62E-Q2 

0.990 

4.78E-04 

-10.16 

Arm 

1.103 

1.73E-02 

1.027 

1 . 301-03 

-  6. 89 

Hand 

0.985 

6.16E-02 

0.921 

4.02E-03 

-  6.50 

Front  Thigh 

0.863 

1. 15E-02 

0.995 

3.91E  *04 

14.63 

Front  Cal f 

1 .018 

2.04E-0? 

0.980 

8.29E-04 

-  3.73 

Rear  Thigh 

0.863 

1.42E-Q2 

0.993 

3.91E-04 

15.06 

Rear  Calf 

0.905 

2.37E-02 

0.932 

3.43C-03 

2.98 

Foot 

G.983 

7. ? IE-03 

0.893 

5.39E-Q3 

-  9.16 

f£  format  notation  is  equivalent  to  scientific  notation;  i.e.,  E-Q2  = 
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Table  15 

Time  Averaged  (Mean)  Normalized  Temperatures  for  Each  of  Body's 
Twelve  Segments  Computed  From  the  Mean  Temperature  Profiles  of 
Experimental  File  DTP3AR1  and  the  Equivalent  Predicted 
Normalized  Temperature  Wnich  was  Estimated  by  Equations  10  and  11. 
The  percent  difference  between  the  segmental  predicted  and 
experimental  time  averaged,  normalized  temperatures  was  computed 
by  Equation  12. 

Time  Averaged  Tskm:  26.64  °C,  Std.  Dev.:  0.6198 

Time  Averaged  Ta:  4.25  °C,  Std.  Dev.:  0.3600 

Normalized  Temperatures 
Experimental  Predicted 


Segment 

Mean 

Std.  Dev.^ 

Mean 

Std.  Dev.’*’ 

%  Diff. 

Head 

0.576 

5.90E-02 

1.047 

8.50E-04 

81.77 

Chest 

1.133 

1.83E-02 

1.036 

7.18E-04 

-  8.56 

Abdomen 

1.091 

5.80E-02 

1.039 

7.87E-04 

-  4.77 

Upper  Back 

1.273 

4.46E-Q2 

1.068 

1.20E-03 

-  16.10 

Lower  Back 

1.154 

2.84E-02 

0.989 

3.21E-04 

-  14.30 

Arm^  t 

Hand 

0.966 

4.08E-02 

0.916 

1.73E-03 

-  5.18 

Front  Thigh 

0.95° 

3.04E-02 

0.995 

7.876-04 

3,86 

Front  Calf 

1.155 

3.25E-02 

0.979 

5.55E-G4 

-  15.24 

Rear  Thigh 

0.845 

5.30E-02 

0.492 

3,-  2  IE- 04 

17.40 

Rear  Calf 

0.930 

5.45E-02 

0.928 

1.40E-03 

-  0.2? 

Foot 

0.283 

1.86E-02 

0.887 

2.20E-03 

213.43 

format  notation  is  equivalent  to  scientific  notation.,  i.e.,  E-02  =»  1G~^. 
^Insufficient  data  recorded  from  the  arm  to  permit  comparative  analysis. 
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of 

to 


MffO) .  ( --^P!1D ' ™(1)Exr ) 

™(1 )  EXP 


100 


where:  %  Oiff(i)  =  Percent  difference  for  segment  i 

7N(i ) pred  =  Time  avera9®d  predicted  temperature  for  segment  1 

I 

TN(i ) exp  3  Time  averaged  experimental  temperature  for  segment  i  j 
Table  12  displays  the  percent  differences  for  all  regions  except  the  j 

j 

head,  hand,  front  thigh,  and  foot  to  be  less  than  1%  for  the  RHL1AR1  data.  ; 

This  result  was  expected  since  the  modified  Kerslake  correlations  were  j 

1 

developed  with  this  data.  Table  13,  which  is  based  upon  experimental  file 

DTP1ART,  shows  percent  differences  within  10%  for  all  regions  except  the  head  i 

i 

(18.73%) ,  hand  (17.56%),  front  thigh  (39.75%),  and  rear  thigh  (13.23%).  The 

! 

larger  error  associated  with  the  front  thigh  is  probably  an  anomaly  since  ! 
considerably  better  agreement  is  found  for  the  front  thigh  in  the  two  other 
submerged  studies  that  were  analyzed  (Tables  14  and  15).  Table  14,  which  is 
derived  from  experimental  file  DTP2AR1,  shows  percent  differences  which  are 
less  than  or  equal  to  15%  for  all  regions. 

During  the  21  MSW  submerged  evaluation  of  the  DTP  garment  (Table  15,  from 
data  file  DTP3AR1),  insufficient  data  was  recorded  from  the  arm  to  permit 
comparative  analyses.  The  percent  differences  associated  with  the  head 
(81.77%)  and  the  foot  (213.43%)  were  obviously  anomalous  when  compared  to  the 
significantly  better  values  for  the  head  and  foot  in  Tables  13  and  14.  A 
review  of  the  0TP3AR1  data  revealed  that  erroneous  data  was  collected  from  the 
head  and  foot  over  much  of  the  75  minute  period  of  interest. ^  Only  those  data 


^The  75  minute  period  of  time,  which  begins  at  the  30th  and  ends  at  the 
105tn  minute  of  each  experimental  file,  was  chosen  to  establish  a 
representative  common  period  of  exposure  among  the  four  experimental  f!les 
considered. 
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I 

points  which  could  be  assumed  to  be  representative  of  the  experimental 

phenomenon  were  considered  in  the  analyses  reflected  in  Table  15.  1 

I 

I 

Unfortunately,  the  small  number  of  usable  data  points  heavily  biases  the  ; 

I 

results  obtained.  Of  all  the  sets  of  experimental  data,  that  recorded  during  : 

I 

the  21  meter  evaluation  of  the  DTP  garment  was  the  least  congruous.  ! 

i 

In  considering  the  results  of  our  analyses,  it  should  be  noted  that  all 

i 

of  the  experimental  data  within  the  selected  time  period  (30  to  105  minutes)  ' 

i 

i 

was  employed  in  computing  the  time  averaged  normalized  temperatures.  J 

Examination  of  the  mean  experimental  data  also  revealed  the  presence  of  i 

j 

questionable  data  points,  particularly  for  the  head,  hand,  and  foot.  i 

i 

Piantadosi  et  al .  [22]  suggests  that  some  of  the  observed  anomalies, 
artificially  high  heat  flux  values,  were  due  to  an  influx  of  cold  water  around  : 
poorly  sealed  areas  of  the  dry  suit's  outer  garment.  Leaks  occurred  between 

j 

the  suit  sleeve  and  glove,  which  compromised  the  suit's  insulation  in  this 
area,  and  air  often  collected  in  the  garment's  hood,  which  precipitated 
leakage  around  the  hood-face  seal  [22j. 

We  concluded  from  these  results  that  the  modified  kerslake  equation 
(Eq.  10  and  Table  10)  can  be  employed  to  predict  segmental  skin  temperatures 
over  a  mean  skin  temperature  range  of  26  to  28  "C.  Extrapolation  of  the 
correlation  to  Tskm  values  above  this  envelope  i$  probably  reasonable,  but  is 
dependent  upon  levels  of  vasoconstriction  and  shivering  -  two  variables  over 
which  we  had  no  control  during  data  acquisition  or  the  numerical  analyses. 
Ideally,  actual  experimental  data  would  confirm  the  validity  of  tne 
extrapolation  to  other  Tskm  values.  Furthermore,  it  is  apparent  fros?  the 
tabulated  comparisons  of  time  averaged,  normalized  predicted  and  experimental 
skin  temperatures  (Tables  13-15}  that  the  modified  Kerslake  correlation 
approximates  the  principal  segments'  temperatures  within  s  155  error. 
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j  The  regional  predicted  skin  temperatures  of  our  model  (T^  of  Eq.  5)  are 
derived  from  these  predicted  segmental  skin  temperatures  (T(i)  of  Eq.  10  with 
Table  10).^  Table  16  displays  each  model  man  region  with  the  corresponding 

i 

t 

;  Hody  segments.  For  the  regions  that  h'-ve  two  corresponding  segments  (See 

I 

Table  16.),  the  predicted  segmental  skin  temperatures  (Eq.  10  with  Table  10) 
j  are  weight  averaged  with  the  corresponding  Hody  segmental  surface  area 
fractions  (Table  1)  to  obtain  the  regional  predicted  skin  temperature. 

:  Equation  13  shows  an  example  (the  torso)  of  the  weight  averaging  equation 

I 

|  which  is  used  to  find  the  regional  predicted  skin  temperature. 

■  T  .  W(2)T(2)  +  W(4)T(4) 

2  W(2)  +  W( 4) 

where:  J2  =  Temperature  of  region  2,  the  torso  (Fig.  1) 

W(i)  =  Hody  surface  area  fraction  for  segment  i  (Table  1) 

T(i)  »  Segmental  predicted  skin  temperatures  for  segment  i; 
computed  from  Equation  10  with  Table  10 

i=2  indicates  segment  2,  chest 

i=4  indicates  segment  4,  upper  back 

The  regional  predicted  skin  temperauture  is  equal  to  the  segmental  predicted 
skin  temperature  for  regions  that  exhibit  a  one  to  one  correspondence  with  a 
segment  (for  example,  the  head;  see  Table  16.) 


*°$ee  Foreword  for  connotation  of  segment  and  region. 
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Table  16 

Model  Man  Region  Names,  Indices,  and  Corresponding  Hody  Segments 


Regional 

Index 

Region 

Name 

Region 

Abbreviation 

Corresponding  Hody 

Segment  Indices  (Table  1) 
(segment  name) 

1 

Head 

HEAD 

1  (head) 

2 

Torso 

TORSO 

2  (chest),  4  (upper  back) 

3 

Abdomen 

ABDOMEN 

3  (abdomen),  5  (lower  back) 

4 

Th  i  gh 

THIGH 

8  (front  thigh),  10  (rear  thigh) 

5 

Calf 

CALF 

9  (front  calf),  11  (rear  calf) 

6 

Foot 

FOOT 

12  (foot) 

7 

Upper  Arm 

UP  ARM 

6  (arm) 

8 

Lower  Arm 

LOW  ARM 

6  (arm) 

9 

Hand 

HAND 

7  (hand) 

CHAPTER  V 


DEVELOPMENT  OF  OVERALL  HEAT  TRANSFER  COEFFICIENT 


As  indicated  in  Chapter  III,  the  overall  heat  transfer  coefficient,  U-j  ( j ) 
(Eq.  6),  has  been  developed  by  using  the  thermal  characteristic  of  the  DTP 
garment  as  bc’ng  representative  of  a  pressure  compensated  dry  suit.  The 
analysis  presented  below  is  also  valid  for  dry  suit  configurations  which 
behave  similarly  to  the  DTP  garment,  but  appropriate  garment  characteristics 
would  have  to  be  substituted  for  those  associated  with  the  DTP  garment. 
Equation  6  is  presented  below  to  facilitate  the  discussion. 

(6) 

1 

Vj)  =  - 

'  t  Sj(j)  •  In  (r2/r1 )  +  S, (j )  •  In  (r3/r2)  +  S, (j) 

"i  0)  2%  •  “  2%  •  «  S„(j)h0(j) 

As  stated  previously,  the  first  term  of  Equation  6  is  neglected  since  we 
assume  that  the  suit  is  always  subjected  to  a  hydrostatic  pressure 
differential,  below  the  shoulder  level,  which  causes  the  garment  to  rest  on 
the  skin  surface.  Thus,  hypothetically,  no  gas  layer  can  exist  between  the 
skin  and  undergarment.  Since  the  OTP  suit's  outer  garment  is  constructed  of  a 
0.16  CM  thick  crushed  neoprene  foam,  the  outer  garment  can  be  modeled  as  solid 

t!32) 
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neoprene  with  the  same  nominal  thickness  of  0.16  CM  (1/16  inch).  Thus,  the 
outer  garment  can  be  treated  as  being  incompressible  and  impermeable  to  the  ' 

j 

breathing  gases  and  water.  Substituting  the  skin  surface  area  of  a  nodal 
element,  S-j(j)  (Eq.  14),  into  the  third  term  of  Equation  6  (Rqq)  and  canceling  1 

|  Si(j)  =  2nr1Az  (14)  | 

where:  S-j(j)  =  Nodal  skin  surface  area  (M2) 

i  I 

rx  =  Radius  from  regional  center  line  to  skin  surface  (M) 

!  Az  =  Nodal  height  ; 

i  | 

!  the  2IIAZ  term  in  both  the  numerator  and  denominator  yields  Equation  15. 


06 


_  rj  •  In  (r3/r2) 


og 


The  natural  logarithm  in  Equation  15  corrects  tor  the  surface  area  variation 
of  a  cylinder. 


Resistance  of  the  Outer  Boundary  Layer 

Substituting  the  nodal  skin  surface  area,  S-j(j)  (Eq.  14),  and  the  outer 
surface  area  of  the  crushed  neoprene  garment,  S0(j)  (Eq,  16),  into  the  last 

S0(j)  =  2ltr3AZ  (16) 

where:  S0(j)  «  Nodal  skin  surface  area  of  outer  surface  of  the  outer 
garment  (M2) 

r3  «  Radius  from  regional  center  line  to  the  outside  surface 
to  the  outer  garment  (M) 

AZ 


=  Nodal  height 


term  of  Equation  6  (Robl)  an<*  canceling  the  2n&z  term  in  both  the  numerator  and 
denominator  yields  Equation  17. 


R0BL  = 


r3ho 


(17) 


Since  the  empirical  film  coefficients  (presented  below)  were  derived  from 
experimental  data  which  was  recorded  at  one  location  per  region,  these  values 
are  assumed  to  be  regional  averages.  Thus,  the  nodal  indicator,  ' j  * ,  of  h0(j) 
in  Equation  6  has  been  dropped. 


Convective  heat  transfer  coefficient  (h0) 

The  experimental  respiratory  heat  loss  data  (RHL1AR1)  was  utilized  to 
obtain  empirical  correlations  for  the  mixed  mode  convective  film  coefficient 
(h0).  The  correlations  are  expressed  in  terms  of  a  dimensionless  Nusselt 
number  as  a  function  of  the  Grashof  number  and  Prandtl  number  product. 

Figures  65-67  are  plots  of  the  logarithm  (base  10)  cf  the  Nusselt  number 
versus  the  logarithm  (base  10)  of  the  Grashof  number  and  Prandtl  number 
product.  Figure  65  displays  the  data  derived  from  the  RHL1AR1  file  for  the 
torso  (TORSO),  abdomen  (ABDOMEN),  upper  arm  (UP  ARM),  lower  arm  (LOW  ARM), 
calf  (CALF),  and  foot  (FOOT)  (Fig.  1.).  Since  these  regions  may  be  taken  as 
vertical  cylinders  for  a  sitting  subject,  Figure  65  also  displays  the  curve 
generated  from  the  classical  empirical  free  convection  correlation  for 
vertical  cylinders  (theory)  "Nu  =  0.59  •  (GrPr)°*25  [39].  Figures  66  and  67 
present  the  data  derived  from  the  RH11AR1  file  for  the  thigh  (THIGH),  and  the 
head  (HEAD)  and  hand  (HAND),  respectively.  For  a  sitting  subject,  the  thigh  is 
treated  as  horizontal  cylinders,  and  the  head  and  hand  are  considered  to  be 
spheres.  Figures  66  and  67  also  display  the  curve  produced  from  the  classical 
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Figures  65-67.  Composite  Experimental,  Classical  Empirical  Correlation  for 
Natural  Convection  (THEORY)  and  Our  Empirically  Derived 
Correlation  (FITTED)  for  Vertical  Cylinders,  Horizontal 
Cylinders  and  Spheres.  Log  of  Nusselt  number  plotted  against 
the  log  of  the  Grashof  number  and  Prandtl  number  product. 


Number  Versus  log  of  Grashof  Number  and  Prandtl  Number  Product 


Figure  66.  log  of  Nusselt  Number  Versus  Log  of  Grashof  Number  and  Prandtl  Number  Product 
for  Horizontal  Cylinder. 


Number  Versus  log  of  Grashof  Number  and  Prandtl  Number  Product 
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empirical  correlation  for  free  convection  (theory)  for  horizontal  cylinders, 

Nu  =  0.53  •  (GrPr)0*25  [39],  and  spheres,  Nu  »  2.0  +  0.43  •  (GrPr)0*25  [40,41], 
respectively. 

Since  the  RHL1AR1  data  was  recorded  from  sitting,  resting  subjects 
exposed  to  a  20  °C  hyperbaric  95$  helium  -  5%  oxygen  environment  at  a  pressure 
equivalent  to  200  MSW,  the  property’s  necessary  to  determine  these 
dimensionless  parameters  were  determined  at  a  constant  depth  of  200  MSW  and 
and  at  an  average  film  temperature  of  22.15  °C.  The  film  temperature  was 
taken  as  the  arithmetic  mean  between  the  regional  skin  and  the  ambient 
temperatures;  this  average  film  temperature  was  the  mean  of  the  highest  and 
lowest  film  temperatures  during  the  120  minute  duration  of  the  experimental 
studies.  The  specific  heat,  Cp,  and  the  density,  p,  were  interpolated  from 
data  tabulated  in  the  U.S.  Navy  Diving  Gas  Manual  [42].  The  viscosity,  u,  and 
thermal  conductivity,  k,  were  determined  for  a  binary  mixture  by  using  the 
Wilke  correlation  [43]  in  the  viscosity  equation  and  the  Lindsay  and 
Bromley  [43]  formulation  in  the  thermal  conductivity  equation.  The 
constituent  properties  of  helium  and  oxygen  were  interpolated  as  a  function  of 
temperature  and  pressure  from  values  in  the  U.S.  Navy  Diving  Gas  Manual  [42]. 
Appendix  C  contains  example  calculations  for  the  viscosity,  u,  thermal 
conductivity,  k,  and  the  other  properties  mentioned  above. 

The  experimental  data  in  file  RHllARl  is  stored  in  the  twelve  segment 
format  of  Hody  (Fig.  Z),  but  the  regions  mentioned  above  refer  to  the  nine 
regions  of  the  model  man  (Fig.  i).J*  To  convert  the  data  in  RH'.lARi  to  the 
nine  region  model  forswt.  the  experimental  segments  corresponding  to  a 
mode!  man’s  region  were  weight  averaged  with  the  Hody  surface  area  fractions 
which  correspond  to  the  equivalent  model  region.  (See  Eq.  13  for  an  example.) 

’^See  Foreword  for  the  connotation  of  segment  and  region. 


140 


The  differences  between  the  Nusselt  correlations  derived  from 
experimental  data  and  the  classical  empirical  curves,  displayed  in 
Figures  65-67,  are  attributed  to  a  forced  convection  component  superimposed  on 
the  natural  convection.  We  could  not  evaluate  the  accuracy  of  accepted  mixed 
mode  convective  heat  transfer  correlations  since  no  information  concerning 
forced  draft  gas  currents  was  recorded  during  the  experiments.  Since  no 
velocity  information  of  any  kind  was  recorded  during  these  experiments, 
empirical  correlations  (Table  17)  were  derived  from  the  composite  experimental 
data  of  Figures  65-67.  The  equation  for  vertical  cylinders  (Table  17)  was 
derived  by  fitting  a  parabola  to  the  data  of  the  torso,  abdomen,  upper  arm, 
lower  arm,  and  foot.  The  data  from  the  calf  was  dropped  because  it  was 
considered  aberrant  when  compared,  in  dimensionless  form,  to  the  other 
vertical  cylindrical  regions  (Fig,  65).  As  a  result,  we  assumed  that  the  calf 
obeyed  our  derived  correlation  for  a  vertical  inder.  This  anomaly  is 
possibly  a  result  of  the  transition  from  a  ve  '  to  a  horizontal  cylinder 
that  occurs  at  the  knee  for  a  sitting  subject.  Figure  66  displays  the  single 
set  of  data  available  for  a  horizontal  cylinder.  Since  this  data  is  too  far 
removed  from  the  classical  empirical  curve  to  simply  use  the  classical 
empirical  equation,  we  decided  that  it  would  be  reasonable  to  4's>ume  that  the 
experimental  data  follows  a  line  parallel  to  the  classical  empi-ical  line. 
Therefore,  the  formula  for  the  horizontal  cylinders  *as  defined  zs  a  line  with 
a  slope  matching  that  of  the  classical  empirical  line,  but  with  a  different 
y-intercept.  The  equation  for  spheres  pre  ’nted  in  Table  17  is  the  result  of 
fitting  a  straight  'me  between  the  tw1  sets  of  data  shown  in  Figure  67. 
Figures  65-67  also  display  the  corresponding  empirical  correlation  (FITTED) 
which  is  defined  in  Table  17. 

Because  of  the  poor  agreement  between  the  experiuaenta!  data  and  the 
classical  e^pirica’  free  convection  equations,  we  consider  that  our  espi^ical 
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Table  17 

Empirical  Convective  Heat  Transfer  Correlations 


Region  Geometrical 
_  Model 

Torso  Vertical  log  "Nil 

Cylinder 

Abdomen 
Upper  An'>! 
j  Lower  Arm 
Foot 
Calf 

Thigh  Horizontal  log  "NiJ 

Cylinder 

i 

i 

Head  Sphere  log  "Nu 

I 

!  Hand 


Correlation 

a2  (log  (GrPr))2  +  al  log  (GrPr)  +  a0 
a2  =  0- 1817 
ai  =  -2.4128 
a0  ■  9.1934 


a:  log  (GrPr)  +  a0 
aL  =  0.2486 

3g  3  -0.4483 

a!  log  (GrPr)  +  a0 
aj  =  0.46486 
a0  3  -1.9726 
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|  correlations  (Table  17)  are  a  more  accurate  representation  of  the  experimental 

i  events  than  the  classical  empirical  equations  for  free  convection.  Therefore, 

] 

I  our  empirical  correlations  (Table  17)  have  been  employed  to  estimate  the 
resistance  due  to  the  outer  boundary  layer  (Rqbl  Eq.  17)  1n  the  overall 
heat  transfer  coefficient. 


Undergarment  Thermal  Resistance 

For  the  case  of  a  diver  who  is  wearing  a  pressure  compensated  dry  suit, 
the  predominant  thermal  resistance  is  provided  by  an  undergarment  which  is 
fabricated  of  either  a  fibrous  material  or  an  insulative  foam.  The 
experimental  data  which  was  used  in  these  developments  was  collected  from 
resting  subjects  who  were  wearing  the  DTP  suit  with  an  undergarment  fabricated 
of  fibrous  material  (M-400  Thinsulate  made  by  3-M  Corporation).  We  assumed 
that  if  the  model  was  correct  for  the  fibrous  material,  it  also  could  be 
utilized  for  other  suit  ensembles  which  satisfy  the  same  basic  assumptions. 

The  thermal  resistance  of  this  fibrous  undergarment  varies  because  of 
permeation  by  the  ambient  environmental  gas  inside  the  dry  suit  and 
compression  (suit  squeeze)  which  results  from  a  local  hydrostatic  pressure 
differential  between  the  interior  of  the  garment  and  the  surrounding  water. 

Although  the  variation  in  the  properties  of  Thinsulate  have  not  been 
evaluated  under  hyperbaric  conditions,  tests  have  been  performed  at  1  ATA  with 
air  and  other  gases  [13].  These  tests  revealed  that  Thinsulate's  insulative 
quality,  as  well  as  that  of  similar  porous  materials,  is  due  primarily  to  its 
ability  to  entrap  gas  without  the  establishment  of  micro-convective  currents 
within  the  material.  For  this  reason,  we  assumed  that  the  principal  mode  of 
heat  loss  through  the  undergarment  is  conduction.  We  further  assumed  that  the 
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experimental  thermal  resistance  data  obtained  in  air  can  be  extrapolated  to 
any  hyperbaric  environment  by  multiplying  the  experimental  air  data  by  the 
ratio  of  the  thermal  conductivities  of  the  air  and  the  ambient  hyperbaric 
environment  (Eq.  18).  The  assumed  method  of  extrapolation  correlates  well 
with  experimental  findings  in  helium  at  a  pressure  of  1  ATA  [13].  This  method 
of  computation  yields  the  baseline  specific  thermal  resistance  (Rc)^  of  the 
undergarment  at  any  depth  of  immersion  (evaluated  at  the  shoulder  crest)  and 
for  any  compensating  suit  gas. 

R't„(0)  .  R*a(  )  (18) 

where:  R'j^(O)  =  Specific  thermal  resistance  of  Thii  te  which  has 
been  corrected  for  effects  of  the  entrapped  gas  at 
ambient  conditions 

R‘a  »  Specific  thermal  resistance  of  Thinsulate  at  sea  level 
in  air  at  21.1  °C  (2.27  CLO/CM)  as  determined  by  Audet, 
Orner,  and  Kupferman  [12] 

ka  =  Thermal  conductivity  of  air  at  sea  level  and  21.1  °C 

kg  =  Thermal  conductivity  of  entrapped  gas  at  ambient 

conditions 

The  immersed  upright  diver  is  subjected  to  a  differential  pressure 
because  of  the  linearly  increasing  hydrostatic  head.  This  differential 
pressure  is  zero  at  the  neckline  and  approximately  13.8  kilopascals  (2  PS t G) 
at  the  feet  of  a  standing  man.  Audet  et  [12]  developed  a  method  of 
predicting  the  change  in  specific  thermal  resistance,  R'  (CIO/CM  of 
thickness),  as  a  function  of  the  chance  in  material  density  because  of  this 
hydrostatic  pressure  di fferential .  They  assumed  that  the  density  change  in 

^Specific  thermal  resistance  (R1)  is  defined  in  terms  of  CIO/CM  of 
thickness.  1  CIO  =  0.18  m2.HR°C/KCAL  or  in  SI  units,  0.155  M?°C/W. 


the  material  was  directly  proportional  to  the  change  in  thickness. 

Furthermore,  they  assumed  that  the  degradation  in  R'  as  a  function  of  the 
change  in  density  (thickness)  in  fibrous  batt  materials  (Thinsulate)  was 
analogous  to  that  which  was  observed  in  foam  garment  insulating  materials. 

The  behavior  that  they  observed  in  these  foam  insulating  materials  (high 
density  foams)  indicated  that,  as  the  density  increased  (thickness  decreased), 
the  specific  thermal  resistance  (R*)  decreased  linearly.  Adopting  this  model 
for  the  change  in  the  specific  thermal  resistance  of  the  Thinsulate 
undergarment.  Figure  68  was  derived  from  experimental  data  obtained  by 
Audet  et_  aU_  [12].  At  a  compressive  pressure  of  2  PSIG,  the  Thinsulate  had  a 
specific  thermal  resistance  of  2.11  CLO/CM  and  the  thickness  had  been  reduced 
by  60%.  The  decimal  percent  change  in  R'  (DPCSR)  which  corresponds  to  the 
2.11  CLO/CM  value  is  calculated  by  Equation  19  to  be  0.07048,  and  the  60% 
change  in  thickness  corresponds  to  a  decimal  percent  change  in  thickness 
(DPCT)  of  0.6C  (Eq.  20). 


DPCSR 


(19) 


where:  DPCSR  3  Decimal  percent  change  in  specific  thermal  resistance 

R'0  =  Uncompressed  specific  thermal  resistance.  This  was  taken 

from  Audet  et  ah  to  be  2.27  CIO/CM  at  an  uncompressed 
thickness  of  l7fT3  CM. 


R' 


a 


Specific  therma1  resistance  at  elevated  pressure;  at 
2  PSIG  it  is  equal  to  2.11  CLO/CM 


Figure  68.  Linear  Relationship  Describing  the  Undergarn*ent 's  Specific  Themal 
Resistance  Change  as  a  Function  of  Change  in  Thickness. 
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DPCT  =  to  I  t  (20) 

t 

o 

where:  OPCT  -  Decimal  percent  change  in  thickness 
t0  =  Uncompressed  thickness  (1.63  CM) 

t  =  Thickness  after  compression 

In  determining  the  gradient  plotted  in  Figure  68,  it  was  assumed  that,  in  the 
uncompressed  state,  the  Thinsulate  has  no  degradation  in  specific  thermal 
resistance  or  in  thickness.  The  slope  (DR)  of  the  line  (Fig.  68)  which  passes 
from  (0.0, 0.0)  to  (0.60,0.07048)  was  determined  to  be  0.1175.  This  line 
illustrates  the  hypothesis  that,  as  the  undergarment  thickness  is  decreased, 
the  specific  thermal  resistance  is  also  decreased  linearly.  Therefore,  the 
specific  thermal  resistance  corrected  for  hydrostatic  suit  squeeze  may  be 
predicted  at  any  depth  below  the  shoulders  between  0.0  CM  and  136.8  CM 
(2.0  PS  I G)  by  the  derived  equation: 

R*TH(Z)  -  R'th(O)  (1  -  [DR  •  OPCT ( 2 ) ] )  (21) 

where:  R'y^Z)  *  Specific  resistance  of  Thinsulate  at  Z  CM  below  the 
shoulders 

R'th(O)  a  Specific  resistance  of  Thinsulate  at  the  shoulder 
crest 

DR  =  0.1175,  which  is  the  slope  of  the  line  (Fig,  68)  that 
describes  the  decimal  percent  change  in  specific 
thermal  resistance  (DPCSR)  versus  the  decimal  percent 
change  in  thickness  (DPCT) 

0PCT(Z)  =  Decimal  percent  change  In  thickness  of  the  undergarment 
at  Z  CM  below  the  shoulders 

The  variation  in  the  decimal  percent  change  in  thickness  (DPCT)  for  the 
fibrous  batt  material  (Thinsulate),  as  a  result  of  the  hydrostatic  head  which 


is  measured  from  the  shoulders,  is  predicted  by  either  of  the  following 
empirical  equations: 

OPCT  =  2.358  x  5.0"2  •  Z  ,  (22) 

For  Z  _<  17.11  CM  (0.25  PSIG). 

DPCT  =  0.343  +  3.784  x  10-3  •  Z  -  1.345  x  10"5  •  Z2  ,  (23) 

For  17.11  CM  (0.25  PSIG)  <  Z  _<  136.84  CM  (2.0  PSIG). 

Z  is  the  depth  in  CM  below  the  shoulder  crest.  Equations  22  and  23  were 
fitted  to  experimental  data  (Fig.  69)  which  was  taken  on  flat  samples  by 
Audet  et  al .  [12].  The  dashed  line  in  Figure  69  indicates  that  a  break  point 
was  assumed  at  0.25  PSIG  as  a  result  of  our  analysis  of  the  experimental  data 
of  Audet  et  al . ,  which  is  reproduced  in  Figure  70.  This  figure  displays  the 
percent  change  in  thermal  resistance  (CL0)  as  a  function  of  suit  compressive 
pressure  (PSIG).  Inspection  of  Figure  70  reveals  a  rather  sharp  gradient 
between  0.25  PSIG  and  the  origin,  but  shows  that  the  "ate  of  percent  change  in 
thermal  resistance  as  a  function  of  the  imposed  pressure  moderates  at  presures 
above  0.25  PSIG.  This  observation  suggests  that  there  should  be  a  similar 
break  point  in  the  curves  plotted  in  Figure  69  for  the  percent  change  in 
thickness  as  a  function  of  the  compressive  load.  Therefore,  we  assumed  a 
linear  rate  of  change  in  the  Thinsulate's  thickness  between  0.0  PSIG  and 
0.25  PSIG  with  a  more  moderate  curvilinear  rate  from  0.25  PSIG  to  0.44  PSIG. 

We  determined  the  break  point  of  Figure  69  as  that  point  of  intersection 
between  an  extension  of  the  line  between  the  experimental  data  points  at 
0.44  PSIG  and  0,50  PSIG  and  the  constant  compressive  pressure  intercept  at 
0.25  PSIG.  Th»  thermal  resistance  values  (R  in  CIO)  obtained  by  using 
cquations  22  or  23,  21,  and  24  for  the  pressure  range  of  0.0  to  0.44  PSIG 
compare  to  within  10%  of  those  displayed  in  Figure  70. 
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PRESSURE  (PSI) 

Figure  69.  Percent  Change  in  Thickness  as  a  Function  of  Compressive  Pressure.  Re-drawn  from  a  graph 
presented  by  Audet,  Orner,  and  Kupferman  [12]. 
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rcent  Change  in  Thermal  Resistance  as  a  Function  of  Compressive  Pressu 
graph  presented  by  Audet,  Orner ,  and  Kupferman  [12]. 
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rTH  “  r‘th(Z)  (ti  -  (ti  *  DPCT)) 
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(24) 


where:  R-fn  s  Thermal  resistance  of  Thinsulate  undergarment  ( CLO) 

R * ( Z )  -  Specific  thermal  resistance  of  Thinsulate 
undergarment  (CLO/CM) 

t^  =  Uncompressed  thickness  of  undergarment  (CM) 

OPCT  =  Decimal  percent  change  in  thickness  from 

Equation  22  or  23 


To  correct  the  specific  thermal  resistance  of  a  nodal  element  for  the 
effects  of  the  suit  entrapped  gas  and  hydrostatic  suit  squeeze,  the  following 
procedure  is  used.  First,  the  undergarment's  uncompressed  R'  value  of 
2.27  CLO/CM  [12]  at  sea  level  in  21.1  °C  air  is  corrected  for  the  effects  of 
the  suit  entrapped  gas  at  the  ambient  conditions  (Eq.  18).  This  corrected  R' 
value  is  then  used  as  the  shoulder  crest  value  in  Equation  21.  To  evaluate 
the  decimal  percent  change  in  thickness  (DPCT)  at  the  node's  depth  below  the 
shoulders,  the  dimensions  of  the  body  regions  (Table  5)  and  the  pasture  of  the 
subject  are  considered.  Depending  upon  the  magnitude  of  this  aepth  (Z), 
either  Equation  22  or  23  ts  utilized.  The  resulting  value  for  the  decimal 
percent  change  in  thickness  (DPCT)  is  then  substituted  into  Equation  21.  This 
value  represents  the  specific  thermal  resistance  a'ter  it  has  been  corrected 
for  the  ?ffect$  of  the  ambient  entrapped  gas  and  the  reduction  in  thickness  of 
the  undergarment  due  to  hydrostatic  suit  squeeze.  The  unut f garment *s  thermal 
resistance  is  no*  computed  by  converting  R'y^(Z)  (CLO/CN)  to  a  thermal 
resistance  (CLO)  through  multiplication  by  the  undergarment's  estimated 
thickness  at  a  depth  of  Z  CM  below  the  shoulders  (Eq,  24).  The  value 

in  CIO's  is  then  converted  to  51  units  (Rtk)  (Se e  footnote  12  on  page  143.) 
and  s-^stituted  into  Equation  25. 
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Rjh  *  1n  (r2/rl) 

Ar 


(25) 


where:  Ryg  =  Thermal  resistance  of  the  undergarment  corrected  for  suit 
entrapped  gas  and  hydrostatic  pressure  (m2°c/W) 


RTH  =  Thermal  resistance  of  undergarment  at  ambient  conditions 


^r  =  r2  -  rL  =  t-j  •  DPCT 

tj  _  1.63  CM  which  is  the  uncompressed  undergarment 
thickness  [12] 

DPCT  =  Decimal  percent  change  in  undergarment  thickness  from 
Equations  22  or  23 


The  natural  logarithm  in  Equation  21  corrects  the  Rt^  *alue  for  the 
effects  of  the  region's  curved  surface. 

Overall  heat  transfer  coefficient 

Neglecting  any  insulation  which  might  have  been  contributed  by  a  gas 
layer  between  the  skin  and  undergarment  because  of  suit  squeere,  Equation  6  is 
rewritten  below  by  substituting  Ryg  (Eq.  25),  Pqq  (Eq.  15),  and  Rygy  (Eq.  1? J 
for  terms  2  through  4,  respectively. 

1 

U^j)  «  _ - 

ri  *  RTH  *  (r2^rl }  ♦ 

Ar 


In  (r3/r2) 


09 


r , 

r  h 
3  o 
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CHAPTER  VI 


THE  INTERACTIVE  MODEL  PROGRAM  • 


The  regional  temperature  prediction  correlations  (Chapter  IV)  and  the 
overall  heat  transfer  coefficient  formulation  (Chapte1-  V)  have  been 
synthesized  into  a  numerical  model  for  predicting  the  regional  heat  loss  and 
required  regional  supplementary  heating  of  a  resting,  immersed  diver.  The 
body  geometry  is  that  described  in  Chapter  III.  The  numerical  model  assumes 
that  the  predicted  regional  temperatures  (Eg.  10  with  Table  10  and  Eq.  13)  are 
uniform  over  the  entire  regional  surface  and  that  t'"«**e  is  no  thermal 
interaction  between  any  of  the  regions.  This  model  program  assumes  that  the 
diver  will  be  in  a  prone,  sitting,  or  standing  posture  while  wearing  the  ft. 5. 
Navy's  prototype  Diver  Th-;nnal  Protection  garment  that  was  described  in 
Chapter  III.  Further,  the  initial  model  is  restricted  to  a  breathing  gas  of 
air.  This  Imitation  was  necessary  because  all  of  the  available  experimental 
data  from  which  the  model  was  developed  was  obtained  with  air  as  the  breathing 
gas;  however,  this  does  not  restrict  the  model  to  only  air  diving.  The  model, 
as  programmed,  should  yield  reasonable  results  with  other  breathing  gases 
under  the  prescription  that  the  thermal  conductivity  of  these  gases  as  a 
function  of  temperature  an<j  pressure  can  be  predicted.  The  prograBi  computes 
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the  regional  supplementary  heating  requirement  by  subtracting  the  regional 
allowable  heat  flux  from  the  corresponding  regional  predicted  heat  fit*  value. 
The  regional  allowable  heat  fluxes  which  are  used  by  the  model  program  are 
those  attributed  to  Burriss  [33,34]  and  are  tabulated  in  Table  18. 

The  model  program  allows  the  user  to  perform  the  heat  loss  computations 
for  the  entire  region  as  a  whole  or  for  many  finite  difference  nodes  when  the 
region  is  vertical.  When  a  region  approximates  a  horizontal  cylinder  or  a 
sphere,  the  program  assumes  that  the  resistance  of  these  regions  may  be 
accurately  estimated  by  using  thermophysical  properties  which  are  determined 
at  the  lateral  midline  (Fig  71). 


HORIZONTAL  SPHERE 

CYLINDER 


Figure  71.  lateral  Mid-lines  for  Horizontal  Cylinder  and  Sphere. 


154 


Table  18 

Regional  Allowable  Heat  Fluxes 
Excerpted  from  Burriss  [33,34] 

Region  Region  Allowable  Heat 

ID  #  Name  Flux  (W/M?) _ 


1 

Head 

23.30 

2 

Torso 

60.00 

3 

Abdomen 

49.70 

4 

Thigh 

42.40 

5 

Calf 

85.00 

6 

Foot 

97.00 

1 

l 

Upper  Arm 

98.90 

8 

Lower  Arm 

125.25 

9 

Hand 

?AC,86 
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i  The  thermal  resistance  of  vertical  regions  when  considered  to  be  single  nodes 

i 

j  are  computed  with  properties  that  are  determined  at  the  midpoint  of  their 
i  height  (Fig.  71).  When  the  vertical  regions  are  divided  into  discrete  nodal 

j 

;  elements,  the  thermal  resistance  is  determined  for  each  node  separately  by 
using  properties  that  are  determined  at  the  respective  nodal  depth  (Fig.  72). 


t 


Mid-point  of 
the  Height, 
L/2 


A2 


Mid-point  of 
aNode,  AZ/2 


Figure  72.  Single  Nodal  Element  for  a  Vertical  Body  Region  and  Multi-nodal 
Element  for  a  Vertical  Body  Region. 
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The  model  program  consists  of  ten  modules  which  are  named  Main,  REINIT, 

OHTC,  F(F),  CHTC,  TR,  PRED,  TC,  PLOTD,  and  PLOTN.  Documentation  describing 

I 

the  use  of  the  interactive  model  program  is. contained  in  Appendix  D;  each 
module  is  flow  charted  in  Appendix  E.  Appendix  F  contains  a  complete  source 
code  listing  of  each  module.  Each  module  is  listed  below  with  a  brief 
description  of  its  function  (See  Appendices  E  and  F  for  details.),  and  an 
explanation  of  the  source  (derivation)  of  any  properties  which  are  introduced 
in  the  mociel.  Appendix  Q  is  a  sample  calculation  which  illustrates  the 
computations  that  are  performed  in  executing  a  simple  production  mode  run  (See 
Appendix  D  for  explanation.).  A  summary  of  all  the  data  that  is  used  by  the 
program  is  included  at  the  front  of  Appendix  G. 


Main  Program 

The  Main  module  acts  as  a  controller  which  governs  the  sequence  of 
execution.  It  collects  all  of  the  information  that  is  necessary  to  define  the 
problem  and  then  proceeds  to  call  the  other  modules,  as  required,  to  compute 
the  nodal  (one  or  many)  heat  loss(es).  The  regional  rate  of  heat  loss  is 
determined  by  summing  the  nodal  rates  of  heat  loss  and  the  regional 
supplementary  heating  requirement  is  then  computed  by  subtracting  a  regional 
rate  of  allowable  heat  loss  from  this  predicted  regional  rate  of  heat  loss. 
This  module  synthesizes  the  information  provided  to  determine  which  mode  of 
execution  to  follow.  If  the  user  indicates  the  production  mode,  then  a  design 
value  of  the  mean  skin  temperature  and  of  the  ambient  temperature  must  be 
entered.  If  a  research  run  is  indicated,  the  program  allows  the  u:ti  to 
compare  predicted  and  experimental  values  of  temperature,  heat  flux,  overall 
heat  transfer  coefficient,  normalized  temperature,  and  Biot  number.  The 
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questions  asked  and  the  consequences  of  the  possible  responses  are  described 
|  in  Appendix  0,  and  the  mechanics  of  the  program's  execution  is  described  in 

i 

;  the  flow  chart  (Appendix  E)  and  the  source  code  listing  (Appendix  F).  All  of 

| 

the  properties  that  are  defined  in  the  Main  module  have  been  previously 
;  described  in  Chapters  III,  IV,  V,  and  VI. 


RE  I  NIT 

This  function  subroutine  is  provided  to  allow  the  operator  to  return  to 
the  program-assigned  values.  The  Main  module  allows  the  DTP  garment 
properties  (Chapter  III)  and  thicknesses  (Chapter  III),  the  regional  allowable 
heat  fluxes  (Chapter  III  and  VI),  and  the  overall  heat  transfer  coefficient 
biasing  factors  (which  will  be  described  ir.  Chapter  VII)  to  be  dynamically 
altered.  This  routine  returns  these  values  to  their  program-assigned  values. 
See  questions  5  thru  8  of  the  Question  Flow  section  of  Appendix  D  for  the 
program-assigned  values  and  the  method  of  altering  these  values. 


;  OHTC 

This  su’.r ,-utine  determines  the  overall  heat  transfer  coefficient  for  a 
specified  nooe  with  the  method  and  equation-  ',hat  are  described  in  Chapter  V. 

i  The  values  for  Rqg  (^*  15)  and  Ryg  (Eq.  25)  are  determined  in  this 
module.  The  value  for  Rqbl  17)  is  found  by  using  the  Secant 
Method  [44,45]  for  estimating  the  temperature  of  the  outside  surface  of  the 
outer  garment  (T3).  The  Secant  algorithm  utilizes  the  module  F(T)  to 

i 

|  ascertain  whether  the  estimated  value  of  T3  is  correct.  Until  the  energy 


balance  of  Function  F(T)  is  approximately  equal  to  zero,  the  Secant  algorithm 
continues  to  iterate  for  T3. 

The  value  for  the  convective  heat  transfer  coefficient  (h0  of  Eq.  17)  is 
computed  at  temperature  T3  with  module  CHTC.  All  properties  that  are  used  in 
this  module  are  defined  in  the  Main  module.  Because  of  the  limits  imposed  on 
air  diving  bu  the  U.S.  Navy  [46],  we  have  restricted  the  model  program  to 
depths  _<  250  FSW  (76.2  MSW,  111.37  PSIG),  and  the  ambient  temperature  must  be 
<  1.0  °C  to  remain  within  the  limits  defined  for  the  boundary  layer  fluid 
properties. 

F(T) 

This  is  a  function  subroutine  which  calculates  the  heat  flux  which  passes 
through  the  complete  garment  ensemble  and  comDares  it  with  the  heat  flux  that 
only  passes  through  the  undergarment  and  outer  garment  (F(T)  in  Eq.  27). 

These  computations  are  carried  out  by  assuming  that  the  temperature  of  the 
outside  surface  of  the  outer  garment  (T)  is  T3  as  passed  from  subroutine 


RUG  +  ROG  + 


r3  •  CHTC(T,r3) 


where: 


T  =  Temperature  T3,  passed  from  the  Secant  Method 

computation  of  Function  OHTC 

T1  =  Temperature  of  skin  surface  (°C) 

Ta  =  Ambient  temperature  (°C) 

RUG  =  Thermal  resistance  of  undergarment 

ROG  =  Thermal  resistance  of  outer  garment 

rL  =  Radius  of  body  region 

r3  =  Radius  to  outside  of  outergarment 

CHTC(T,r3)  =  Function  subroutine  which  estimates  the  convective 
heat  transfer  coefficient  at  temperature  T  and 
radius  r3 


All  of  the  properties  needed  by  Equation  27  sre  passed  from  module  OHTC.  If 
the  estimated  temperature  T  is  correct,  F(T)  will  equal  zero. 


CHTC 

The  module  CHTC  determines  the  convective  heat  transfer  coefficient  for 
the  outer  (water)  boundary  layer  (h0)  by  using  our  derived  empirical  equations 
(Table  17)  which  are  described  in  Chapter  V.  The  subroutine  first  determines 
whether  the  region  is  a  vertical  or  horizontal  cylinder,  or  a  sphere.  For 
both  the  horizontal  cylinders  and  spheres,  the  characteristic  length  (CL)  is 
defined  to  be  twice  the  radius  of  the  outer  garment  (R3).  The  vertical 
cylinders  use  the  height  as  the  characteristic  length  as  they  are  defined  in 
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i 

I 

Table  5  of  Chapter  III.  A  value  for  the  Grashof  number  and  Prandtl  number  j 

j 

product  (GrPr)  is  now  determined  with  Equation  28.  ' 

GrPr  =  a  •  CL3  •  (T3  -  Ta)  (28)  j 


where: 


13  = 

^  = 
CL  ~ 


Temperature  of  outside  surface  of  outer  garment  (°C) 
Ambient  temperature  (°C) 

9ep.£p  (  _  )  ,  described  below 


yk 


M3°C 


The  property  grouping,  a  (Eq.  ~8),  is  determined  with  Equation  29. 


where:  a 


p 

CP 

u 

k 

Tf 

T3 


a  *  log-1  [9.80714  (log  (Tf ) )0,12733] 
_  g8p2CD 


(29) 


yk 

=  Acceleration  of  gravity,  9.81  M/SEC^ 

=  Volumetric  coefficient  of  expansion  (°K"1)  for  water 
=  Density  of  water  (K&/M-1) 

=  Specific  heat  at  constant  pressure  (J/KG°K)  for  water 

=  Dynamic  viscosity  (KG/M. SEC)  for  water 

=  Thermal  conductivity  of  water  (U/M°K) 

=  Film  temperature  (°K)  =  T3  £ 

2.0 

=  Outside  surface  temperature  of  outer  garment  (aK) 
a  Ambient  temperature  (CK) 


log*1  x  »  10* 


The  constants  9,80714  and  0.12733  of  Equation  29  are  curve  fitting 
coefficients  that  were  determined  when  a  power  curve  was  fitted  to  the  data  of 
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I  Table  19  which  was  excerpted  from  Holman  [40].  Log*3  T  and  log  a  of  Table  19 


correspond  to  x  and  y,  respectively,  in  Equation  30. 

y  =  axb 


Table  19 

a  as  a  Function  of  Temperature 
Excerpted  from  Holman  [40] 


im 

log  T 

a  (1/M3oC) 

log  a 

4.444 

0.6478 

1.91  x  10+9 

9.281 

10.000 

1.0000 

6.34  x  10+9 

9.80209 

15.555 

1.1919 

1.08  x  10+1° 

10.033 

21.111 

1.3245 

1.46  x  10+1° 

10.164 

This  curve  fit  was  performed  with  the  HP-65  Stat  PAC  1  routine  Stat  1-24A  [47] 
for  power  curve  fitting.  The  determined  value  or  a  (Eq.  29)  is  substituted 
into  Equation  28,  which  yields  an  estimate  for  the  product  of  the  Grashof  and 
Prandtl  numbers.  This  value  of  GrPr  is  substituted  into  the  appropriate 
equation  of  Table  17  to  produce  an  estimate  for  the  Nusselt  number  which  is 
associated  with  the  garment  water  interface. 


nu  ,  iaEL 


where:  Nu  -  Nusselt  number  (dimensionless) 

h0  =  Convective  heat  transfer  coefficient 
Cl  o  Characteristic  length  (H) 

k  =  Thermal  conductivity  of  water  (U/N°;<)  at  the  film 
temperature  (Tf) 


l3log  is  to  base  10. 


Solving  Equation  31  for  the  convective  heat  transfer  coefficient  (h0)  requires 
a  value  for  the  thermal  conductivity  of  water.  This  is  determined  by 
Equation  32  which  is  a  curve  that  is  fitted  to  tabulated  thermal  conductivity 
data  that  was  excerpted  from  Holman  [40]  (Table  20), 

k  =  0.56662  +  T  *  1.7977  x  10'3  (32) 

where:  k  =  Thermal  conductivity  of  water  (W/M°C) 

T  =  Film  temperature  (°C) 

This  fresh  water  thermal  conductivity  is  used  because  all  of  the  available 
immersed  experimental  data  was  recorded  from  divers  who  were  submerged  in 
fresh  water  and  because  the  thermal  conductivity  of  resting  or  slowly  moving 
sea  water  may  be  approximated  by  that  of  fresh  water  [48], 

Table  20 

Thermal  Conductivity  of  Water  as  a  Function  of  Temperature 
Excerpted  from  Holman  [40] 


T  (°C) 

k  (W/M°C) 

0.000 

0.566 

4.444 

0.575 

10.000 

0.585 

15.555 

0.595 

21.111 

0.604 

Equation  32  was  fitted  to  the  data  of  Table  20  with  the  HP-65  Stat  PAC  1 
routine  Stat  1 - 22A  [47]  for  linear  regression.  Substituting  the  thermal 
conductivity  value  that  is  estimated  by  Equation  32  into  Equation  31  allows 
the  convective  heat  transfer  coefficient  (h0)  to  be  determined. 
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TR 

This  subroutine  predicts  an  estimate  of  the  model  man's  regional  skin 
temperature  by  weight  averaging  the  corresponding  predicted  segmental 
temperatures  (Eq.  13  and  Eq.  10  with  Table  10  of  Chapter  IV).  Table  16  of 
Chapter  IV  displays  the  correspondence  between  the  model  man's  regions  and  the 
experimental  segments  which  were  used  in  developing  the  predicted  temperature 
correlations  of  Equations  13  and  10  with  Table  10. 


PREO 


This  subroutine  predicts  a  segmental  skin  temperature  by  using 
Equation  10  with  the  constants  of  Table  10  (See  Chapter  IV.}.  These  segmental 
skin  temperatures  are  the  values  which  are  weight  averaged  in  Function  TR, 
above. 


TC 

Module  TC  determines  the  thermal  conductivity  of  air  for  use  in 
Equation  IS  of  Chapter  V.  The  value  is  interpolated  from  tabulated  air  data, 
which  was  excerpted  from  the  U.S.  Navy  Diving  Gas  Manual  [423 *  as  3  function 
of  temperature  and  pressure.  The  two-dimensional  interpolation  routine  was 
taken  from  Carnahan,  Luther,  and  Wilkes  [<*9],  Because  of  the  limits  imposed 
on  air  diving  by  the  U.S.  Navy  [46],  we  have  restricted  the  model  program  to 
depths  <  250  FSW  (76.2  MSW,  111.37  P5IG),  and  the  •‘.nbient  temperature  must  be 


164 

<  1.0  °C  to  remain  within  the  limits  defined  for  the  boundary  layer  fluid 
properties. 


PLOTD 

This  subroutine  takes  dimensional  temperatures,  heat  fluxes,  and  overall 
heat  transfer  coefficients  (both  experimental  and  predicted)  and  graphs  them 
as  a  function  of  time.  This  data  is  generated  whenever  predicted  and 
experimental  data  are  compared.  These  graphs  are  produced  with  the  Digital 
Equipment  Corporation's  (DEC)  PLXY-11M  plotting  routines  [50]  and  are 
printed  on  a  DEC  LXY11  Printer/Plotter  (Printronix  P3Q0). 


PLGTN 

The  PIOTN  module  uses  the  experimental  and  predicted  temperatures,  heat 
fluxes,  and  overall  heat  transfer  coefficients  to  compute  dimensionless, 
normalized  temperatures,  TN  (Eq.  33),  and  Biot  numbers  [40],  81  (tq.  34). 

TN  *  .  Ti  .~-.V  (33) 

Tskm  -  T 

3 

where;  TN  =  Regional  normalized  experimental  or  predicted  skin 
temperature 

Tj  =■  Regional  expei^mental  or  predicted  skin  temperature 
Ta  =  Experimental  ambient  temperature 

Tskm  =  Hody  mean  skin  temperature  (Eq.  8);  determined  with  the 
experimental  data 
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BJ  =  Hill  (34) 

k 

9 

where:  81  =  Experimental  or  predicted  Biot  number  [40] 

Ui  =  Experimental  or  predicted  overall  heat  transfer  coefficient 

rj  =  Radius  from  the  region's  center  line  to  the  skin  surface 

kg  =  Thermal  conductivity  of  the  entrapped  breathing  gas  (air); 
determined  at  the  experimental  or  predicted  film  temperature 

These  experimental  and  Predicted  normalized  temperatures  and  Biot  numbers  are 
then  graphed  as  a  function  of  the  time.  Each  plot  (normalized  temperature  or 
Biot  number)  shows  both  the  experimental  and  the  predicted  values.  These 
plots  are  produced  in  a  manner  simila-  to  that  described  in  the  PLOTO 
section. 


CHAPTER  VII 


MODEL  VERIFICATION,  CONCLUSIONS,  AND  RECOMMENDATION 


Empirical  Correlation  of  the  Model 

Inspection  of  f.quation  5  (Chapter  IIIl  indicates  that  the  rate  of  heat 
loss  from  each  of  the  hypothetical  nodal  elements  (Qi ( j ) )  is  computed  as  the 
product  of  the  respective  nodal  surface  area  ( S-j  (j } ) ,  overall  heat  transfer 
coefficient  per  unit  of  skin  area  (Ui(j)),  and  the  temperature  differential 
between  the  regional  skin  surface  (T^)  and  the  ambient  environment  (Ta).  The 
data  analysis  performed  in  Chapter  IV  revealed  that  Equation  10  would  predict 
the  segmental  ?kin  temperatures  (T(i))  with  an  accu"’.y  of  15%  for  the  mean 
skin  temperature  range  of  26  to  29  °C  and  that  the  regional  values  (T^)  are 
predicted  by  weight  averaging  the  segmental  values.  Table  6  displays  the  good 
agreement  which  exists  between  the  surface  area  percentages  of  our  model  man 
(Fig.  ’.)  and  the  accepted  experimental  surface  area  percentages  of  Hody,  thus 
indicating  the  accuracy  of  the  model  men's  dimensions  (Table  5).  To  complete 
the  model,  we  investigated  the  accuracy  of  the  empirical  Uj(j)  terra  of 
Equation  26  by  comparing  its  predictions  for  each  region  (Fig.  1}  with  a 
representative  experimental  overall  heat  transfer  coefficient  for  the  same 
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region.  This  comparison  is  exoressed  in  the  form  of  a  dimensionless  Biot 
number  (BI),  Equation  34,  which  Is  repeated  below. 

B!  =  Hill  (34) 

k 

9 

Expressing  the  overall  heat  transfer  coefficient  in  this  dimensionless  form 
neutralizes  the  influence  of  the  ambient  temperature  and  pressure  on  the 
ensemble's  thermal  conductance  and  thus,  permits  the  comparison  of 
experimental  data  collected  at  various  ambient  depths  and  temperatures.  As 
may  be  seen  in  Equations  26,  18,  21,  and  24,  the  effects  of  the  ambient 
temperature  and  pressure  are  restricted  to  the  first  term  of  the  denominator 
of  our  empirical  expression  for  U-j(j)  (Eq.  26).  Inspection  of  Equations  18 
and  21  indicates  that  any  variation  in  the  undergarment's  thermal  resistance 
due  to  the  ambient  pressure  and  temperature  is  primarily  caused  by  changes  in 
the  thermal  conductivity  of  the  gas  with  which  the  suit  is  pressure 
compensated.  The  other  terns  of  the  denominator  of  Equation  26  ar*  considered 
essentially  invariant  over  the  ambient  temperature  and  pressure  range  in  which 
the  immersed  experiments  (0TP1ART,  DTP2AR1,  DTP3A31;  ser  Chapter  II.)  were 
conducted. 

To  determine  regional  uiasing  factors  for  correcting  the  hypothetical 
overaH  conductance  of  Equation  26,  we  employed  experimental  data  that  was 
collected  during  three  evaluations  of  the  DTP  gansent  (DTPiAR?,  i/TPjAftl , 
QTP3AR1).  The  i r strumentation  system  which  was  utilized  for  data  collection 
was  analogous  to  that  described  by  Zgsnrtck  [30]  and  was  calibrated  with  the 
technique  which  was  developed  by  Nucfco’s  [19].  Two  of  the  evaluation  studies 
(OTP  1 ART  and  0TP2ASI)  were  conducted  with  divers  isrsersed  in  water  at  a 
pressure  equivalent  to  3  and  the  third  (DTP3ARI)  was  conducted  at  a 
pressure  equivalent  of  21  MS'*.  The  n-^ber  of  esperi-wntal  subjects  **  each  of 
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I 

t 

the  three  studies  (DTP1ART,  DTP2AR1,  and  DTP3AR1)  was  4,  3,  and  3,  j 

j 

respectively.  For  each  experimental  study,  representative  segmental  i 

i 

I 

temperature  and  rate  of  heat  loss  profiles  (stored  in  DTP1ART,  QTP2AR1,  and  ; 

I 

DTP3AR1.)  were  prepared  by  numerically  averaging  the  data  from  the  diver 

1 

! 

subjects  at  each  time  interval.  This  process  is  described  in  detail  in 


Chapter  II.  The  ambient  temperature  of  the  water  ranged  between  2  and  5  °C. 


Except  for  the  head  and  feet,  the  garment  ensemble  worn  by  these  prone,  1 

i 

I 

resting  divers  consisted  of  a  fibrous  batt  (Thinsulate)  undergarment  of  j 

t 

1.63  CM  nominal  thickness  at  1  ATA  and  an  outer  garment  of  crushed  foam 


neoprene  which  has  a  nominal  thickness  of  0.16  CM.  The  head  was  covered  by  a 


conventional  neoprene  hood  and  the  feet  were  insulated  by  a  double  thickness 
of  the  fibrous  batt  material.  The  passive  protection  garment  was  pressure 
compensated  with  air  at  a  pressure  equivalent  to  the  ambient  hydrostatic 
pressure  at  the  shoulder  crest.  All  three  of  the  test  series  are  documented 


in  Chapter  II. 


To  obtain  experimental  values  that  are  equivalent  to  the  nine  region  1 

model  (Fig.  1)  from  the  segmental  data  of  the  composite  data  records  (DTP1ART,  j 

i 

DTP2AR1,  DTP3AR1),  experimental  values  of  temperature  and  rate  of  heat  loss  | 


from  the  corresponding  Hody  segments  were  averaged  by  using  the  corresponding 

I 

I 

Hody  weighting  factors  (Table  1)  (See  Table  16.).  From  these  profiles  and  a  i 

knowledge  of  the  ambient  temperature,  experimental  overall  conductance  (U^(j))  i 

! 

profiles  were  derived  as  a  function  of  time.  These  conductance  profiles  were  ' 

j 

then  converted  to  Biot  number  profiles  with  Equation  34.  To  discount  any 
transitory  effects  of  peripheral  vasoconstriction,  each  of  the  experimental 
Biot  number  profiles  was  examined  for  a  period  of  75  minutes  beginning  at  the 
30th  minute,  after  the  suggestion  of  Piantadosi  et_  a]^  [22].  Similarly, 
predicted  conductances,  as  a  function  of  time,  were  calculated  for  the  same 
time  period  by  means  of  Equation  26  and  plotted  as  predicted  Biot  numbers.  i 


j 
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Both  the  experimental  and  the  predicted  regional  Biot  number  profiles  were 
averaged  over  time  and  are  tabulated  in  Tables  21,  22,  and  23  for  DTP1ART, 
DTP2AR1,  and  DTP3AR1  data  files,  respectively.  The  standard  deviations  that 
were  determined  when  the  profiles  were  time  averaged  are  also  listed  for  each  ! 

i 

region.  The  head  data  is  not  displayed  in  these  tables  because  of  the  j 

anomalies  found  in  the  experimental  profiles  for  this  region.  To  determine 

the  accuracy  of  the  predictive  correlation  (Eq.  26),  the  time  averaged 

_  1 

predicted  and  experimental  Biot  numbers  (BI)  were  used  in  Equation  35  to  } 

j 

calculate  a  percent  difference.  These  regional  percent  difference  values  for  j 
the  three  DTP  studies  are  also  tabulated  in  Tables  21,  22,  and  23.  Negative 
values  indicate  that  the  experimental  Biot  number  is  greater  than  the  I 

predicted.  j 


BIipRED  -  BIiEXp 

%  DIFFi  =  (  - — -  )  •  100 

BIi  EXP 


(35) 


i 

where:  %  DIFF-j  =  Percent  difference  for  region  i  ! 

J 

Blipncn  =  Time  averaged  predicted  Biot  number  for  region  i  j 

"BTjEXP  =  Time  averaged  experimental  Biot  number  for  region  i 


.  Because  of  the  large  percent  differences  (%  DIFF)  that  are  displayed  in 
Tables  21  through  23,  biasing  factors  were  obtained  for  each  region  by 
averaging  the  time  averaged  experimental  and  predicted  Biot  numbers  of 
Tables  21,  22,  and  23;  calculating  a  decimal  percent  difference  (DP)  in  a 
manner  analogous  to  Equation  35;  and  determining  the  biasing  factor  by  means 
of  Equation  36.  The  study-averaged  values  of  the  experimental  and  predicted 
Biot  numbers  for  all  model  regions,  except  the  head,  and  their  respective 
biasing  factors  (BF)  are  tabulated  in  Table  24. 


1  +  DP 
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Table  21 

Time  Averaged  (Mean)  Biot  Numbers  for  Each  Model  Region  (Fig.  1) 
Computed  From  the  Composite  Profiles  of  Experimental  File  DTP1ART 
and  the  Equivalent  Predicted  Biot  Number.  Percent  difference 
between  the  respective  regional  experimental  and  predicted  tirr.o 
averaged  Biot  numbers. 

Biot  Numbers 

Experimental  Predicted 


ID 

Region 

Mean 

Std.  Dev.'*’ 

Mean 

Std.  Dev."*" 

%  Diff 

1 

Head 

2 

Torso 

3.65E+01 

4.20E+00 

1.01E+01 

1.00E-03 

-72.3 

3 

Abdomen 

2. 30E+01 

1.43E+00 

1.01E+01 

2.93E-03 

-56.1 

4 

Thigh 

1.93E+01 

2.08E+00 

5.32E+00 

1.16E-03 

-72.4 

5 

Calf 

1.42E+01 

8.15E-01 

4.36E+00 

8.15E-04 

-69.3 

6 

Foot 

2.16E+01 

7.80E-01 

2.96E+00 

2.13E-04 

-86.3 

7 

Upper  Arm 

7.32F+00 

8.18E-01 

3.59E+00 

4. 56E-04 

-51.0 

8 

Lower  Arm 

5.34E+00 

5.96E-01 

2.75E+00 

4.90E-04 

-48.5 

9 

Hand 

1.17E+01 

1.00E+00 

4.44E+00 

1.54E-03 

-62.1 

+E  format  notation  is  equivalent  to  scientific  notation;  i.e., 


E-02  =  10'2 
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Table  22 

Time  Averaged  (Mean)  Biot  Numbers  for  Each  Model  Region  (Fig.  1) 
Computed  From  the  Composite  Profiles  of  Experimental  File  DTP2AR1 
and  the  Equivalent  Predicted  Biot  Number.  Percent  difference 
between  the  respective  regional  experimental  and  predicted  time 
averaged  Biot  numbers. 

Biot  Numbers 


Experimental  Predicted 


ID 

Region 

Mean 

Std.  Dev.+ 

Mean 

Std.  Dev.1" 

%  Dlff. 

1 

Head 

2 

Torso 

2.07E+01 

3.82E+00 

1.01E+01 

2.39E-03 

-51.2 

3 

Abdomen 

2.27E+01 

1.38E+00 

1.01E+01 

2.51E-03 

-55.5 

4 

Thigh 

1.58E+01 

1.61E+00 

5.33E+00 

1.02E-03 

-66.3 

5 

Calf 

1.43E+01 

1.38E+00 

4.37E+00 

9.D5E-04 

-69.4 

I 

6 

Foot 

1.45E+01 

1.32E+00 

2.96E+00 

3.73E-04 

-79.6  i 

7 

Upper  Arm 

9.14E+00 

6.62E-01 

3.60E+00 

8.29E-04 

-60.6  ! 

1 

8 

Lower  Arm 

6.67E+0Q 

4.82E-01 

2.76E+00 

8.27E-04 

-58.6 

9 

Hand 

1.67E+01 

2.95E+00 

4.45E+00 

2.53E-03 

-73.4 

*E  format  notation  is  equivalent  to  scientific  notation;  i.e. ,  E-02  =  10~2. 


Time  Averaged  (Mean)  Blot  Numbers  for  Each  Model  Region  (Fig.  1) 
Computed  From  the  Composite  Profiles  of  Experimental  File  DTP3AR1 
and  the  Equivalent  Predicted  Biot  Number.  Percent  difference 
between  the  respective  regional  experimental  and  predicted  time 
averaged  Blot  numbers. 

Biot  Numbers 


Experimental  Predicted 


JD 

Region 

Mean 

Std.  Dev.1,  Mean 

Std.  Dev.* 

%  Diff 

1 

Head 

2 

Torso 

2.25E+01 

3.88E+Q0  1.01E+01 

8.95E-03 

-55.1 

3 

Abdomen 

2.30E+01 

3.45E+00  1.01E+01 

9.50E-03 

-56.1 

4 

Thigh 

2.06E+01 

1.48E+00  5.32E+00 

4.27E-03 

-74.2 

5 

Calf 

1.35E+01 

2.02E+00  4.36E+00 

3.31E-03 

-67.7 

6 

Foot 

7.67E+01 

2.26E+01  2.96E+00 

1.56E-03 

-96.1 

7 

Upper  Armt* 

3.59E+00 

2.69E-03 

8 

Lower  Arm** 

2.75E+00 

1.93E-03 

9 

Hand 

1.82E+01 

1.65E+01  4.44E+00 

6.43E-03 

-75.6 

format  notation  Is  equivalent  to  scientific  notation;  i.e.,  E-02  *  10"2. 

^Insufficient  data  recorded  from  the  arm  to  permit  comparative  analysis. 
Mean  experimental  values  of  Table  24  were  obtained  by  averaging  the  time 
averaged  Biot  numbers  of  records  DTP1ART  and  DTP2AR1,  contained  In 
Tables  21  and  22,  respectively. 


Table  24 

Mean  Experimental  and  Predicted  Biot  Numbers  Averaged  From 
the  Respective  Biot  Numbers  of  Tables  21,  22,  and  23  From 
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Files  DTP1ART,  DTP2AR1,  and  DTP3AR1,  Respectively. 
Recommended  biasing  factors  computed  by  Equation  36. 


Biot  Numbers 


ID 

Regi on 

Mean  Experimental  * 

Mean  Predicted* 

Biasing 

Factor 

1 

Head 

Same  as  ID  2 

2 

Torso 

2.66E+01 

1.01E+01 

2,63 

3 

Abdomen 

2.29E+01 

1.01E+01 

2.27 

4 

Thigh 

1.86E+01 

5.32E+00 

3.50 

5 

Calf 

1.40E+01 

4.36E+00 

3.21 

6 

Foot 

3.76E+01 

2.96E+00 

12.70 

7 

Upper  Arm** 

8.23E+00 

3.59E+00 

2.29 

8 

Lower  Arm'*''*' 

6.01E+00 

2.75E+00 

2.19 

9 

Hand 

1.55E+01 

4.44E+00 

3.49 

*E  format  notation  Is  equivalent  to  scientific  notation;  i.e.»  E-02  »  10"^. 

**Mean  experimental  values  obtained  by  averaging  the  time  averaged  Blot 
numbers  of  records  0TP1ART  and  DTP2AR1,  contained  in  Tables  21  and  22, 
respectively. 
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The  hypothetical  U-j ( j )  value  obtained  from  Equation  26  should  be 
multiplied  by  its  respective  regional  biasing  factor  prior  to  being  used  in 
estimating  the  aetual  heat  loss  from  a  resting  diver  (Eq.  5).  For  the  case  of 
the  head,  we  suggest  that  the  biasing  factor  associated  with  the  torso  be 
applied  since  the  vasomotor  reflex  of  the  torso  and  abdomen  better  approximate 
that  of  the  head  and  neck  than  does  that  of  the  limbs.  Because  it  was  highly 
Incongruous,  none  of  the  experimental  head  data  collected  during  the  three 
Immersed  DTP  evaluations  was  considered  in  our  analysis.  Subjects  often 
reported  that  air  collected  In  the  hood  during  the  experiments,  thus 
compromising  the  Inherent  thermal  insulation  and  precipitating  leakage  around 
the  hood-face  seal.  The  large  biasing  factor  associated  with  the  foot  in 
Table  24  is  a  result  of  the  large  differences  in  the  experimental  values  among 
the  three  OTP  studies.  The  wide  differences  In  the  experimental  data  is 
considered  to  be  anomalous. 


Discussion  and  Recommendation 

The  three  DTP  garment  tests  which  were  used  for  determining  the  overall 
heat  transfer  coefficient  biasing  factors  of  Table  24  were  conducted  as 
pragmatic  evaluations  of  the  garment  ensemble.  These  tests  were  never 
Intended  to  serve  the  dual  role  of  providing  physiological  heat  loss  and  skin 
temperature  data  for  thermal  modeling.  The  •  'raiments  In  total  Involved  only 
ten  subjects,  hardly  a  reliable  statistical  base  from  which  to  draw  definitive 
conclusions.  Furthermore,  several  problems  cited  In  reference  22  appear  to 
have  plagued  the  subjects  throughout  all  of  the  suit  tests.  Inadequate 
sealing  of  the  dry  suit's  outer  garment,  In  particular  a  leaking  joint  between 
the  suit  and  glove,  was  a  common  complaint  of  the  experimental  subjects. 


175 


Another  source  of  concern  was  the  hood  worn  with  the  breathing  mask.  It 
collected  air  and,  as  a  result,  compromised  the  integrity  of  the  hood-face 
seal  by  permitting  leakage  into  the  suit.  Our  discussions  with  the  principal 
investigators  of  these  tests  indicated  that  the  dry  suit  may  have  been  under 
Inflated  for  much  of  the  experimental  immersions  because  of  problems 
associated  with  the  breathing  mask  configuration.  The  principal  investigators 
thought  that  this  condition,  commonly  called  suit  squeeze,  may  have  reduced 
the  overall  insulation  of  the  garment  ensemble  by  hydrostatically  compressing 
the  insulative  undergarment.  Their  conjecture  seems  to  be  supported  by  the 
results  of  more  recent  evaluations  of  an  improved  version  of  the  same  diving 
suit.  Naval  Coastal  Systems  Center  investigators  found  that  the  total  body 
insulation  provided  by  the  suit  was  substantially  increased  when  the  problems 
noted  above  were  corrected  [51]. 

Recognizing  the  limitations  and  constraints  associated  with  the  cited 
data,  we  suggest  that  the  biasing  factors  of  Table  24  be  employed  with  the 
model  presented  until  better  correlating  factors  can  be  obtained  from  data 
collected  during  controlled  experiments.  Although  the  biasing  factors  were 
developed  by  using  the  results  of  studies  wherein  the  subjects  remained 
primarily  prone,  we  recommmend  their  application  for  a  resting  diver, 
regardless  of  his  posture.  The  variation  of  the  local  overall  conductance  due 
to  posture  Is  considered  when  the  hypothetical  value  of  U^{j)  Is  computed  by 
Equation  26.  Our  predictive  model  may  be  easily  transferred  to  other  garment 
configurations  if  the  characteristics  of  the  suit  materials  are  known  and  if 
they  satisfy  the  same  assumptions  as  developed  for  the  OTP  ensemble. 

Additional  experiments  need  to  be  conducted  not  only  for  the  prone, 
resting  diver,  but  also  for  the  sitting  and  standing  diver  while  at  rest. 
Special  care  should  be  taken  during  these  experiments  to  obtain  reproducible 
conditions  among  the  experiments,  particularly  with  respect  to  posture  and 
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ambient  environmental  conditions.  Such  experiments  will  result  in 
statistically  reliable  data  and  better  correlations  for  the  overall  heat 
transfer  coefficient  from  divers  in  conventional  attitudes.  As  well,  basic 
experiments  should  be  performed  to  Investigate  the  velocity  field  which  exists 
around  the  diver  for  the  various  postures  cited  to  permit  a  more  definitive 
development  of  the  outer  garment  film  coefficient.  Although  this 
Investigation  deals  only  with  a  model  for  a  resting  diver,  we  recommend  the 
conduct  of  similar  experiments  at  various  work  rates  in  order  to  develop  well- 
authenticated  predictive  models  for  working  divers.  We  decided  (In 
consultation  with  our  Navy  Medical  Research  and  Development  Command  sponsor) 
to  defer  the  development  of  the  working  model  until  more  experimental  data  was 
j  available.  This  decision  was  reached  because  of  the  lack  of  available  working 
j  experimental  data  recorded  in  similar  ambient  environments  and  at  the  same 
rate  of  work. 

Since  it  can  be  argued  that  some  body  regions  are  more  like  truncated 
cones  than  cylinders,  we  have  examined  the  effects  of  using  a  truncated  cone 
rather  than  a  cylinder  as  a  body  region.  Based  on  Information  that  is 
presented  by  Webb  [23]  and  on  personal  observation,  it  seems  reasonable  to 
assume  that  the  thigh  is  the  most  probable  body  region  to  be  defined  as  a 
truncated  cone.  Table  25  displays  median  data  which  was  taken  from 
reference  23  for  the  circumference  of  the  upper  and  lower  thigh.  The  radii 

I 

that  correspond  to  each  of  these  circumferences  is  also  presented. 


Table  25 


Geometrical  Data  for  the  Thigh 
Adapted  from  Data  Presented  by  Webb  [23] 


Circumference 

Radius  (R)  in  CM 

Inches 

CM 

R  *  Cl rcumference/2n 

Upper  Thigh 

25.1 

63.75 

10.15 

Lower  Thigh 

17.0 

43.18 

6.87 

Length  of  Thigh  =38.3  CM,  from  Table  5 

To  compare  the  improvement  in  the  model  that  might  be  attained  by  using  a 
cone  instead  of  a  cylinder  for  the  thigh,  the  effects  of  these  geometries  on 
the  predictive  heat  loss  model  must  be  examined.  Inspection  of  Equation  26 
indicates  that,  except  for  the  resistance  of  the  Thinsulate  undergarment 
( Rjh ) •  the  thermal  conductivity  of  outer  garment  (k0g),  and  the  convective 
heat  transfer  coefficient  (h0);  all  of  the  terms  in  the  denominator  are  solely 
functions  of  the  various  radii.  The  thermal  conductivity  of  the  outer 
garment,  k0g,  is  assumed  to  be  constant  and  is  only  a  function  of  the 
thickness.  The  thermal  resistance  of  the  undergarment,  was  empirically 
derived  from  experiments  which  were  performed  with  flat  samples,  and  the 
logarithm  of  the  ratio  of  radii  (Eq.  26)  is  intended  to  correct  for  the 
curvature  associated  with  its  application  to  a  cylinder.  The  convective  heat 
transfer  coefficients,  h0,  were  derived  from  geometry  independent  Information 
except  for  the  characteristic  length,  which  appears  in  both  the  Grashof  and 
Nusselt  numbers.  For  vertical  regions,  this  length  could  be  arbitrarily 
defined  as  either  the  vertical  height,  L,  or  the  slant  height,  S,  of 
Figure  73.  For  horizontal  regions,  this  length  is  a  characteristic  diameter. 
From  these  findings,  it  is  apparent  that  the  effects  of  the  change  in  geometry 
are  restricted  to  changes  in  the  radii  and  the  associated  slant  height. 


Figure  73.  Representative  Region  Presented  as  a  Truncated  Cone. 


To  examine  the  possibility  of  approximating  a  cone  with  a  cylinder,  we 
assuaed  that  the  surface  area  and  the  vertical  height  must  remain  the  same  as 
cited  by  Webb  [23].  This  assumption  results  from  the  hypothesis  that, 
irregardless  of  the  shapes  assuaed  for  the  various  regions,  the  model  men  aost 
continue  to  approximate  an  average  fiavy  diver.  Therefore,  the  height  and 
surface  area  of  the  regions  must  remain  unchanged.  Equation  37  is  written  ^ 
equating  the  surface  area  of  a  cylinder  and  a  truncated  cone-  For  the 
purposes  of  this  comparison,  end  effects  are  neglected. 
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2nRcL  *  n(Rj  +  r2)  #r  -  RX)2  +  L2  (37)  ! 

where:  Rc  =  Radius  of  a  cylinder  j 

L  =  Vertical  height 

Rj  =  Sma1Te,‘  radius  (Fig.  73)  t 

?>2  =  Larger  rad  jus  fig.  73) 

Dividing  both  sides  of  Equation  37  by  ;>JL  yields  Equation  38. 

r  .  "i  »  «2  .  V(R;  -  R1)2  *  i*  (38, 

C  2  L 


Equation  38  indicates  that  a  cone  may  be  approximated  by  a  cylinder  of  radius 

RV*  if  the  square  root  term  divided  by  L  (in  Eq.  39)  is  approximately 
«. 

equal  to  one.  Substituting  the  upper  and  lower  thigh  radii  (Table  25)  for  R2 
and  Rj,  respectively,  into  Equation  39  yields  a  value  of  cos  a  equal 
to  0.9964. 


where:  cos  o  *  Cosine  of  angle  a  in  Figure  73 
R2  •  10.15  CH  (Table  25) 

R,  -  6.87  CH  (Table  25) 

t  «  38.3  CH  (Table  25) 


The  value  of  0.9964  from  Equation  39  corresponds  to  an  angle  a  of 
approximately  4.86*.  Since  this  angle  is  small,  it  is  realisable  to  assume 
that  the  thigh  may  be  treated  as  a  cylinder  without  introducing  any 
significant  error  in  our  predictive  model.  The  thigh  is  orcfeebly  the  region 
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that  Is  least  like  a  cylinder,  therefore,  we  also  assume  that  all  other 
possible  conical  regions  may  be  treated  as  cylinders  without  causing 
significant  errors  In  our  predictive  model.  Since  pressure  is  considered  to 
act  uniformly  In  all  directions  at  a  particular  depth,  we  can  assume  that  the 
pressure  would  act  on  the  inclined  surface  of  a  cone  normal  to  the  unit  area. 
Thus,  the  hydrostatic  head  at  a  given  depth  below  the  shoulder  crest  would 
result  In  the  same  net  compression  on  the  Inclined  surface  of  the  cone  as  the 
cylinder.  This  indicates  that  the  ratios  of  the  radii  would  be  the  same  with 
either  geometry.  From  this  additional  consideration,  we  conclude  that  use  of 
conical  geometries  for  certain  regions  would  be  of  little  merit  in  increasing 
the  accuracy  of  our  model  as  it  is  now  defined. 


